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Abstract 

The properties of the Volume operator in Loop Quantum Gravity, as constructed by Ashtekar and Lewandowski, 
are analyzed for the first time at generic vertices of valence greater than four. The present analysis benefits from 
the general simplified formula for matrix elements of the Volume operator derived in [24], making it feasible to 
implement it on a computer as a matrix which is then diagonalized numerically. The resulting eigenvalues serve 
as a database to investigate the spectral properties of the volume operator. Analytical results on the spectrum at 
4-valent vertices are included. This is a companion paper to [25], providing details of the analysis presented there. 
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1 Introduction 

During the last 15 years major achievements have been obtained towards constructing a quantum theory of gravity in 
a background independent way. Within the framework of Loop Quantum Gravity (LQG) [U [21 [3] it has been possible 
to reformulate General Relativity as a constrainecJl gauge field theory which is then quantized using a refined versio 
of canonical quantization which is based on [5]. In order to complete this program several conceptual and technical 
problems have to be solved. One major task is to understand and evaluate the action of the operators in the quantum 
theory which correspond to the classical constraints. Following the program of Refined Algebraic Quantization these 
operators have to be imposed on the up to now kinematical theory. A key ingredient [71 [8] for this is the understanding 
of the quantum operator corresponding to the classical volume of a region in three dimensional Riemannian space. 

Note that there are different versions of this operator due to Ashtekar and Lewandowski 16J and due to Rovelli and 
Smolin [13j resulting from two different rcgularization schemes. We will focus here on the first version which considers 
diffeomorphism invariant properties of the graphs, namely the relative orientation of the graph edges underlying the 
kinematical basis states encoded in sign factors, whereas the latter is insensitive to this information. Moreover the 
latter rcgularization [T3] seems to lead to difficulties when one attempts to reconstruct the action of the electric flux 
operators from the volume operator [221 123j . In what follows we will always imply the operator due to Ashtekar and 
Lewandowski when using the phrase 'volume operator'. This paper contains the details on the results presented in the 
companion paper [25j , on the properties of the volume operator in Loop Quantum Gravity. We recommend the reader 
to have a look at f25] first in order to get an overview before getting into the details. Note that the basic notation 
which is used throughout this paper is introduced in [25 . Also an introduction to details of Loop Quantum Gravity 
essential for our analysis is provided in [25j and will not be reproduced here. 



^ Upon rewriting General Relativity in Hamilton's formalism, the background independence of the underlying theory is manifested in the 
occurrence of constraints which must be fulfilled by physical configurations in the thus constructed model. By 'background independence' 
we refer to independence from a choice of fixed background geometry. 

^This quantization program is referred to as Refined Algebraic Quantization (RAQ), see [T, for a pedagogical introduction. 
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This paper is organized as follows. In section 2 we recall the construction (regularization scheme) of the volume 
operator from the classical volume expression of a spatial region in order to clarify the origin of the combinatorial 
sign factor and the recoupling in the resulting operator expression. Here we closely follow 15 . Starting from |24| 
we then present in section 3 a more elaborate discussion of the computation of the matrix elements of the volume 
operator, in particular all cases of special arguments of the general matrix element formula obtained in [24] are worked 
out in detail. Subsequently in section 4 we show how considerations of gauge invariance impose conditions on the 
combinatorial sign factors, which reveals a kind of 'self-regulating' property of this operator with respect to high 
valent vertices, providing a suppression of the volume contribution of certain high valent vertices. Next we present a 
detailed analysis on the set of diffeomorphism inequivalent embeddings which exist for the edges of an A^- valent vertex 
into a spatial three dimensional Riemannian manifold. This analysis is accompanied by a numerical Monte Carlo 
sprinkling computation. Section 5 then describes the technical details of the computational implementation. In the 
following section 6 the results of the detailed numerical analysis of the volume spectrum is presented. The eigenvalues 
are computed via singular value decomposition, utilizing the LAPACK library [35]. The numerical work employs the 
Cactus computational framework |36J, which provides automatic parallelism, facilitates use of the LAPACK library 
(including managing inter-language procedure calls), and places the code in a modular context from which it can 
easily be used by others for future investigations. Finally in section 7 we present new analytic results on properties 
of the volume operator at four-valent gauge invariant vertices. We show that its spectrum is simple and give an 
expression of the volume eigenstates in terms of its matrix elements and its eigenvalues. Moreover we show how one 
can (under mild assumptions on the four spins at the vertex) obtain an analytic lower bound of the smallest non-zero 
eigenvalue, solving an outstanding problem of 24J. Our results are then briefly summarized in section 8; for a detailed 
summary we refer the reader to the companion paper [25j . An appendix contains background material regarding spin 
networks and angular momentum theory, in order to make the present series of papers self contained and accessible 
to non-specialists. 

2 Definition and Derivation of the Volume Operator 

Let us briefly introduce our notation first. Note that a basic introduction to LQG is given in the companion paper 
[25] and will not be reproduced here. For more details the reader is referred to [HEllS]. 

As General Relativity is treated in Hamiltonian formalism as an S'f7(2)-gauge field theory, four dimensional space 
time is foliated into three dimensional spatial slices S with induced metric qab{x) and orthogonal foliation direction 
parametrized by a real foliation parameter. The canonical variables are then densitized triads Ef{x) and a connection 
A^(t/), at a;,y € S. Usually i,j = 1,2,3 are su(2)-directions, and a,b = 1,2,3 are spatial tensor indices. The set of 
classical (smooth, continuous) connection configurations is denoted by A, the set of distributional connections is called 
A. The Poisson brackets of the Ef (x) and the connection Aj^ (y) are then given by 



Here k = SttGn, Gn being Newton's constant. Electric fluxes Ei{S) are constructed by integration of the dual electric 
field over two dimensional orientable surfaces S* C S and the integral of the connection along one dimensional oriented 
piecewise analytic edges e C E giving the holonomy he{A) for each A ^ A. In this way (|2.ip is regularized. A collection 
of edges is called a graph 7. Its set of edges is called E{'j). It is constructed such that two edges mutually intersect at 
most in their beginning/end point, called a vertex v in the vertex set ^^(7). Continuous maps f = o are called 
cylindrical functions, the set of cyhndrical functions is denoted by CYL. Here : A {/ie(^)}eeB(7) G SU{2)^^^'^^^ 
and : {/ie(^)}eG_E(7) G 5'?7(2)l^'^^)l (D. Sometimes we will also refer to as cylindrical functions. Moreover 
one finds that spin network functions (SNF) r^j^^(A) — nee-E(7) V^Je + 1 [^Tje ('*e(^))] ^ „ serve as a basis for CYL. 
Here [ttj^ {he{A))^ ^ ^ is a representation matrix element function of an irreducible SU (2) representation of weight je 
which is associated to every e € E{'-f), and TOg, He — —je, ■ ■ - je denote the matrix element. 

2.1 Classical Starting Point 

In this section we will closely follow the construction presented in [15] . 

Let us review the most important steps in order to write down a well defined operator V acting on the kinematical 
Hilbert space of LQG. The classical expression for the volume of an open, connected three dimensional spatial region 
R given by: 



where we choose the Riemannian signature det(g) > Vx and have used the fact that the square root of the determinant 
of the spatial metric qab{x) can be re-expressed in terms of the densitized triads (a;), as introduced in [5S]. 



{Ef{x),Ai{y)]^K5t5l5{x,y) 



(2.1) 




(2.2) 
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2.2 Regularization Scheme 

The regularization procedure displayed here will be similar to the regularization of (|2.ip by smearing the dual of 
electric fields Ef{x) over two dimensional surfaces S and the connection along one dimensional edges e in order 

to obtain its holonomy he{A), however the idea is here to define the smearing surfaces somewhat intrinsically, adapted 
to a graph 7 being the support for a cylindrical function 

Let us introduce the characteristic function XA {p, x) of a cube in a coordinate frame x centered at a point p and 
spanned by the three vectors — Aifii, the fii being normal vectors in the frame x. The cube A has coordinate 
volume vol{A) = 2'^AiA2A3 det(ni, ri2, rJs). Then xa{p,x) can be defined as 

3 

Xa{p,x) = 1[q(^X-\ {n\{x-py) I) (2.3) 
1=1 

and ( • I • ) is the Euclidean inner product, Q{z) is the usual unit step function with Q{z) = if z < 0, Q{z) = ^ if 
z = and e{z) = 1 if z > 0. 





Figure 1 : The regulator XA [p, x) 
If we take the limit A ^ (by decreasing A^ ^ for i = 1, 2, 3) we realize that 



lim 



XA 



{p, x) = S'^^^ (p, x) 



A~*a vol(A) 

where S^^'' indicates the delta distribution in three dimensions. Now upon introducing the smeared quantity 



(2.4) 



one notices that 



vol(A) vol(A') vol(A")y ^'^/^'^/^'"^^^^'^^ XA'ip,^^) XA"{P,^^^) 

R Ft R 

1 



x-e^^'eatcEfix)E'^{y)E^,{z) 



(2.5) 



(i) Taking the limits A^ , A,^ , A" ^ in any combination and at any rate with respect to each other we come back 
to ^e'^''eabcEf{p)E]{p)El{p) due to jMl). 



(ii) (i) holds for any choice of linearly independent normal vectors ni,n[,n'/. 
(iii) 

V(R) = lim hm lim / d^p./|£;(», A, A', A")| 

^ ' A^0A'^0A"^0 7 J-yi v-i-' ' n 



(2.6) 



The classical Poisson bracket (|2.ip contains functional derivatives with respect to the electric field E!^{x) and the 
connection Al{y). Upon quantization the electric field Ef{x) acts as a functional derivative 



El{x) — > E^{x) -ihK 



SAlix) 



(2.7) 
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and we have (symbolically), for the action of this functional derivative on a holonomy hi,{A) of an edge e with x = e(tx) 
being a point of e, (here = —iai denotes a basis for the Lie algebra su{2), Gi are the Pauli matrices) : 



(5fae(A) 



lim + lim 


I" 













^e[0,l] 



(e(i), a;)e°(t)/ie[o,t] i7/ie[t,i] 



a; = e(i)| 

e°(ti:) ^e[o,t]Y^e[ta] a; = e(t)| 
i e"(t2:) /ie[04]"l ^ = e(t)| 

Then we can formulate the action of a smeared version of (j2.7p contained in the regulated expression (j2.5p 

1 



lo<t<i 
lt=i 



(2., 



d'a: XA(2:)i?r(x) 



vol(A) 

on a function /-y(/iei , • ■ • , ^ejv) cylindrical with respect to the graph 7, e.g. for x — 6(0|o<t<i 

1 
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(2.9) 



i?f(p,A)/^ = 



vol(A) 



;GB(7) 



^ / XA(p,e(i)) e"(0 tr 



e[0,l] 



(/lei 7 ■ ■ ■ I ^ejv ) 



(2.10) 



where the trace is taken with respect to the defining representation of SU{2) and we have used the chain rule and 
Leibnitz rule in order to evaluate the functional derivative on /-y. 

Now obviously the regulated expression ()2.5p contains three expressions of the form ()2.10p , which have to be applied 
successively if one wants to evaluate the action of p.Sp on a cylindrical function. This involves a careful analysis of 
the action (|2.8p of the individual functional derivatives, in particular if they act on one edge simultaneously, which 
results in expressions where different r-matrices have to be inserted into the holonomy at different values of the curve 
parameter t. 

The action of E{p, A, A', A") on a cylindrical functions then schematically reads 



Eip,A,A',A")f,= 



^ M ■ eabc X 



X I dtj dt' I dt" XA{p,e{t)) XA'ip,'-^^^^) XA»(p /^*^+^'^*2+^"^*"^ ) e{tre'{t'fe"{t'r Oee'e" (t, i', t") 



(2.11) 



Here M = 31 vo1(A) voi(A') voi(A") ^^"^ extended to all triples of (not necessarily distinct) edges (e e'e"). 

By successive application of the three functional derivatives and Leibnitz rule, one obtains expressions Oee'e"{t-,t' ^t") 
which are combinations of traces of holonomies containing inserted r matrices and partial derivatives with respect to 
holonomies of the appropriate edges. 



In order to evaluate the characteristic functions xa{p, x), xa'{p, ^^), XA"{p, ^±|±^) contained in expression (|2.5p 
one considers a coordinate transformation in the integrals: for a given triple (e, e', e") of (not necessarily) distinct 
edges contained in the edge set E{'j), the vector valued function 



^ee'e" (^i ^ i ^ ) 



eit)+e'it') + e"it") 



(2.12) 



whose Jacobian is given by 



det 



d{t,t',t") 



= dot 



/ 9^1 dx^ dx^ \ 

I dt dt dt \ 

dx^ dx'^ dx^ 

dt' dt' dt' 



33 



eabc eitr e'ity e"(i") 



(2.13) 



\ dx^ dx^ dx^ I 
^ dt" dt" dt" ' 

equals, up to a numerical prefactor, the tangent vector^ contracted with the totally antisymmetric tahc in (|2.5p . which 
is the reason for the choice of arguments in the characteristic function xa" ■ 



'As implied by inserting 112.71 1. 112.81 1 into the classical expression 112.51 1. 
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The observation is now that taking the limit A" 0, and simultaneously demanding that 

XA"ip, '^(*)+'^ ^V''"^ ^ ) = 1, implies that e{tint) = e'(i-„() = e"(i = p, that is the edges e, e', e" have to intersect 
at p at the according curve parameters tint, tint^^int- Therefore the point p is a vertex v in the vertex set V(/y) of the 
graph 7. Moreover due to the convention that edges can at most intersect at their beginning and endpointsp we can 
without loss of generality assume that (e, e', e") are outgoing from p = v, which consequently serves as their beginning 
point and hence tint = t'int = ^'Lt = 0- 

We can thus pull the remaining characteristic functions Xa(p, e(t)) — xa'{p, ^^^^^y-^) — 1 out of the integral. 
The remaining part of the integral then is schematically given as 



E 

e e' e 



^ '''t ""tt"^ uZi'^"^^ de.AOAO) jdtjdt'jdt"Si^\p,.,,^„)dei(^ 

3! VOl(A) VOl(A') eee y , , J J J J y^, eee J \d{t,t' 



Bt"- , ,, 



(2.14) 



In order to get a non vanishing expression from (j2.14p the triple (e,e',e") must consist of three distinct edges with 
linearly independent tangents in p = v. 



Now perform a change of variables 



dH 



det 



d{t,t',t") 



d t sgn 



det 



d{t,t',t" 



dH 



e(ee'e") 



■det 



■ det 



\d{t,t',t") 

B{t,t',t") 



(2.15) 



in order to evaluate the distribution in (j2.14p . Note that it is here where the sign factor e(e e'e") enters. For this 
purpose we have to insert a 1 = [e(ee'e")]^ and, observing that integrating the J^^-* (p, Xee'e") over the positive octant 
only gives an additional prefactor of i, we arrive at 



^lim^£;(p,A,A',A")A 



E 

e e'e"e-B(7) 



i ip 33 Xa{p,v) XA'{p,v) , , n n n^ 

3!8.8 vol(A) vol(AO O Qeev» (0,0,0) 



(2.16) 



where 



(2.17) 



and XI :— tr (ri/ie[o,i] g^^^ ) denotes the right invariant vector fields of SU(2) resulting from the action of the 



functional derivatives (|2.8p . Note that we can choose the order of the right invariant vector fields arbitrarily in (|2.17p . 
because they commute if they act on distinct edges, that is distinct copies of SU(2). 

The second prefactor i in (|2.16p stems from the fact that we only evaluate one sided functional derivatives according 
to ()2.8p : For a vertex v in the vertex set ^(7) of the graph 7, the point p is the beginning point of the outgoing edges 
e,e',e", that is p = e(i)|t=o - e'(i')|t'=o = e"{t")\t"=o. 

Inserting (|2.16p back into (|2.6p we may synchronize the remaining limits A, A' ^ by choosing A = A' and take 

Xa(p, f) — [xa(p, i')]^ out of the square root. Taking the limit A —> then results in the operator describing the 
volume of a spatial region R, namely the volume operator V{R)-y acting on the cylindrical function /-y over a graph 7 



where 



ViR)^f^^ j d'p^det{q){p)^ = j d'pV{p)^ 
R R 



\ 



\i-Z ^ €{ei,ej,eK) €ijkX}X^jX^ 



(2.18) 

(2.19) 
(2.20) 



ej. = .7,e^eE(7) 



^This can always be achieved by subdividing and redirecting the edges of a graph 7, see |25| . 
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Here Z is a constant depending on the regularization procedure. Its numerical value according to [TS] and above is 
found to be Z = ■ {j)^ , however we will keep it unspecified in our calculations, since the Immirzi parameter P can 
be freely chosen. It only contributes an overall constant scaling to the spectrum of the volume operator. Mostly we 
will set Z ~ 1 below and reinsert it when necessary. 

The sum has to be taken over all vertices v € V{'y) of the graph 7 and at each vertex v over all possible triples 
(e/, ej, ex) of outgoing edge^ of the graph 7. Here e{ei, ej, ck) is the sign of the cross product of the three tangent 
vectors of the edges (e/, ej, ck) at the vertex v. 

Again the are the right invariant vector fields on SU{2) fulfilling the commutator relation 

[X}, X^j] = -2 Sjj e'^^ X^ (2.21) 



2.3 Right Invariant Vector Fields as Angular Momentum Operators 

One can introduce the self-adjoint right invariant vector fields Yj := —^X-'j fulfilling the usual angular momentum 
commutation relations [F/,yj] — i Sjj e^^^ Yj . 

Their action on cylindrical functions and hence spin network functions T^j^^-. (A) can be associated to the action 
of ordinary angular momentum operators J| acting on a (recoupled) abstract spin system | a J M ;n ) as analyzed 
in [22!, 23 . This correspondence is summarized in the companion paper [25j and recalled in more detail in section fP. 41 
in the appendix of this paper. We can therefore equivalently replace: 



2X3 



1 



,X}X^jX^ ^ e,jkJ}J^JK \SJM;n) (2.22) 



Using furthermore the antisymmetry of eijk and the fact that [J|, Jj] = whenever I ^ J we can restrict the 
summation in (|2.20p to / < J < A' if we simultaneously write a factor 3! in front of the sum. The result is: 



K,7 = J^-3!-23 e{ei,ej,eK) e^,u J\J'jJ'k (2.23) 

Y I<J<K 

We then make use of the identity 

J^J'jJ'k = \[{JufAJjKf] (2.24) 

where {Juf = ELi('^/ This relat ion can be derived by writing down every commutator as [( J/,/)^, ( Jjk)^] = 
3 

S ~^ Jj)^ ' i"^] + Jk)"^] ' using the identity [a, be] — [a,b\c + b [a, c] for the commutator, the angular momentum 

commutation relations ([J% J-'] = ie^^'^j'^) and the fact that [j], Jj] = whenever I ^ J . 
We may then summarize: 



^vn = J\Z- 4ei,ej,eK) qijK\^- \Z\-Q" ^ \Z\^ -Q^Q" (2.25) 

Y I<J<K 

where qijK '■= \{Ji.jY, {JjkY\ and Z ^ Z ■ j ■ il ■ 2^ = i ■ 0^ ■ j3 again denotes the Immirzi parameter. Note that 



the actual numerical value of Z varies in different regularizations [161, 115pl . However [22] shows how it can be fixed to 
2 — 0i . 'ih . Creg , with Creg = sTg ' demanding consistency of the volume quantization to the usual quantization 
of the electric fluxes. We will set Z = 1 va our subsequent analysis, as it only gives an overall numerical scaling of the 
volume spectrum. 

We have introduced the shorthand Q^Q as a totally symmetric real matrix. Its eigenvalues AgtQ are real and 

come in pairs Agtg = |Aq| . Equation (I2.25P can be understood as follows: The Volume operator V has the same 

I - I - 

eigenstates as Q or Q'^Q, but its eigenvalues are defined as Ay := |Aq| ^ — \^q\q\ " ■ 

Our task is then to calculate the spectra of totally antisymmetric real matrices of the form: 

Q:= <IJK)qijK (2.26) 

/<,/<i<'<Ar 

where e{IJK) = sgn( det (e/(w), ej(w), e^f (w))) denotes the sign of the determinant of the tangents of the three edges 
e/, ej, ck intersecting at the vertex -y, and N is the valence of v. 
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As mentioned before, one can without loss of generality always redirect edges such that there are only outgoing edges at each vertex. 
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exactly reproduces the value of Creg obtained in 
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2.4 Matrix Elements in Terms of 3nj-Symbols 

Now we can apply the recoupling theory of n angular momenta as introduced in the companion paper [25] . Further 
details are provided in appendix section [B] The aim is to represent qijk in a recoupling scheme basis. We will do this 
with respect to the standard basif0 a(12), where we can now easily restrict our calculations to gauge invariant spin 
network states by demanding the total angular momentum J and the total magnetic quantum number M to vanish, 
i.e. we will take into account only recoupling schemes which couple the outgoing spins at the vertex v to resulting 
angular momentum 0. In terms of the recoupling schemes these states are given by: 

I giU) J J = M = ) I g{IJ) ) (2.27) 

where we have introduced an abbreviation, since the quantum numbers j , J — , M = are the same for every 
gauge invariant spin network state with respect to a fixed A^-valent vertex v with edge spins j := (ji, . . . , jat). 

We will now represent qijk ■= [{JuY, (J./k)"^] in the standard recoupling scheme basis where | a ) := | a(12) ), 
I a' ) := I a'(12) ). The point is that by construction a recoupling scheme basis | g{IJ) ) diagonalizes the operator 
{G2f = {Jij? = {Ji + Jj?: 

(G2f\ g{IJ) ) = g2{IJ){g2{IJ) + l)| g{IJ) ) (2-28) 

Furthermore every recoupling scheme | g{IJ) ) can be expanded in terms of the standard basis via its expansion 
coefficients, the 3nj-symbols. So it is possible to express jlS] 

{a{l2)\qijK\(i'{l2)) = 



{d{12)\[{J I j)\{JjKf]\d' {12)) 

{a{l2)\{Jjjf{JjKf]\a' {12)) - {a{\2)\{J jK)\Jij)^\d' {12)) 

g2{IJ){92{IJ) + l)[{d{l2)\g{IJ)){g{IJ)\{JjK)^\d'{l2)) - {d{l2)\{J.jK)^\g{IJ)){g{IJ)\d' {12))] 



9{IJ) 



J2 g2{IJ){g2{IJ) + l)g2{JK){g2{JK) + I) {g{IJ)\g" {12)) {g{JK)\g" {12)) x 

g{IJ),g{JK),g" (12) 

x[{g{IJ)\d{12)){g{JK)\d\l2)) - {d{12){JK)\g){g{IJ)\^^ 



E 

9"(12) 



J2 92{IJ){g2{IJ) + l){g{IJ)\g"{12)){g{IJ)\d{12)) 
aiu) 



X 92{JK){g2{JK) + l){g{JK)\g"{12)){g{JK)\d'{12)) 

g{JK) 



[ a(12) ^ a'(12) 
Y [f^^(5(12),5"(12)) xF-^^(a'(12),g"(12)) 



5(12) a'(12) 



9" (12) 

where we have introduced the shorthand 



i^^-'(a(I2),5"(12)) 92{IJ) [g2{IJ) + l) { g{IJ) \ g"{12) ){ g{IJ) \d) 

g(iJ) 



(2.29) 
(2.30) 

(2.31) 



We can nicely see from ()2.29p that we obtain a real antisymmetric matrix possessing purely imaginary eigenvalues 
(we could alternatively consider the purely complex version by multiplying all matrix elements by the imaginary unit 
i). We have inserted the suitable recoupling schemes | g{IJ) ), | g{JK) ) diagonalizing {Jij)^ and {Jjk)^ and their 
expansion in terms of the standard basis | g{12) ) by using the completeness of the recoupling schemes | g{IJ) ) for 
arbitrary / 7^ J H 

1= Y \9{IJ)){9{IJ)\ 
So we have as a first step expressed the matrix elements of qijK in terms of 3nj-symbols. 



(2.32) 



'^See appendix section for the definition. 

*The summation has to be extended over all possible intermediate recoupling steps §2, . . . , g-n—l, that is [jV -~ 3q\ < 9kUq^Jr) < jq + jr, 
allowed by the Clebsch-Gordan Theorem. 
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2.5 Removal of the Arbitrariness of the Edge Labehng 



Wc have described how one can represent the volume operator as a matrix acting on a linear vector space whose basis 
states are labelled by the multilabel a of the standard recoupling schemes. As can be seen from definition IB. 31 in the 
appendix there exits a unitary basis transformation between any two different orders of recoupling, e.g. g{IJ), a(12): 

I g{IJ) J M ; n) = ^ { a(12) J M ; n \ g{IJ) J M ; n ) ■ \ a(12) J M ; n ) (2.33) 

3(12) 

where the matrix elements ( a(12) J M ; n \ g{IJ) J M ; n ) =: Usi^ of this basis transformation are the so called 
3nj-symbols. Now this basis transformation is unitary as can easily be seen because the recoupling schemes are 
orthonornial 

( a(12) J M ; n \ a'(12) j' M' ; n' ) = Uss' = SaS'Sjj'SMM'Snn' (2.34) 
and provide a complete orthonormal basis 

< a(12) J M ; n\ a'(12) J M ; n) = 

_ .a(12) 
~ "a'{12) 

= { a(12) J M ; n \ g{IJ) J M ; n){ g{IJ) J M ; n \ a (12) J M ; n) 
9{IJ) 

= E U-sUga' (2.35) 
Moreover the 3nj-symbols are real (the overline denotes complex conjugation) and symmetric 



( a(12) J M ; n\ g{IJ) J M ■ n ) ^ Ug^ ^ Us^ = { a(12) J M ; n \ g{IJ) J M ■ n) 

= C^gS = ( g{IJ) J M ] n\ a(12) J M ; n) (2.36) 

Thus we have 

f/C/t = UU'^ = t/^ = 1 (2.37) 

Because the change of the recoupling order can be implemented as a unitary transformation, we have for the matrix 
representation Q of the volume operator of ()2.26p . after a transformation from the standard basis a(12) to a new basis 

Qs'g ( 5' I Q I 5 ) = ( E ^s'a' «' I Q I E ^9-^ ^ ) = E ^la'^rai «' I Q I « ) = E ul,^,UraQa'a (2.38) 

W a a'a a'a 

The matrix Q is being transformed as Q i-^ U~^QU because ~ U~^. Such a unitary transformation does not 
change the spectrum of Q, hence the spectrum of the volume operator is independent of the chosen order of recoupling. 

We use this property to drastically decrease the number of assignments of spins to edges we must use in the 
numerical analysis, because we can always choose a particular edge labelling ei, . . . , e^r in which the spins are sorted 

ji < h < ■ ■ ■ < .In = jmax (2.39) 

Such a labeling then corresponds to 

arbitrarily labeled spin assignments, where Ni is number of elements in one of p sets of mutually identical spins, 
^i^i Ni — N. Throughout the remainder of the paper we refer to such a sorted assignment of spins to vertex edges 
as a spin configuration (or j-configuration). 



3 Explicit Matrix Elements of the Volume Operator 
3.1 Starting Point for Matrix Element Implementation 

The simplified closed expression for (|2.29p . which has been derived in [23], serves as a starting point for our compu- 
tations. Its derivation will not be reproduced here. However, in order to make it implementable on the computer, we 
will recall the crucial intermediate results of [Hj and give the detailed outcome of the main formula for the matrix 
elements for all combinations of its arguments. 
Let us first state the complete results on 

i^^'^(a(12),5"(12)) E .92 (/J) (.92 (/J) + l){g{IJ) \ g"{l2) )( g{IJ) \ a(12) ) (3.1) 
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as contained in (|2.29p . 



1=1 J ^2 



N 



(a2 + 1) n ^9>. 



1=1 J = 3 



1=1 J > 3 



1=2 J=3 



1 = 2 J > 3 



{'I 


31 


f } ("1)''^ < 


' 0.3 


j3 






02 


I 1 


02 


js J 



§2 <12 



AT 



J2 Jl 52 
02 jl 



J-1 



n^(5>«)(-i) 



n— 3 



1 a„ a„_i 



X (-1)" 



aj J OJ 

1 Oj_.i j,7 



J-1 



ft=2 



k=J 



' N 
. fc=3 



31 32 92 

0.2 32 



/ iNa, J 03 J3 .92 
^ ^ W a2 J3 



Jl J2 52 
02 32 



n A(.g>„)(-l)-^"+^"-+--+i ^'f-i /" j X (-1)- \ Y 5}-! 

^\ I i On ttn-l I 11 aj-i J J 



n— 3 



/ > 2 J = / + 1 



j-i 



/-I 
-2 E J„ 



k=J 



3i+i)^-A{g'{ ai) 



oi-i 31 9 1 
1 o/ j/ 
/ 

+c(j/ 3i+i) n '^9k«'= 

fc=2 J 1=2 



X (-l)'^-f + i 

/-I 



0/+1 31+1 9i 

1 O/ 

N 



k=I+l 



I>2 J > I + l 



I-l 
-2 E J„ 



.7-1 



i(-i)^ "(-1) "-+K-i)^'^'-^'^(-i)°^-^+H-i)'^^-^"^0/jj)^^(5y o/) 



(3.2) 



oi-i 31 9i 
oi ji 



J-I 



n ^(5>n)(-l) 



g„_i+an-l+l J 3n 5n-l 5n 



n=/+l 



1 On 0„_i 



X (-1)" 



a., J Oj J J 5j_] 
1 Oj_i jj 



J-1 



fe=2 



/-I AT 

n ^9i'-i n 

;=2 fe=,/ 



Where we have used the shorthands 

C(a, b) = a{a + 1) + 6(6 + 1) X{a. h) = 2a{2a + l)(2a + 2)26(26 + 1)(26 + 2) 

(3.3) 

A{a, 6) = v/(2o+ 1)(26+ 1) 
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The explicit expressions F^'^(a(12), 5"(12)) of serve as the building blocks in order to calculate the matrix 
elements ( a(12) | qijK \ a'(12) ) given in (|2.30p as: 



( a(12) I qijK I S'{12) ) = J2 

9" (12) 



F"{a{12),g"{12)) x F'^^ (a'(12), r'(12)) 



a(12) ^ a'(12) 



One can evaluate (|2.30p as written down above by using the appropriate expressions from (|3.2p . However there is still 
some simplification possible [24 . 

When looking at (j3.2p . obviously each expression comes as a sum of two terms. So multiplying two expressions 
F-''^ {■) X F"^^ {■) as in the first term on the right hand side of (|2.30p would result in four terms. We will now show that 
by symmetry all product terms containing a symmetric C(-) factor vanish, as we antisymmetrize with respect to the 
interchange [a(12) ^ a'(12)]. For every term containing C(-) we effectively get a prefactor of nA;=2 ^a'^a,, or nfe=2 ^9'k<^'k 
if we appropriately concatenate the products of (5-factors in p.2p . Therefore the product terms C{ji,jj) x C{jj,jK) 
are obviously symmetric under [a(12) ^ a'(12)] and vanish under antisymmetrization. 

Let us discuss the J-factors in (|3.2p in order to exclude mixed terms containing one C(-). We will present 
here only the discussion for products containing C(jj,jx) in F''^ (a'(12), g"(12)) . The discussion of C{ji,jj) in 
F^'^(a(12), g"(12)) is completely analogous. 

For every term containing C{jj,jK) we effectively get a prefactor of 11^=2 ^g',ia-'^ concatenate the products of 5- 
factors. Therefore when considering a product of the form (|2.30p . with a non-C(://, jj)-term contributed by F^'' {■) and 



^ C'Ojj j-fc)-term coming from F'^^(-), every g'l in the F^'^(-)-term can be replaced by a'^. That is (f"(12) = a'(12) 



But then the mixed terms containing one C{jj,jK) become symmetric wrt. [a(12) ^ a'(12)] for the following reasons: 
• (-l)"2-«2 = (-i)4-'^2 because (4 - 02) £ Z 



• A(a'f,,ak) is symmetric by construction 



all (but 2) 6j-symbols occurring in definition (|3.2p for F^'^(-) are of the form 

f jl 9n-l 9n ! = an-l 1 

\ 1 a„ a„_i J \ 1 5;' J 

by the symmetry propertiecl of the 6j-symbols. But .g"(12) = a'(12). 



I-l N 



• due to the Y[ ^g'^au W ^g',ia„ in the non-C(j/, jj) terms of F^'^{-): 

k=2 n=J 

the remaining 6j-symbols | '^^^^ -^^^ j, g^j^j _f jj 9j~i then symmetric 

wrt. a(12) ^ a'(12), because g"{12) = a'(12) 

and hence, gj_i = a'j_^ , a/_i = g'j'_^ = a'j_^ , g'j' = a'j, aj ^ g'} ^ a'j . 
Thus when calculating (|2.30p using (j3.2p we are always left with the single product not containing a symmetric C{ji,jj) 



or C(jj, j/s-)-part. Note, however, that for cases containing / = 1 J = 2 in (|3.2p we get a contribution 



Let us finally take a closer look at the (5-products in the remaining non-symmetric product terms of (|2.30p . In 
these products we have for I < J < K: 



7-1 N 

li^g'^-^Il^g^i 

,fc=2 l=J 



J-1 N 

n '^9;>™ n ^9>'n 
,m— 2 71— K 



I-l J-1 K-1 N 

n ^->^< n n n ^---'^^ (^.4) 

k=2 m=I l = J n=K 



We have thus to look at the symmetries between / ... if — 1 only. Here we find: 

J-1 K-l 

For the first term in (^30]) 5(12) a'(12) J| ^giL^^'^ JI ^al'^-i 



m=I l=J 
J-1 K-1 



(3.5) 



and the second term in (|2.30p a'(12) a(12) ^s^a™ Y\. ^g'l'^'i 

( a(12) ^ a'(12), a'(12) ^ a(12) ) m=/ " i=j ' ' 



Using again symmetry properties of the 6:/-symbols, and integer-arguments for the (— l)-exponents, together with ()3.5p . 
we can show, as in the case of symmetric terms, that almost all terms are symmetric with respect to the interchange 
a(12) ^ a'(12). 

So we are left with the discussion of the following terms: One easily sees that for a product F^'^(-) x F'^^{-) (of 
two terms from p.2p using the according arguments given in (|2.30p ) one has: 



^6j-symbols are invariant under the interchange of 2 columns and the simultaneous flip of two columns, see appendix IC. 31 
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/-,\ai-i(-,\ai-aj/^\aj\ '^J J J ^,7-1 



^ ^ ^ ^ 11 JJ a'j 



(3.6) 



Here a/_i = by P-Sp . (a/ — a^) G and the aj exponents cancel. Thus the remaining part of the product terms 
in (|2.30p not symmetric wrt. the interchange a(12) a' (12), a' (12) — > a(12) is given by 



aj Jj aj_^ 
1 a J -I jj 



1 jj a'j 



We win give the exphcit matrix element expressions ()2.30|) in the following section. 



(3.7) 



3.2 Explicit Matrix Element Expressions 

The discussion above leads us to consider the following 14 special cases for the matrix element ( a(12) | qijk \ a'(12) ): 



1=1 J = 2 


K = 


3 




A > 


q 
O 


1=1 J = 3 


K = 


J + 




K > 


J + 


1 J > 3 


K = 


J + 




K > 


J + 


J = 2 J = 3 


K = 


4 




K > 


4 


/ = 2 J > 3 


K = 


J + 




K > 


J + 


J>2 J>H 


-1 K = 


J + 




K > 


J + 


J>2 J>H 


-1 K = 


J + 




K > 


J + 



(3.8) 



Using combinations of the according cases (|3.2p as written down in (|2.30p one obtainj^: 
3.2.1 



1=1 J = 2 K = 3 



{ d{12) I 9123 I a'(12) ) 



^ [F''{a(12),g"{12))F''{a'(12),g"{12)) 



9"(12) 



a(12) ^ a'(12) 



1 1 , , 1 (3-9) 

a2(a2 + 1) - 4(4 + 1)J -(-l)--''^-^M-lP+'(-l)"^""H-l)"^'^(j2 J3)^A(4 a2) x 

N 





32 


0,2 1 


f 03 h 02 1 




a'2 


J2 J 


I 1 a'2 J3 J 



71 — 3 



^''For completeness all Kronecker-(5's at the end of the following equations range from K . . . N, however if we consider gauge invariant 
states, besides = ajv = J = 0, we automatically have a^_i = ijv— 1 ^ ''^ every case. 
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3.2.2 1=1 J = 2 K>2, 



( a(12) I qi2K I a'(12) ) 



^ [Fi2(a(12),/'(12))F2^(a'(12),5"(12)) 

S"(12) 



a(12) ^ a'(12) 



\ 1 a; a:,_i / ^ \ 1 a'^_, jk ] 



a2{a2 + 1) - a^(4 + 1)] ^ (-l)-^-^^-!)^'"^^-!)"^-''^ ^(j2 jx)^ A(a2 4) x | 
X II A{an a'J (-l)-J"+«"-i+«n-i+i 

n=3 
AT 



X I I an 



(3.10) 



3.2.3 7=1 J = 3 K = J+1 



( a(12) I gi34 I a'(12) ) 



[F'^{a{12),g"{12))F^^{a'{12),g"{12)) 



a" (12) 



a(12) a'(12) 



(_1)J1+J2(_l)j4(_l)a^ X^j^ J3)2X{J, J4)^ A(4 a2) A{ 



^ ^ I 1 02 J t 1 ^3 J4 J 



(-l)H-i)3<j 1 1 ^.^ 



(-l)°2(-_l)a3 



f "3 .?3 


0.2 1 


J a.2 


is 


41' 


I 1 "2 


ja J 


I 1 


as 


is / 



n=4 



(3.11) 



3.2.4 



7=1 J = 3 7f> J+1 



( 0(12) I q^sK I 5'(12) ) 



^ [Fi3(a(12),5"(12))F3^(a'(12),5"(12)) 



9" (12) 



5(12) <^ a'(12) 



K-l 

1 (_i)aW x{j, js)ix{js jK)i Aia'2 02) ^(o^ as) 



K-l 



X < 1 z I X n ^c- oi-D-— <-« »;,-' 



1 a'^_i Jk 



X 



r 1^a2^^ 1 J "'^ «2 I I ^2 h as 
1 a2j3/l 1 4J3 



(-1)"2(-1)«3 j "3 .73 0-2 ^ j 0-2 js as 



1 "2 is J 1 1 as is 



X n "^"n «n 



(3.12) 
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3.2.5 1=1 J>3 K = J+1 



( a(12) I gijj+i I a'(12) ) 



^ [Fi^(a(12),5"(12))F^^+i(a'(12),5"(12)) 

S"(12) 



a(12) ^ a'(12) 



j-i 



J2 Ji 

1 02 jl 



J-1 



n— 3 



a^_i+a„_i + l J 3n O'n-i '^n 

1 o„ a„_i 



a,7+i Jj+i a,/ 
1 a J JJ+i 



jj 



' ' ^ ' \1 aj., jj j\ 1 a'j jjf ^ ^ \ 1 a'j_, 

N 



aj-1 Jj aj 
1 aj 3 J 



n=J+l 



(3.13) 



3.2.6 1=1 J>3 K>J + l 



( a(12) \qxjK\ a' {12) ) = 



[F^\a{12)J'{12))F'^{a'{12),g"{12)) 

S"(12) 



a(12) ^ a'(12) 



K-1 



(_l)n=i n=.+i " (_i)iK(_i)< x{n jj)'^X{jjjKf^A{a'^ a2)A{aj a'j) 



J2 Jl 02 

1 a2 jl 



.7-1 



X n A(a;a„)(-1) 



n=J+l 



a^_i+o„_i + l J ''■n-l 

1 a„ dn-l 



Ik Jk ax-i 
1 cik-i Jk 



1 / aj-i jj a'j \ 
/ I 1 aj jj] 



N 



n=is: 



(3.14) 



3.2.7 1 = 2 J = 3 K = A 



< 5(12) I 9234 I a'{12) ) 



^ [F^\a{12)J'{12))F^\a'{12)J'{12)) 



S"(12) 



a(12) ^ a'(12) 



(_l)..+..+.4(_i)a;^(_^.^ j3)^^(j3 J4)^A(4 «2)A(a3 4) ^3 f 
4 I i ^2 J2 



/I f ^/ 



(-1)"^-!)''^ \ 'J f 

J- ^2 J3 



^3 .73 a2 

1 a2 J3 



JV 



n=4 



0-4 .74 aa 

02 j3 03 

1 03 h 



(3.15) 
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3.2.8 1 = 2 J = 3 K>A 



( a(12) I g23K I a'(12) ) 



^ [F^\d{l2)J'{l2))F^^{d'{12)J'{l2)) 



S"(12) 



o(12) <J4> a'(12) 



X ^ "2 ^ ^ 

1 02 J2 



n=4 



a„_i+o^_i + l J Jn dn-l In 



1 O'n-l 



1 a'^-i Jif 



' 1 I „' 



(_l)a2(_i)a3 , 03 J3 a2 ^2 J3 Og _ I a 

^ ^ I 1 02 J3 j i 1 as J3 ' ^ ^ ^ ^ 



^3 i3 02 1 r a2 j3 a 
1 02 i3 j 1 1 03 i: 



n=K 



3.2.9 7 = 2 J>3 7^ = 7+1 



( o(12) I q^jj+i I a'(12) ) 



^ [F2^(a(12),5"(12))J^^^+i(a'(12),5"(12)) 



S"(12) 



a(12) <^ a'(12) 



(3.16) 



'I z t } ^ n/("» <)(-i)--"-« { ^ } x { »■'« '7, 



Jj+1 aj 

j.J+1 



N 

X n ^<«n 



3.2.10 



/ = 2 j>3 is:>j + i 



( 0(12) I 92 7 7+1 I 0'(12) ) 



^ [F2-^(a(12),5"(12))F-^^+i(a'(12),5"(12)) 



S"(12) 



07-1 J7 Oj 
1 Oj ij 



a(12) <^ a'(12) 



(3.17) 



X < ^2 "2 ^ X 

1 02 J2 



{-lY^-'^{-iy-{-l)<X{n jj)-^X{jj jKV^A{a'^ a2)A{aj a'j) 
n A{an <)(-l)«''-+<-+i I y f" 



7-1 



n—3 
K-1 



X n A(a„ o;)(-i)""-+<-+N -^r ''"r' >x 



n=J+l 

^ ^ ^ ^ 11 Oj_i jj j\ 1 O'^ 

AT 



O^ Jif Oif-l 

1 a'^_i jx 



Jj 



1 < <-i 

^ ^ ^> \ 1 a'j_, jj f\ 1 aj 3 J ] 



(3.18) 
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3.2.11 J>2 J = I+1 K = J+1 



( a(12) I qii+ii+2 I a (12) ) 



J2 [F''+'{a{12),9"{12))F'+''+\a'il2),g"{12)) 



Ls"(12) 



a(12) ^ a'(12) 



h-iy'+j'+^{-l)-'-^{-lfr+2X{jjjj+,)ix{jj+^ ji+2V^A{a'j ai)A{ai+, a'j^,) 



( ai-i ji «/ 1 X f a/+2 J/+2 



ai+i 

jl+2 



/-I 



n'^<«'. n ^< 

n=2 n=I+2 



ai JI+1 ai+1 
1 <+i jj+l 



(3.19) 



3.2.12 



I>2 J = I+1 K> J+1 



( a(12) I qii+iK I 5'(12) ) 



5^ [F^^+i(a(12),5"(12))F^+i^(a'(12),9"(12)) 



Ls"(12) 



a(12) ^ a'{12) 



K-1 

E j 



i-iy'i-l) ''(_i)«.-i(_i)^--(_i)aWx(j, i,+i)^^(i/+i jKf^A{a'j ai)A{ai+i a'j+,) 



1 0/ J/ 



(-l)"^(-l)"^+i 



A'-l 



n=7+2 n n 1 

a/+i 



aK-1 

1 o'if_i 



^+1 ji+1 a'j 1 f < ji 
1 a/ j \ 1 a'l 



aj+i 

7+1 



/ 1^aW iNaj + i I Oj+l aj 1 I O/ a/+i 



7-1 



TV 



X n ^< n 



n=2 



(3.20) 



3.2.13 



I>2 J>I+1 K = J+1 



( a(12) I qijj+i I a'(12) > 



[f' ' {S{12),g"{l2))FJ (a'(12), 5"(12)) 



S"(12) 



a(12) <^ a'(12) 



J~l 

1 (_!)„?/" (_!)«.-! (_l)a.-ai^_l),.+,^_l)a',+,^(^-^ j7)^X(jy j7+i)^A(a^ a,)A(aj a'j) 



X < "^-1 > X 

1 0/ J/ 



J-1 



n A{a'^an){-ir 



n=/+l 



a„_i+aj,_i + l J Jn O-n—l On 
1 < «n-l 



1 a'j jj+i 



^ ^ ^ 1^ 1 1 aj., Jj j\ 1 a'j jjj^^> ^ ^> \1 a'j_, jj j [ 1 aj jj 



7-1 



N 



X n "^"n «n n ^< 



n=2 



n=J+l 



(3.21) 
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3.2.14 I>2 J>I+1 K>J+1 



{ a(12) I qijK I a (12) ) = 



5] [F"(a(12),/(12))F^^^(a'(12),/(12)) 



Lg"(12) 



a(12) ^ a (12) 



i(-l)'S^-'"(-l)""-^+i-'"(-l)'^-n-l)''^""^(-l)^"(-l)''-^(j/ JJ)^^(JJ Ji^)^^(«/ a/)^(aj a'j) 



^ n— i+1 ^ 
K-l 



n— J4-1 



7-1 



N 



n=2 n=K 



(3.22) 



3.3 The Arguments: Special 6j-Symbols 

In all cases discussed, very special 6j'-symbols occur: All of them contain a 1 as an argument. 6j-symbols of this kind 
possess short closed expressions and they restrict the arguments allowed (see [29]). 

3.3.1 Different Cases and their Classification 

The following four types of 6j-symbols occur in (I3.9p .. . .■ P.22p : 



II[b,c] 



III[b, c] 



a b c 
1 c- 1 fe- 1 



a b c 
1 c-1 b 



a b c 
1 c- 1 6+ 1 



= (-1) 



(-1) 



(-1) 



a+b+c 



{a + b + 


c){a + b + c + 


!)(-« + 


6 + c - 1) (-a + 6- 






(26- 1)25(25 


+ l)(2c- 


- l)2c(2c+ 1) 




2{a + b 


+ c+l){-a + 


b + c){a 


-6 + c)(a + 5- c-(- 


-1)' 




2b {2b +1) (26 


+ 2) (2c- 


- l)2c(2c+ 1) 




(a - 6 + 


c - 1) (a - 6 + 


c) (a + 6 


-c+l){a + b~c- 


f2)' 



(26 + 1) (26 + 2) (26 + 3) (2c - 1) 2c (2c + 1) 



jy[6,c] 



a b c 
1 c 6 



= (-1) 



a+b+c+l 



6(6+l) + c(c+l)-a(a+l) 



26 (26 + 1) (26 + 2) 2c (2c + 1) (2c + 2) 



So we always have Gj-symbols of the general form 



a p q 
Irs 



This leads us to consider nine possible cases of argument combinations 

p — s 



(3.23) 



(3.24) 



+ 1 

-1 





q = r + I 


q = r 


q = r — 1 


p — s + 1 


(1) 


(2) 


(3) 


p — s 


(4) 


(5) 


(6) 


p = s — 1 


(7) 


(8) 


(9) 



(3.25) 



Using the symmetry properties (|C.4p . (|C.5p of the 6j-symbols, these 9 combinations are covered by the 4 types of 
6j-symbols contained in p.23p fwith the appropriate arguments in brackets): 
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(1) 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 



a p q 
q q-1 p-1 



a p q 
1 q p-1 



a p q 
1 q+1 p-1 



a p q 
1 q-1 p 



a p q 
1 q p 



a p q 
1 q+1 p 



a p q 
1 q-1 p+1 



a p q 
1 q p+1 



a p q 
1 q+1 p+1 



a q p 
1 p-1 q 



a q p 
1 p-1 q+1 



a q+1 p 
1 p q 



a q p+1 
1 p q 



a p q+1 
1 q p 



a q + 1 p + 1 
1 p q 



I[b ^p^c = q] 
II[b g,c = p] 
III[b^q,c^p] 
II[b = p,c=q] 
IV[b = p,c=q] 
II[b = p,c = q+ 1] 
III[b^p,c^q] 

I[b = q + l,c = p+l] 



(3.26) 



3.3.2 Inequalities to be Fulfilled for the Arguments of I, II, III, IV 

Following the general definition (|C.ip of the 6j-symbols certain inequalities for their arguments must hold. Due to 
the special form of the arguments these 12 inequalities reduce to a set of 5 inequalities for each case: 



I[b, 



a b c 
1 c-lb-1 



Aia,b,c) : 

A{a,c-l,b- 1) 
A(l,6,6-1): 
A(l,c-l,c) : 



a+b-c> 
a-b+c> 
-a + c + b > 2 
b > 1 
c > 1 



II[b,c] 



a b c 
1 c-lb 



A{a,b,c) : 
Aia,c-l,b) 

A{l,b,b) : 
A(l,c-l,c) 



a+b-c> 

a-b + c > 1 

-a + b + c> 1 

c > 1 



III[b, c] 



a b c 
1 c- 1 6+ 1 



A{a,b,c) : 

A{a,c- 1,6+1) 
A{l,b,b+1): 
A{l,c-l,c) : 



a+b-c> 
-a + b + c > 

a-b+c > 2 
6 > 
c > 1 



IV[b,c] = 



a b c 
1 c b 



A{a,b, c) 



A{a,c, b) 
A(l,6,6) 
A(l,c,c) 



a+b-c> 
a-b+c > 
-a + b + c> 



b > \ 
c > i 



Moreover we have an integer condition, namely in each case a + 5 + c G ]N must be fulfilled. All other 6j-symbols 



containing a 1 but different argument ranges as given above or not fulfilling the integer condition vanish identically. 
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4 Mathematical Preparations for Implementation on a Computer 



Here we address the remaining mathematical issues necessary for numerical computation of the volume spectrum. 
Recall from (|2.25p that we have to calculate the spectra of totally antisymmetric real matrices of the form: 

Q:= J2 <IJK)qijK (4.1) 

I<.]<K<N 

where again e{IJK) = sgn( det {ei{v),ej{v),eK{v))) denotes the sign of the determinant of the tangents of the three 
edges ei,ej,eK intersecting at the vertex v. The actual computation is twofold: firstly we have to compute the 
recoupling scheme basis in order to apply the matrix element expressions of section 13.21 for the computation of the 
constituent matrices Qijk- Secondly the possible sign factors e{IJK) have to be computed. 

4.1 Recoupling Scheme Basis 

Let us first describe the computation of the recoupling basis. As it turns out, the number of gauge invariant recoupling 
schemes contained in all possible recouplings of a set of spins can be surprisingly small. To our knowledge an exact 
result for this number is known only for the special case where all spins equal j — \ |21j . 

By definition IB. 21 in the appendix a standard recoupling scheme is the successive recoupling of all N angular 
momenta ji , . . . , at the vertex v: Fix a labelling of the edges. Then recouple ji , j2 to a resulting angular momentum 
o-2{ji 32)- We have due to the Clebsch-Gordan-Theorem |ji —32] < 02(^1 ji) < ji +J2- Now couple 02, js to a resulting 
angular momentum 03(02 J3) where again |a2 — jal < 03(02 J3) < 02 + js. This continues until we arrive at the last 

recoupling step where J = aN^aN^i jpf) and thus {apj^i — j^l < J < aN-i + jN- If we work with gauge invariant 

recoupling schemes we have J — and therefore must have jjv o^v-i. 

The challenge is now to find, for a given spin configuration ji, . . . , jjv, all possible gauge invariant states, that is 
all recoupling sequences as illustrated below for which 

4"™) lofc-i - Jfel < akiak-i Jk) < {ak~i + Jk) =: 0^^"^ (2 < fc < TV - 1) (4.2) 

(with oi — ji) and additionally o„_i = jj^. Let us moreover define the dimension for each intermediate recoupling 
step as 

/ N {max) {min) , r A o\ 

dim{ak) ■■= a}. -a}. ' + 1 (4.3) 

This expression has to be understood iteratively, since for the recoupling step the {k — 1)^^ step provides one 
initial spin as illustrated in figure [H 

4.1.1 Estimate for the Number of Eigenvalues 

As described in the previous section the construction of a recoupling state | 02(^1 ^2) • ■ • OAr_i(o7v-2 3n-i) J{(1n-i Jn) M ) 
is a successive process in which the possibilities at the k^^ step are given by the inequality (|4.2p . Moreover due to 
(|2.39p in section [2.51 we may assume ji < ^2 < • • ■ < Jn < jmax- We can iteratively estimate di'm{ak) as 

dim(o2) < 2jN + 1 dim(a3) < 3jN + 1 ^ ■■■ ^ dim(oAr_2) < (N - 2)]^ + 1 (4.4) 

Note that due to gauge invariance dim(oAr_i) — dim(ojv — J) — 1. So the dimension of the gauge invariant subspace 
contained in the [rifc=i(2j'c + 1)] -dimensional tensor space is bounded from above by 

7V-2 N-2 

dim (|o2(ji n) ■ ■ ■ J(ojv-i jn) M)) < n • ^'^ + 1] ^ n • + ^] = + 1]"^"' • - 2)' (4-5) 

1=2 1=2 

which is at the same time a bound for the maximum number of eigenvalues contributed by a fixed spin configuration. 
The number of spin configurations ji < J2 < • ■ • ^ jN can be calculated as follows: Assume an iV-dimensional 
orthonormal lattice with cubes of edge length 1 as elementary cells. Choose a lattice origin, and associate each jk 
{k — 1 . . . N) with a direction of the cube at the origin. If we exclude jk — then we can get the number of spin 
configurations A^j_configs ^"-"^ ^ particular jtv by counting the number of lattice points as: 

^j-con«. = E E • • • E - tn^^m^^ + 1) • . . . • (^. + ^ - 2) = (^^; ;^ - (4.6) 

6jV-l = l ftN-2 = l b2 = l ^ ^ 
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{max) 
'5 



a 



(max) 



a 



(max) 
3 



4 



(max) 



(min) 



(min) 



(max) 



a 



{min) 



{min) 



a, 



{min) 



02 



0-3 



04 



05 



0-6 



a 



(max) 
N-3 



(max) 



(max) 



(min) 



(min) 



(max) 



(min) 



aN~3 aw-i un — J — 



Figure 2: Schematic presentation of a gauge invariant recoupling scheme 



where we have defined bk '■— 2jk and used that X]fc=i ^ ~ + Thus we can estimate the number A^evaisj^- of 

eigenvalues contributed by spin configurations Ji < J2 < • ■ • < Jn with fixed Jn- 



cvals 



2)!- 



bN + N-2 
1 



N 



1 W-l 



TV- 1 



[Zjn]^ \i + ... + N-2) + 



1 N-3' 



O 



1 Af-5' 



JN 



1 2N-4 



2'^-\N ~3) + -{N ~1){N - 



2) ) [jn] 



2N-5 



o 



JN 



1 2N-6 



(4.7) 



Here we have computed the coefficients of the two highest orders in ■ We conchide that the number of eigenvalues 
resulting from spin configurations ji < j2 ^ • ■ ■ ^ jN, jN fixed will be given by a polynomial in jjsi of degree 27V — 4. 
This result will be used in order to give a polynomial fit to the number of eigenvalues obtained in the eigenvalue 
computations of section [S) 



4.2 Sign Factor Combinatorics 

Beside the recoupling part in the definition (|2.25p . the sign factors e{IJK) provide sensitivity of the volume operator 
to the diffeomorphism invariant structure of the edge embeddings at a particular vertex. The following discussion 
summarizes this up to now not fully understood problem and sketches the starting point for a detailed analysis. We 
find that there is a fascinating interplay between gauge invariance and the spatial embedding of a given graph, showing 
a possible link between edge spins and the (spatial) diffeomorphism invariant information encoded in the edge structure 
of a graph. 

4.2.1 Gauge Invariance at the Operator Level 

At the operator level gauge invariance 

3 N 2 

iJ(tot)?--=ll{llJK) -0 (4.8) 

k=l K=l 
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implies that 



N-l 



J% = -Y.Jl Vfc = 1,2,3 



(4.9) 



L=l 



This can be used in (|2.25p in order to ehminate all qijN containing N as an edge labels as follows: Because qijx 
const • ^ijkJjJjJK have 



(3 = const- ^ e{IJK) ei-jkJ}JjJK 

I<J<K<N 

= const •[ ^ e{IJK) e^jkJ}J'jJK+ J2 <IJN) €^,kJ}J'jJ\ 

I<J<K<N I<J<(K=N) 

Using (|4.9p we can rewrite the second sum term on the right hand side of (|4.10p as 

Af-l 

I<J<{K=N) 



(4.10) 



I<J<N L=l 



N-l 



I<J<N L=l 



L=I+1 



L=J+1 



Configurations with 2 identical edge labels L — I or L ^ J will not contribute 
to the sum due to the antisymmetry of eijk 



L<I<J<N 



I<L<J<N 



I<J<L<N 



J},Jl] = for I ^L. 



L<I<J<N I<L<J<N I<J<L<N 



L<I<J<N 



I<L<J<N 



I<.J<L<N 



(4.11) 



Upon introducing a suitable relabelling of the sum variables in (|4.1ip , taking each sum term into the form of the first 
line in (|4.10p . we conclude that the volume operator, if acting on a gauge invariant spin network function, can be 
rewritten as 



V := 



V =: 



Z ■ Y WJK) - <JKN) + t{IKN) - e(/JiV)] quK 

I<J<K<N 



Z ■ Y (^{IJK) qijK with a{IJK) := e{IJK) - e{IJN) + t{IKN) - e{JKN) (4.12) 

/<J<A'<Af 



which is the final expression for the volume operator acting on gauge invariant recoupling states. We would like to 
mention here that (|4.12p can be interpreted as a kind of 'self-regulating' property of the volume operator when applied 
to certain 'pathological' edge configurations, as illustrated in the following paragraphs. 

Consider the following explicit construction of a configuration of N edges outgoing from a vertex v, where each 
ordered edge triple ei,ej,eK, I < J < K contributes with a negative sign factor e(I,J,K). Since we are only 
interested in the sign factor we can make simplifying assumptions, in particular we may arbitrarily choose certain 
numerical values. Consider the vertex v as the origin of a 3 dimensional coordinate system with axes x,y,z. Now 
consider a circle with radius r — I centered at ?/ = 1, parallel to the x-z plane. Let every edge tangent vector ex 
end at a point on the circle with coordinates (cos</)if, IjSmcj)^) and (f>K — 27r-^. Now one may check that for each 
ordered triple e/, ej, eK with I < J < K < N we have: 

, /• • ■ N . • r K-Ii . r K-Ji . r J - /i 

det (e/, ej, e^j = -4 • sm ^tt ^ J sm [^tt ^ J sm ^ J (4-13) 

Since for all arguments x of the sine functions we have < x < tt, all of these functions are > and therefore we get 

e(/, J, K) = sgn ^ det (e/ej, 6;^)^ = —1 for all ordered edge triples e/, ej^ex, I < J < K at v. 

As one can see, this edge configuration is quite special — it is as if all edges lie in one octant only (if we rotate the 
coordinate system) or on one 2 dimensional surface as the hull of a cone. Such an edge configuration would appear 1 
or 2, rather than 3 dimensional. However, upon imposing gauge invariance, we conclude from (|4.12p that such an edge 
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Figure 3: 'Pathological' edge configuration 



configuration would effectively contribute zero volume, coming from only the orientation of the edges and independent 
of the spin confiauratioi^^. 

In |11|.|12| an upper bound for the coherent state expectation value of the operator corresponding to the inverse 
scale factor in a Friedmann-Robertson- Walker universe was derived. This upper bound provided boundedness of the 
expectation value, independently of whether the coherent states are peaked on classically singular data. However, 
the authors were unable to exclude the growing of the expectation value as the valence of the vertex on which it is 
evaluated is increased. They supposed this to be an effect of a too rough estimate, assuming all e{IJK) factors to 
be identical 1 or -1. However by the explicit construction of the edge configuration shown in figure [3 one can see 
that there indeed exist graphs leading to an unbounded coherent state expectation value of the inverse scale factor 
operator. This problem is now cured upon imposing gauge invariance via (|4.12p . which is why we refer to this issue 
as a 'self-regulating' property of the volume operator. 

A similar problem arises with respect to a state as presented in figure [31 in principle we can produce volume 
eigenvalues which increase with the number of edge triples at a vertex. See the discussion in section 6 of |25| . 

4.2.2 Numerical Computation of Sign-Factor Configurations e Using a Monte Carlo Process 

Consider an 7V-valent vertex v. For each of the (^) triples of edges l</<J<i4r<iVwe have a sign fac- 
tor e{IJK). Hereafter we refer to an assignment of the (^) such sign factors at a vertex as a sign configuration 



e := {e{IJK)}\i<j<K<N - (e(123), e(124), . . . , e(12iV), e(134), . . . , e(13iV), . . . , e(iV - 2, iV- 1, A^)) . The remainin] 



question in the definition (|2.25p is what sign configurations can be realized in the tangent space at the N-veiten^^ 
ti g E at all. This turns out to be a highly nontrivial problem. 

Our first attempt to solve this problem is via a Monte Carlo simulation. We note that the set of directions of tangent 
vectors to the embedded edges at a vertex is equivalent to the set of points of a unit 2-sphere. Thus the embedding 
of an 7V-valent vertex is equivalent to an TV-tuple of points on the sphere. We can explore the collection of such 
realizable embeddings by randomly selecting sets of N points on the sphere, via a Poisson process (by which we mean 
simply that the points are 'sprinkled' at unit density with respect to the area element on the sphere). For each such 
'sprinkling' of the N points, we can compute the sign factors e{IJK), for each of the (^) triples \ <I<J<K<N. 
The particular algorithm we employ is to sprinkle N points, compute the resulting sign configuration, and then store 
in a database (hash table) the number of times this sign configuration has arisen so far. We continue this process 
until the smallest such exceeds a fixed threshold. (In practice this threshold can be set at 1, because the number 
of sign configurations is enormous.) 

One should note that the numbers computed by the above algorithm have no intrinsic meaning in and of 
themselves, since we might do a coordinate transformation on the unit sphere, which changes the probability measure 
used during the sprinkling process. Moreover this method is unable to give sign configurations e containing coplanar 
triples {UK) with e{IJK) = 0, because they form a set of measure zero with respect to area on the 2-sphere. This 
is the technical reason why the present work will not consider these coplanar sign configurations, and is restricted to 



Nevertheless, using this method we can compute all the sign configurations e containing only linearly independent 
triples, as demonstrated in section [4.2.31 below. Using (j4.12p it is then straightforward to compute the effective signs 
a{IJK) — 0, ±2, ±4 resulting from gauge invariance and the restriction e{IJK) = ±1: 

The following table details our results: 

In the next to last column we have listed the number of realizable a -configurations, which is the number of 
assignments of a{IJK) ~ ±4, ±2, or to edge triples ei,ej,eK with edge labels \ <I<J<K<N according 

^^The spin configuration must allow for gauge invariance at v of course. 

'^^Recall that S denotes the 3-dimensional spatial foliation hypersurfaces used in the (3+l)-formulation of General Relativity; see the 
companion paper |25) for details. 




e{IJK) = ±1. 
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valence N 
of V 


number 
of triples: 

(?) 


maximal 
number of 
theoretically 
possible e-sign 
configs 
T„,„^ = 2(3') 

max ^ 


number of re- 
alized e*-sign 
configs 

Tr 


fraction 
Tr 
T 

J- max 


total number 
of realized 
(7-configs 


multiplicity 

(7 = 

-config 


4 


4 


16 


16 


1 


5 


6 


5 


10 


1024 


384 


0.375 


171 


24 


6 


20 




23,808 


0.023 


16,413 


120 


7 


35 


235 


3,486,720 


1.015 • 10"^ 


3,079,875 


720 



Table 1: Possible sign factor combinations a vertices of valence 4-7. 



to gUl). They are defined as a := {a{IJK)}\i<j<K<N = ((t(123), cr(124), . . . , cr(12iV - 1), cr(134), . . . , cr(13iV - 
l),...,a{N-3,N-2,N~l)). 

Note that there is a unique map S : e a where a :~ {a{IJK)}\ ,ej,cj^eE(v) . However S is not invertible, in 

i<j<K 

particular it will happen that for different e*, e* '-configurations we have S{e) = S{e') ~ a. We denote the number of 
e-configurations giving the same (i-configuration by Xff and call Xff the redundancy of a. 

In the last column we have written down the number Xa=d of distinct sign configurations e which lead to a 
(T-configuration for which (t{IJK) = yi<J<K<N. Note that this number seems to be given by {N — 1)1. 

4.2.3 Analytical Computation of the Sign Configurations e 

The results from the previous section are somewhat unsatisfactory, in that one may be concerned that the Monte Carlo 
process does not reveal every realizable sign configuration. For example how do we know that we do not terminate the 
sprinklings too early? Additionally, the Monte Carlo technique will never find linearly dependent triples of tangent 
vectors. Therefore it is highly desirable to have an analytic proposal at hand in order to compute the number of sign 
configurations directly, without using Monte Carlo methods. Here we present such a proposal, which will be worked 
out in more detail in a forthcoming paper [26j . 

In what follows we are going to derive a set of inequalities that must be fulfilled by a sign configuration in order for it 
to correspond to a realizable configuration of tangent vectors, meaning it must be embeddable in 3-dimensional Rieman- 
nian space. Let us slightly change the notation in the remainder of this section for the sake of clarity: Denote the valence 
iV of a vertex w by n in order to not confuse it with an edge labelling. Moreover let us use s[IJK] := e{IJK) in order 
to clearly distinguish the original signs e{IJK) from new signs which are introduced below. By its definition we have 
s[IJK] = sgn( det (e/(w), ej(D), ej<-(t;))) and for the tangent vectors ei{v) we introduce the shorthands ej{v) =: vj = 
'v}\ 

Vj , which are vectors in the tangent space to v. Then the determinants det (e/(w), ej{v), e/f (u)) can be expressed 
as 

det {ei{v), ejiv), eK{v)) =: s[IJK] ■ Q[IJK] = ^ e-ijk i^/WjWif (4-14) 

i j k 

where Q[IJK] :— \det{ei{v),ej{v),eK{v))\ > by construction. For an n-valent vertex v we have (3) = 
^n{n — l){n — 2) ordered triples with according sign factors. 

Let us for simplicity assume linear independence of all tangent vectors at the vertex v, that is Q[IJK] 7^ 
V / < J < < 71 in what follows, which implies that s[IJK] ^ and s[IJK] = ±1 for all ordered triples UK. 
This ignores of course a certain class of configurations, however the calculation presented here can be extended to the 
linearly dependent case in an obvious way, c.f. section [4.2.51 

Without loss of generality we may then introduce a basis in the tangent space of v given by the first three tangent 
vectors vi,V2, V3 as 

vi=(^Q^ 1/2= 1^1 j 1/3= j (4.15) 

Note that (|4.15p implies that s[123] = 1. This is a restriction since it fixes a positive orientation of this triple. However 
we will express all the other vectors with respect to that choice. Configurations in which ■i/i,i72,'y3 have negative 
orientation s[123] = —1 will give a global minus sign to all s[IJK]. We will work with the choice s[123] = 1 in what 
follows but we keep in mind that for every sign configuration we are going to discuss with respect to s[123] = 1, 
there exists a corresponding second one with inverted signs and s[123] = — 1. It will be sufficient to discuss only the 
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s[123] = 1 case, apply the found selection rules for realizable sign configurations, and multiply the obtained number of 
realizable configurations by 2. The second restriction Q[123] = 1 implied by the choice (|4.15p is just the normalization 
of the coordinate volume to 1. This only gives a positive global numerical factor which does not affect the signs s[IJK]. 

Employing all these conventions we can cast (|4.14p in the following form: 

^ global factor 2 





s[f23] Q[123] = 


1 






s[12L] Q[12i] = 


vl 




(II) 


s[f3L] Q[f3i] = 






(III) 


s[23L] Q[23L] = 


vl 




(1) 


s[lLM] Q[1LM] = 


-vhli 


+ vlvlj 


(2) 


s[2LM] Q[2LM] = 




- vlvli 


(3) 


s[3LM] Q[3LM] = 


-vhli 


- vivli 


(4) 


s[LMN] Q[LMN] = 


-vhli 





3 < L <n 



3 < L <n 



3 < L <n 



3 < L < M <n 



3 < L < M <n 



3 < L < M <n 



(4.16) 



'^If - vlvl,v% + vlvfivlj 



vivlivjf + v]vl,vl 



3<L<M<N<n 



The relations (|4.16p may be rewritten by eliminating the vector components in (1) . . . (4) by plugging in (/), (//), {HI)- 
Then 



(1) 
(2) 
(3) 

(4) 



s[lLM] Q[1LM] = -vIvIj+vIvIj = s[12i] s[13A/] g[12L] Q[13A'f] - s[13L] s[12M] Q[13L] 0[12Af] 
s[2LM] Q[2LM] = vlv\^ - vlvlj = s[12i] s[23A/] g[12L] Q[23Af] - s[23L] s[12A/] Q[23L] Q[12Af] 
s[3LM] Q[3LM] = -vlvlj - vlvfj = s[13i] s[23Af] g[13L] Q[23Af] - s[23L] s[13A/] Q[23i] g[13Af] 



s[LMN] Q[LMN] = -vlvj^^ + vlvljV% - vlvl^vf^ + - vlvljV% + vlvj^v^ 

= s[12L] s[13M] s[23N] Q[12L] Q[13M] Q[23N] - s[12L] s[23M] s[13N] Q[12L] Q[23M] Q[13N] 
+ s[13L] s[23M] s[12N] Q[13L] Q[23Af] Q[12iV] - s[13L] s[12A-f] s[23N] Q[13L] Q[12M] Q[23N] 
+ s[23L] s[12M] s[13N] Q[23L] Q[12M] Q[13N] - s[23L] s[13M] s[12N] Q[23L] Q[13M] Q[12N] 



(4.17) 



A further simplification can be brought to (|4.17p if we realize that our assumption s[IJK] ^ \/I<J<K<n 
implies that ^ V K — 1 . . .n, k = 1,2,3 due to equations {I), {II), {III) in (|4.16p . We arc therefore free 

to set = (3[12i] — 1 which only scales the individual vector!^ v^ and thus the moduli Q[IJK] but docs not 
change the direction v^ points to. That is setting = 1 does not affect the signs s[IJK]. Additionally we define 



XL ■= = Q[23L] and tjl 



Q[13L], so we can finally rewrite (|4.17p as 



(4.18) 



(1) Q[1LA4"] = s[lLM] [s[12L] s[13M] yu - s[l3L\ s[l2M] vl] 

(2) 0[2LAf] = s[2iAf] [s[12i] s[23A/] xm - s[23L] s[12M] xl] 

(3) Q[3LM] = s[3iAf] [s[13i] s[23A./] yLXu - .s[23i] s[13A/] xlVm] 

(4) Q[LMN] = s[LMN] [ s[12L] s[13M] s[23N] yM xn - s[12L] s[23Af] s[13iV] xm Vn 

+ s[13i] s[23Af] s[12A^] yL xm ~ s[l3L\ s[l2M] s[23N] yL xn 
+ s[23i]s[12Af]s[13A^] a;L?/Ar - s[23L] s[l3M] s[l2N] XLyM ] 

where we have used the fact that ( s[IJK] )^ = 1 in order to bring all sign factors to the right hand side of (|4.18p . 



^^This can be equivalently thought of as scaling the vector vl by dividing all its components by giving a new vector v'^ pointing 
in the same direction as but with components = t— ^ ■ 
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Now since by construction Q[IJK] is a manifestly positive quantity we can express l|4.18p as a system of inequalities: 

(1) < s[lLM] [s[l2L\ s[liM] yu ~ s[13i] s[12M] t/l] 

(2) < s[2LM] [s[12i] s[23Af] xm - s[23L] s[12Af] xl] 

(3) < s[iLM] [s[13L] s[23M] yLXu - s[23L] s[13M] xlVm] 



(4) < s[LMN] [ s[l2L\ s[13M] s[237V] yM xn - s[12L] s[23M] s[13iV] xm Vn 
+ s[13i] s[23Af] s[l2N] yL xm - s[13i] s[12M] s[23N] yL xn 
+ s[23L]s[12M]s[13N] XLVN - s[23L] s[13M] s[12N] xl yM ] 
Now we divided (1) by yn > 0, (2) by xm > 0, (3) by {xm ■ 2/m) > and (4) by {xn ■ Vn) > 

(1) < s\lLM] \s\12L] s\lMI] - s\13L] s\12M] —1 

VM 

(2) < s\2LM] \s\12L] s\23M] - s\23L] s\12M] —1 

XM 

(3) < s[3LM] [s[13L] s[23M] — ~ s[23L] s[13M] —1 

Vm Xm 

(4) < sVLMN] \ s[l2L] s[l3M] s[22,N] — - s\12L] s\23M] s\13N] — 

Vn Xn 

+ s[13L] s[23M] s[12N] - s[l3L\ s[l2M] s[23N] — 

UN Xn yN 

+ s\23L]s\\2M]s\l3N]— - s[23L] s[13Af] s[12iV] — ^ 1 

XN XN VN 

As last step we introduce the quotients 



(3) < s[3LM] [s[l3L\ s[23M] Yii - s[23L] s[13M] X 



(4) < s[LMN] [ s[\2L\ s[\3M] s[23N] Y^' - s[12L] s[23M] s[13N] X^^ 



N 



(4.19) 



(4.20) 



^ = (Yn^" := ^ = [XfY^ (4.21) 

Vm Xm 



with the obvious property that for L < M < N we have transitivity, that is 

= Yti ■ and X^ = X^j ■ X%' (4.22) 

We can then bring (|4.20p into its final form 

(1) < s[liAf] [s[12L] s[13Af] - s[l3L] s[l2M] Y^j] 

(2) < s[2iAf] [s[12L] s[23Af] - s[23L] s[12Af] Xf^] 



(4.23) 



+ s[13i] s[23Af] s[l2N] YI^jY^' X%' - s[13i] s[12Af] s[23A^] YI^jY^ 
+ s[23L] s[12Af] s[l3N] X'j;,jX^' - s[23L\ s[l3M] s[l2N] X^,iX^'Y^' ] 

We will subsequently discuss the properties of first the set (1), (2), (3) and secondly (4) as contained in (|4.23p in 
two separate subsections. 

The system (1), (2), (3) 

If we introduce the short hand signs 



= s[lLA/] s[12L] s[13Af] 4^^'^' = s[2LM] s[12L] s[23M] pf^^^ = s[3LM] s[l3L] s[23M] 
= -s[lLA/] s[13L] s[12Af] pf"^^ = -s[2LM] s[23L] s[l2M] pf^^^ = -s[3iAf] s[23i] s[13Af] 



*Note again that due to our conventions yK > and > \/K = 1 . . .n. 
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we can rewrite the equations (1), (2), (3) of (|4.23p as the simple system 



(1) 


< 




r P2 '^M 


(2) 


< 




r Pi 


(3) 


< 


Pb ^ M 


, ^LM] 

^ Pe 



(4.25) 



Each of the 3 inequahties in (|4.25p has 4 different sign-combinations: 



Solutions to inequality (1) 



Solutions to inequality (2) 



Solutions to inequality (3) 



case 


(LM) 
Pi 




condition 




case 


ALM) 

Ps 


ALM) 
Pi 


condition 


1.1 


+ 


+ 






2.1 


+ 


+ 


0<^M 


1.2 


+ 




< < 1 




2.2 


+ 




< Xf^ < 1 


1.3 




+ 






2.3 




+ 




1.4 






no solution 




2.4 






no solution 



case 


(LM) 

pI 




condition 


3.1 


+ 


+ 




3.2 


+ 




vL ^ vL 
^ ^ M 


3.3 




+ 


^ti < ^ii 


3.4 






no solution 



(4.26) 



It remains to discuss the complete set of solution^ for the coupled inequality system (1), (2), (3) in (|4.25p . where 
we will exclude those sign configurations which contain the cases 1.4 or 2.4 or 3.4 of ()4.26p since they do not contribute 
a solution. One finds the following table of solutions. It has to be read as follows: Each row gives a combination of 
3 solutions contained in (|4.26|) (for example the first row indicates a combination of solution 1.1 for inequality (1), 
2.1 for inequality (2), 3.1 for inequality (3) as given in (j4.26|) ). Then the resulting sign combinations are written as 

= 1 and '— ' = —1. Finally the conditions implied by the given set of solutions are written down. 



^^That is, consistent combinations of the signs. 
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{'-) 






[LA/] 
Pi 


[LM] 


[LM] 

Pa 


JLM] 
Pi 


Pb 


JLM] 

Pe 


conditrions given by (^}!(^) 


finSyl condition 


1 


1 


1 


+ 


+ 


+ 


+ 


+ 


+ 


o<y^^|o<xf, 




1 


1 


2 


+ 


+ 


+ 


+ 


+ 




0<YL\0<X'^\Xj^<Y^ 




1 


1 


3 


+ 


+ 


+ 


+ 




+ 


Q<Yti\G<Xlj\Yl;i<Xj;j 


0<Y^<X'^j 


1 


2 


1 


+ 


+ 


+ 




+ 


+ 


o<yj,^|o< x^<i 


o<y^^^|o< xj^Ki 


1 


2 


2 


+ 


+ 


+ 




+ 




^<ytl\^<^M<A^M<ytl 


'^<yk\^<^M<A^M<yti 


1 


2 


3 


+ 


+ 


+ 






+ 


^<^tl\^<^M<A^tl<^M 


o<yM<^M<^ 


1 


3 


1 


+ 


+ 




+ 


+ 


+ 


^<yti\^<^M 


0<Yii\l<X^j 


1 


3 


2 


+ 


+ 




+ 


+ 




^<ytiV<^M\^M<yti 


1 < ^M < ^M 


1 


3 


3 


+ 


+ 




+ 




+ 


^ <^m\^<^m\^M <^M 


0<nf<^M|l<^M 


2 


1 


1 


+ 




+ 


+ 


+ 


+ 


o<y^^,<i| o<xf, 


o<yj^<i| o<x|-. 


2 


1 


2 


+ 




+ 


+ 


+ 




0<Yti<l\ 0<Xlj\X'^i<Ytj 


0<X^j<Yii<l 


2 


1 


3 


+ 




+ 


+ 




+ 


o<yij<i| o<xf,|y^^<xf, 




2 


2 


1 


+ 




+ 




+ 


+ 


< y^^ < 1 < x^j < 1 


o<yA^<i| o<x|:,<i 


2 


2 


2 


+ 




+ 




+ 




o<yij<i| o<x^<i\x^<Yli 


0<^M<^Af<l 


2 


2 


3 


+ 




+ 






+ 


o<yij<i| o<x^<i\Yii<x^ 


0<Ym<^m<^ 


2 


3 


1 


+ 






+ 


+ 


+ 


0<Y^i<l\ KX^f 


0<n/<l<^M 


2 


3 


2 


+ 






+ 


+ 






no solution 


2 


3 


3 


+ 






+ 




+ 


o<y^^,<i| i<x^\ Y^j<x'^j 


0<Y^^<l<X^j 


3 


1 


1 




+ 


+ 


+ 


+ 


+ 


1<Y^,\ 0<X^j 


l<Y^j\ 0<X^ 


3 


1 


2 




+ 


+ 


+ 


+ 




l<Yii\ 0<X^,\ X^jKY^j 


0<X^.j<l<Yii 


3 


1 


3 




+ 


+ 


+ 




+ 




^<yti<^M 


3 


2 


1 




+ 


+ 




+ 


+ 


KY^^l 0<Xf,<l 


0<^M<l<nf 


3 


2 


2 




+ 


+ 




+ 




KY^jI 0<X^j<1\ X^jKYii 


0<^M<l<^Af 


3 


2 


3 




+ 


+ 






+ 


KY^jI 0<X^,<1| Y^j<Xj;j 


no solution 


3 


3 


1 




+ 




+ 


+ 


+ 






3 


3 


2 




+ 




+ 


+ 




^<Ym\ K X^j X^j < Y^j 




3 


3 


3 




+ 




+ 




+ 


l<Yt.,\ l<Xl,\ Y^,<Xi, 





By inspection of the definition of the p-signs in (I4.24p due to the arguments of sign factors we notice that 

pf''^-pf''^-...-pf''^=-l (4.28) 

and therefore in our discussion we wiU only have 12 cases contained in (j4.27p with an odd number of minus signs 
among the p's. Moreover we ignore the cases 2 — 3 — 2 and 3 — 2 — 3 since they possess no solution. The remaining 
solutions are given in the following table: 
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case 


(1) 




(3j 


Pi 


JLM] 


JLM] 

Pa 


Pi 


JLM] 

Ps 


Pa 


final condition 


1 


1 


1 


2 


+ 


+ 


+ 


+ 


+ 






2 


1 


1 


3 


+ 


+ 


+ 


+ 




+ 


<>V < 


3 


1 


2 


1 


+ 


+ 


+ 




+ 


+ 


o<yj^^|o<xf,<i 


4 


1 


3 


1 


+ 


+ 




+ 


+ 


+ 


0<^Af|l<^M 


5 


2 


1 


1 


+ 




+ 


+ 


+ 


+ 




6 


2 


2 


2 


+ 




+ 




+ 




0<^M<^M<1 


7 


2 


2 


3 


+ 




+ 






+ 




8 


2 


3 


3 


+ 






+ 




+ 


0<Yii<l<X^j 


9 


3 


1 


1 




+ 


+ 


+ 


+ 


+ 


KYiil 0<Xj;,, 


10 


3 


2 


2 




+ 


+ 




+ 




0<Xl,<l<Yt, 


11 


3 


3 


2 




+ 




+ 


+ 






12 


3 


3 


3 




+ 




+ 




+ 


1 < ^if < 



Note that for a 5-valent vertex [n = 5) we only need to consider the solutions 1 — 12 as provided by ()4.29p for the pair 
(L, M) = (4, 5), however for n > 5 we have for every triple L < M < N < n three relevant copies of the system (|4.25p 
which read as 



(1) 


LM] 


< 


[LMl 

Pi 


(2) 


LM] 


< 


JLM] 

Ps 


(3) 


LM] 


< 


Pb ^ M 


(1) 


[LN] 


< 


p'r 


(2) 


[LN] 


< 


P'r 


(3) 


[LN] 


< 


Ps " M 



(4.30) 



+ P2 ^M^N 

+ Pi ^M^N 

I Ji^^] vL vM 
N + Pe ^M^N 

wherein for the pair {LN) we have again used that = Y^^Y^^ and X^ = X^jX^ by definition. Using definition 
()4.25p of the p-signs we can write: 

p[^^-^l = s[lLM]s[l2L\s[liM] pf^^^l = s[lMN]s[l2Ad]s[liN] pf^^ = s[lLN]s[l2L\s[liN] 
pif ^'-'l = -s[lLM] s[13i] s[12M] p\^^^^ = -s[lMN] s[13M] s[12iV] p^^^' = -s[lLN] s[13L] s[12A^] 

p^^^] = s[2LM] s[12L] s[23M] p!,^^^' = s[2MN] s[l2M] s[23A^] pf^^ = s[2LN] s[12L] s[2?,N] 
pli^^^l = -s[2LM] s[23i] s[12Af] pj/^^^' = -s[2MN] s[23M] s[12iV] pf = -s[2LN] s[23L] s[12A^] ^^-^^^ 

pf^^^ = s[iLM] s[liL\ s[2iM] p[f^^l = s[?,MN] s[13Af] s[23A^] pf ^' = s[iLN] s[liL\ s[2iN] 
p|f = -s[dLM] s[2iL] s[liM] pjf^^l = -s[iMN] s[23M] .s[13iV] pf^^ = -s[2,LN] s[22,L] s[l?,N] 

which can be seen as a coordinate transformation from the sign factors s[LMN] to the p's. However, the definitions 
in (j4.3ip are redundant as the following three identities are obvious: 

^[LM]^[LM]^[MN]^[MN] ^ ,[i2L] .[13L] .[13iV] .[127V] ^ -pl^^lp^^"^' 

pf'%[^''^p[''''^pf"'^ = s[m s[2m s[2m s[i2n] = -pi^^Vf ^' 

^[LM]^[LM]^[Mm^[MN] ^ ^[^3^] ^[23L] s[23iV] s[13iV] = -pl^^lpl^^' (4.32) 
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A short MATHEMATICA calculation reveals that the coupled inequality system (|4.30l) has the following 312 
solutions if we impose the restrictions (|4.32p : 



LM 


MN 


LN 


1 


1 


6 


1 


1 


10 


1 


1 


11 


1 


2 


6 


1 


2 


7 


1 


2 


8 


1 


2 


10 


1 


2 


11 


1 


2 


12 


1 


3 


5 


1 


3 


9 


1 


4 


5 


1 


4 


9 


1 


5 


3 


1 


5 


4 


1 


6 


1 


1 


7 


1 


1 


7 


2 


1 


8 


1 


1 


8 


2 


1 


9 


3 


1 


9 


4 


1 


10 


1 


1 


11 


1 


1 


12 


1 


1 


12 


2 


2 


1 


6 


2 


1 


7 


2 


1 


8 


2 


1 


10 


2 


1 


11 


2 


1 


12 


2 


2 


7 


2 


2 


8 


2 


2 


12 


2 


3 


5 


2 


3 


9 


2 


4 


5 


2 


4 


9 


2 


5 


3 


2 


5 


4 


2 


6 


1 


2 


6 


2 


2 


7 


2 


2 


8 


2 


2 


9 


3 


2 


9 


4 


2 


10 


1 


2 


10 


2 


2 


11 


1 


2 


11 


2 


2 


12 


2 



LM 


MN 


LN 


3 


1 


5 


3 


1 


9 


3 


2 


5 


3 


2 


9 


3 


3 


6 


3 


3 


7 


3 


3 


10 


3 


4 


6 


3 


4 


7 


3 


4 


8 


3 


4 


10 


3 


4 


11 


3 


4 


12 


3 


5 


1 


3 


5 


2 


3 


6 


3 


3 


7 


3 


3 


8 


3 


3 


8 


4 


3 


9 


1 


3 


9 


2 


3 


10 


3 


3 


11 


3 


3 


11 


4 


3 


12 


3 


3 


12 


4 


4 


1 


5 


4 


1 


9 


4 


2 


5 


4 


2 


9 


4 


3 


6 


4 


3 


7 


4 


3 


8 


4 


3 


10 


4 


3 


11 


4 


3 


12 


4 


4 


8 


4 


4 


11 


4 


4 


12 


4 


5 


1 


4 


5 


2 


4 


6 


3 


4 


6 


4 


4 


7 


3 


4 


7 


4 


4 


8 


4 


4 


9 


1 


4 


9 


2 


4 


10 


3 


4 


10 


4 


4 


11 


4 


4 


12 


4 



LM 


MN 


LN 


5 


1 


3 


5 


1 


4 


5 


2 


3 


5 


2 


4 


5 


3 


1 


5 


3 


2 


5 


4 


1 


5 


4 


2 


5 


5 


6 


5 


5 


7 


5 


5 


8 


5 


6 


5 


5 


7 


5 


5 


8 


5 


5 


9 


6 


5 


9 


7 


5 


9 


8 


5 


9 


10 


5 


9 


11 


5 


9 


12 


5 


10 


5 


5 


10 


9 


5 


11 


5 


5 


11 


9 


5 


12 


5 


5 


12 


9 


6 


1 


1 


6 


2 


1 


6 


2 


2 


6 


3 


3 


6 


4 


3 


6 


4 


4 


6 


5 


5 


6 


6 


6 


6 


7 


6 


6 


7 


7 


6 


8 


6 


6 


8 


7 


6 


8 


8 


6 


9 


5 


6 


9 


9 


6 


10 


6 


6 


10 


10 


6 


11 


6 


6 


11 


10 


6 


11 


11 


6 


12 


6 


6 


12 


7 


6 


12 


8 


6 


12 


10 


6 


12 


11 


6 


12 


12 



LM 


MN 


LN 


7 


1 


1 


7 


1 


2 


7 


2 


2 


7 


3 


3 


7 


4 


3 


7 


4 


4 


7 


5 


5 


7 


6 


6 


7 


6 


7 


7 


7 


7 


7 


8 


7 


7 


8 


8 


7 


9 


5 


7 


9 


9 


7 


10 


6 


7 


10 


7 


7 


10 


10 


7 


11 


6 


7 


11 


7 


7 


11 


8 


7 


11 


10 


7 


11 


11 


7 


11 


12 


7 


12 


7 


7 


12 


8 


7 


12 


12 


8 


1 


1 


8 


1 


2 


8 


2 


2 


8 


3 


3 


8 


3 


4 


8 


4 


4 


8 


5 


5 


8 


6 


6 


8 


6 


7 


8 


6 


8 


8 


7 


7 


8 


7 


8 


8 
8 


8 
9 


8 
5 


8 


9 


9 


8 


10 


6 


8 


10 


7 


8 


10 


8 


8 


10 


10 


8 


10 


11 


8 


10 


12 


8 


11 


8 


8 


11 


11 


8 


11 


12 


8 


12 


8 


8 


12 


12 



LM 


MN 


LN 


9 


1 


3 


9 


1 


4 


9 


2 


3 


9 


2 


4 


9 


3 


1 


9 


3 


2 


9 


4 


1 


9 


4 


2 


9 


5 


6 


9 


5 


7 


9 


5 


8 


9 


5 


10 


9 


5 


11 


9 


5 


12 


9 


6 


5 


9 


6 


9 


9 


7 


5 


9 


7 


9 


9 


8 


5 


9 


8 


9 


9 


9 


10 


9 


9 


11 


9 


9 


12 


9 


10 


9 


9 


11 


9 


9 


12 


9 


10 


1 


1 


10 


2 


1 


10 


2 


2 


10 


3 


3 


10 


4 


3 


10 


4 


4 


10 


5 


5 


10 


5 


9 


10 


6 


6 


10 


6 


10 


10 


7 


6 


10 


7 


7 


10 


7 


10 


10 


8 


6 


10 


8 


7 


10 


8 


8 


10 


8 


10 


10 


8 


11 


10 


8 


12 


10 


9 


9 


10 


10 


10 


10 


11 


10 


10 


11 


11 


10 


12 


10 


10 


12 


11 


10 


12 


12 



LM 


MN 


LN 


11 


1 


1 


11 


2 


1 


11 


2 


2 


11 


3 


3 


11 


3 


4 


11 


4 


4 


11 


5 


5 


11 


5 


9 


11 


6 


6 


11 


6 


10 


11 


6 


11 


11 


7 


6 


11 


7 


7 


11 


7 


8 


11 


7 


10 


11 


7 


11 


11 


7 


12 


11 


8 


8 


11 


8 


11 


11 


8 


12 


11 


9 


9 


11 


10 


10 


11 


10 


11 


11 


11 


11 


11 


12 


11 


11 


12 


12 


12 


1 


1 


12 


1 


2 


12 


2 


2 


12 


3 


3 


12 


3 


4 


12 


4 


4 


12 


5 


5 


12 


5 


9 


12 


6 


6 


12 


6 


7 


12 


6 


8 


12 


6 


10 


12 


6 


11 


12 


6 


12 


12 


7 


7 


12 


7 


8 


12 


7 


12 


12 


8 


8 


12 


8 


12 


12 


9 


9 


12 


10 


10 


12 


10 


11 


12 


10 


12 


12 


11 


11 


12 


11 


12 


12 


12 


12 



In this table each solution for the coupled inequality system ()4.30p is given by a consistent combination of solutions 
for each of the three subsystems [LAI], [MN], [LN]: Each such solution of the subsystems is labelled according to 
the cases of the first column in (|4.29p . For example, the first row in the left table states that one solution of (|4.30p 
is obtained if we combine the following cases given bv (|4.29p : case 1 for the [LM]-subsystem, case 1 for the [MA^]- 
subsystem and case 6 for the [LA^J-subsystem. 
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The inequality (4) 

Let us now additionally consider inequality number (4) contained in (|4.23p which is given as: 
< s[LMN] [ s[12L] s[13Af] s[23iV] Y^' - s[12L] s[23Af] s[13iV] X^' 

+ s[13L] s[23A/] s[l2N] Y^jY^^ X^^ - s[l3L] s[l2M] s[23N] Y^jY^^ (4.34) 
+ s[23i] s[12A/] s[13N] X'^jX^' - s[23i] s[13A/] .s[12A^] X'^jX^'Y^' ] 
This can be rewritten, upon using the identities (|4.3ip . as follows: 

[MN] , „nori .rQ)i,f^n ^[^-f^l yM 



< s[LMN] [ s[12i] s[3AfiV] p^'"' Y^^ + s[12L] s[UdN] p^^'""^ X 



N 



- s[lZL] s[2MN] ^1 YkY^'X^' - s[13L] s[2MN] ^1 Y^Y^' (4.35) 

+ s[23L] s[lMN] pf'''^ XijX^' + s[23L] s[lMN] p[''^^ XijXffY^" ] 
Collecting terms with equal prefactors we arrive at: 

< s[LMN] [ s[l2L]s[iMN]{pf"'^ Y^' + pf""^ X^) 

- ,s[13L].[2Af7V]yM'(pf ^1 X^^ +pf ^1 ) (4.36) 
+ .[23L] s[lMN] Xi,X{f (pf ^V^^ + pi''"^ ) 

Here we notice that the terms in the round brackets precisely resemble the terms on the right hand side of the 
inequalities (1)[*^^1, (2)1^^^!, (3)[^^^1 contained in the inequality system (|4.30p . and are thus positive by construction 
if we want to obtain a solution. Let us abbreviate the sign-prefactors in (|4.36p by 

^[LMN] _ ^j^2i] s[3MN] = s[12L] s[l3M] s[23N] p'^^^^l 

^[LMN] _^[;l3^j s[2MN] = -s[13i] s[12M] s[23A^] 4*^^' (4-37) 
^[LMN] ^j23i] s[\MN] = s[23L] s[l2M] s[13N] p'*^^' 

where we have again used ()4.3ip on the right hand side. Let us take a closer look at products of sign factors on the 
right hand sides of (|4.37p . According to (|4.3ip the following identities hold: 

^[LN] ^[LN] ^ _5[i32,] s[23N] s[23L] s[13N] ^ s[23L] s[l3N] = -s[13L] s[23N] pf^^ pf^^ 

(4.38) 

^[LM] ^[LM] ^ _^[^2i] s[13Af] s[13L] s[12A/] ^ s[13i] s[12Af] = -s[12i] s[13Af] p'^^'^' p'^*^' 
Let us define the new sign-factor 

^[LMN\ ._ gj23i] s[l2M] s[liN] = s[12A/] s[23L] s[137V] 

= -s[13L] s[12Af] s[237V] pf^^ pf^^ 
= -s[2-iN] pf^^ pf^^ s[13L] s[l2M] 

= s[12i] s[13A/] s[2iN] pl^^l plf^l (4.39) 

Using (|4.39p we can now express the k's of (|4.37p completely in terms of p-signfactors and the overall sign A^^^^^': 

[LMN] _ .[LMN] MIN] 

— p^ 

^[LMN] ^ ^^LMN] p[LN] ^[LN] ^[MN] (4 

[LMN] _ .[LMN] JLM] [LM] [LN] [LN] [MN] 
"■3 ^ ^ Pi P2 P5 Pe Pb 

In the last step we use ()4.37p . ()4.40p to finally express ()4.36p in terms of p-signfactors and the overall sign A^^^^^l: 
< s[LMN] A[^*^^l [ pf''^ pf''^ ^1 (pf ^1 + pf ^1 

(4.41) 

[LN] [LN] [MN] I [MN] , WN] n [MN] ^L yM ( [MN]^M , JMN] 

+ Ps Pe Ps [Pi ^N + Pa ) + Pi ^M^N [P2 -'at + Pi 
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The advantage of reformulating (|4.36p as (|4.4ip is that we can give now a general tree of solutions to the total 
system (|4.23p for an arbitrary edge triple L < M < N < n: Using the p-signs we can label each such solution again 
by a solution to the coupled triple inequality system (|4.30p . contained in (|4.33p and additionally by an overall sign 
S := s[LMN]X^^'^^y For every solution of the triple system (IQU]) we seek compatible signs S := s[LMN]\^^^''^^ = 
±1. If a solution exists the sign 5* := s[LMN]\^'"^'^'^^ is written down in the subsequent table in the fourth column. 
We find that the total inequality system (|4.23p obeys 372 solutions given in the table below. 
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4.2.4 Results 



Tabic (|4.42p can be taken as a starting point to compute all possible sign configurations e for a given n-valent vertex 
V. Here it will be of particular interest whether one can get conclusions for high valent vertices from the consistency 
of all the triple solutions, such that as more and more e signs must be equal most of the a[IJK\ prefactors must 
vanish. We expect that it will be possible to develop an effective algorithm starting from (|4.42p A preliminary 
MATHEMATICA calculation reveals, for the numbers of linearly independent triples at the valence n = 4, 5, 6: 



n 


Number of triples 
(3) 


Number of potential 
sign configurations 

T — 2(3) 

J- max — ^ 


Number of 
realizable sign 
configurations 


Tr 
T 

J- max 


4 


4 


16 


16 


1 


5 


10 


1024 


384 


0.375 


6 


20 


220 


23,808 


0.023 



Table 2: Analytic check of the sprinkling results of table [T] using the inequality system method of this section. 

Comparing this table with the results of the Monte Carlo sprinkling process reveals that, in the case of 4, 5, and 
6- valent vertices, we have found all possible (non-coplanar) sign configurations from the sprinkling. 

4.2.5 Remark on Linearly Dependent Tangents 

In this section we sketch how edges having tangents which are linearly dependent can be treated. Given the set of 
tangent vectors of the edges at a vertex w, consider first only the maximal set of edges having linearly independent 
tangents using the procedure described above0- Calculate the possible sign configurations for this set of edges. Then 
consider the remaining edges having linearly dependent tangents for each of the obtained sign configurations as follows: 
For three linearly dependent tangent vectors vk,vl,vm we find constants aMK, olml S K, such that vm — olmk • 
vk + oiML ■ vl- In particular we have for the components v\j, i = 1,2,3 

vIj = aMK ■ Vr + o^ML ■ vl 

I 

sgn[vlf] ■ \vIj\ = sgn[aMK]sgn[v]i] ■ \aMKVK \ + sgn[aML] sgn[i;^] • |aMLwi| (4.43) 

These conditions then have to be imposed in the procedure above, giving further restrictions to possible sign configurations 
in (|4.42p on the complete set of triples. A detailed numerical computation based on the analysis presented here will 
be given in [55]. 



Such a set will always exist, since in order to have non- vanishing volume we must always have at least one linearly independent triple. 
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5 Computational Details 



In this section we discuss the details of our implementation of the numerical computation of the volume spectrum. 

Our code takes the valence iV of a vertex and the maximum spin jmax as input parameters. It computes eigenvalues 
for the A'^-vertex for all spin configurations ji < j2 < ■ • ■ < Jn < Jmax, subject to the constraint that J* 
integer, since otherwise there is no gauge invariant recoupling scheme. Since the edge spins take on a countable set of 
values, we store doubled spins as (1 byte) integers on the computer. 

5.1 Recoupling Scheme Implementation 

As discussed in section 14. 1[ we must compute the set of all possible recouplings of the edge spins to total angular 
momentum J — 0. We do this with a recursive procedure, as follows. Each basis vector a begins with ji as its 
first component. The second component 02 can potentially take any value |ji — J2I < 0-2 < ji + j2 which differs 
from those bounds by an integer. For each of these potential values, the third component 03 can take any value 
1^2 — is I < 02 < 02 + js (which also differs by an integer from the bounds), and so on. Thus the recursive procedure 
adds a value for the next component of a potential basis vector a. When it reaches aN-2 it considers only values for 
which it is possible to arrive at aN-i ~ Jn, which is required to get aN — J — 0. Any instances of the recursion for 
which this is not possible are abandoned at this stage. 

The particular algorithm employed stores each of the components of a obtained so far for each instance of the 
recursion. Since the majority of these instances will be abandoned at 'stage aM-2\ this is somewhat inefficient, and 
can cause problems when dealing with high valences (or particularly large spins). A more efficient use of memory 
would be to employ some tree-like algorithm which only stores, for each instance of the recursion, the potential values 
of the 'newest' component of a. 

Another issue with this algorithm is that it traverses every potential sequence of values for the spin recouplings, 
and only very close to the end does it abandon sequences which cannot possibly recouple to total angular momentum 
zero. It may be more efficient to make this determination as soon as possible, so the computer does not waste cycles 
traversing 'not possibly gauge invariant' recoupling schemes. 

It turns out that in practice neither of these issues make any difference, because the computable values of valence 
and spin are much more tightly bounded by consideration of the vast number of sign configurations e. The CPU 
time required to compute eigenvalues, and disk space required to store them, vastly outweigh any concerns about the 
efficiency of the recoupling algorithm. 

5.2 Realizable Sign Factors via Poisson Sprinkling Process 

The Monte Carlo computation of the set of realizable sign configurations is described in section l4.2.2l Here we comment 
on the Poisson sprinkling into S^. 

We describe points on S'^ via the usual {9, (/>) coordinates, and wish to select N points at random such that, for 
any region of the sphere of area A, the number of points sprinkled into that region is Poisson distributed with mean 
Given random variables x,y which are uniformly distributed in the unit interval, we select N points from 

9 = arccos(2a; — 1) , . 

This gives 9 G [0, n] and G [0, 2tt]. If we form an area element by taking the product of their differentials, we get 
{—) sm9d9d(j) = Andx dy, indicating that this gives the correct area element on the sphere. (The sign clearly has no 
significance) [lZ| 

5.3 Eigenvalue Computation and Numerical Errors 

As explained in section [2T3l we have a real antisymmetric matrix Q, or a real symmetric matrix Q^Q, whose eigenvalues 
we wish to compute. LAPACK ^35) is an extensive, efficient linear algebra library written in Fortran 77, available from 
Netlib [37] . One of the computations it performs is a singular value decomposition of a real matrix A into USV^, 
where U and V are orthogonal, and S is diagonal. The diagonal elements of S are the singular values. The singular 
decomposition of Q is 

Q = USV^ (5.3) 

To say this another way, we can regard 115.11 1 as a coordinate transformation on the 2-sphere. The area element in the {9, <f>) coordinates 
is related to that in the Cartesian coordinates by 

dx Ady = Jd9 A d4> (5.2) 
where J is the Jacobian of the coordinate transformation. Inverting 1)5. ip and expanding leads to the same equality as above. 
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Multiplying on the right with gives 

QQ'^ = USV^VS'^U'^ 
= USS'^U^ 

which is the diagonalization of QQ^ = QQ^ . Thus the eigenvalues of QQ^ are simply the square of the singular values 
of Q. As discussed in section [^751 these eigenvalues are simply the square of those of Q (which are pure imaginary and 
come in complex conjugate pairs). Thus the singular values of Q are non-positive real numbers and come in pairs, 
and their modulus is the square of that of the eigenvalues of V. (Note that due to the definition of V we are only 
concerned with the modulus of its eigenvalues.) 

LAPACK employs an efficient recursive algorithm to compute the singular values, which it returns with fixed 
relative accuracy, (e.g. as compared to the 'symmetric eigenvalue problem' (c.f. section TG.G.Sp . for which the errors 
scale with the 1-norm of the input matrix.) The absolute error of the singular values is bounded by 

(5.4) 

where e is machine epsilon (around 10"^^), s^ax is the largest singular value for that matrix, and c is a constant which 
governs the convergence criterion of the algorithm. For our calculations c ~ 90. 

The effect of these numerical errors on our dataset is extremely small, and is described in detail in section [6.61 

5.4 The Cactus Framework 

The Cactus framework |36| is a modular, open source, problem solving environment for high performance computing. 
It was originally developed to deal with the complex and computationally intensive problem of numerical relativity, 
where one seeks to extract the detailed form of gravitational waves which are emitted from astrophysical events such 
as binary star or black hole coalescence. The extensive needs of that computational task have led to what is now a very 
generic and extensible framework, which enables a wide class of applications to be easily placed on supercomputers. 
For this reason it seemed a natural framework to use for the development of code to explore Loop Quantum Gravity, 
and in particular to perform this volume spectrum computation. 

The name 'Cactus' is meant to describe the framework metaphorically as a collection of modules, called 'thorns', 
which are bound together by the 'flesh'. The thorns declare their characteristics, e.g. public variables called 'grid 
functions', to the flesh, which then generates code to allow the thorns to communicate and interoperate cleanly. To 
maintain portability Cactus defines its own datatypes, so it can guarantee their size and presence on all architectures 
it supports. 

Perhaps the greatest benefit which Cactus affords is automatic parallclization. Cactus contains a software layer 
called a driver, which hides the details of implementing parallel code from the application programmer. For our 
computation it allowed us to easily run on a distributed memory supercomputer, without much more difficulty than 
that required to develop the code to run on a single processor. 

The volume spectrum computation described above is composed of three stages. First, given a spin configuration, 
one must compute the recoupling basis, as described in sections |4 . II and [5.11 Next, with the recoupling basis in hand, 
one computes the {^^^) sub-matrices qijK of Q for each triple 1<I<J<K<N, shown in (|2.26p . These 
submatrices are shown in full detail in equations (|3.9p through (|3.22p . The final task is to combine the resulting 
submatrices, according to (|2.26p . for each (T-configuration (c.f. section l4.2.2p . This last task is easily parallelized, 
by splitting the set of a-configurations among the processors. This is implemented within Cactus by storing a a- 
configuration within the 64 bit integer type 'CCTK_INT8@ The thorn which implements the volume computation then 
declares a one dimensional CCTK_INT8 grid function, which is automatically distributed among the processors by the 
standard Cactus driver PUGH. 

A second benefit of using Cactus is that it provides a degree of language independence. The LAPACK library 
is written in Fortran 77, a language which is now quite old and inconvenient to work with. With Cactus one can 
write thorn code in C, C-l — h, or any variant of Fortran, and Cactus takes care of all issues regarding inter-language 
communication. All of the code developed for this Volume computation is written in C, save one routine in Fortran 
90 which performs all the calls to the LAPACK library. 

Additionally Cactus is 'familiar with' external packages such as LAPACK and BLAS (upon which LAPACK 
depends). It provides thorns which automatically detect and link to these libraries, which makes it extremely easy to 
port the code to a new machine. 

Last, but not least, is the modularity and extensibility of Cactus. The underlying philosophy behind its design 
is to facilitate code sharing and building a 'community code base' upon which all workers in a field can draw. A 
user is not required to release his or her thorns, of course, but the modularity of the thorn-ficsh design makes it 
easy to incorporate code developed by others in one's own code, so that one can build on the developments of others 

^*It turns out 64 bits is enough to store the tr-configuration for a 7- vertex, but not for an 8- vertex. The code could easily be generaUzed 
to handle higher valences. 
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rather than each user having to start from scratch. We intend to release the code which generated the results of this 
paper[^ It is hoped that this software will be a first step in constructing a community code base for performing high 
performance computations within Loop Quantum Gravity. 

6 Numerical Results on Volume Spectrum 

We will now show the results of our computations. It is in general remarkable that the (T-configuration will determine 
the shape of the spectrum as well as the behavior of the smallest non-zero eigenvalues (in particular the presence of a 
volume gap), and the largest eigenvalues as we increase the maximum spin jmax- 

6.1 General Remarks 

Let us begin by briefly reviewing the particular computation we have performed. Recall that given a graph 7 we have 
a set of edges E{-^) and a set of vertices 1^(7). Each edge e € E{-^) carries a representation of SU{2) with weight 
je, that is we have a set j := {je}e£E(j) of spins. In the gauge invariant setup we furthermore have at every vertex 
V G V{'^) an intertwiner ly projecting on to one trivial representation contained in the reduction of the tensor product 
of the 5[/(2)-representations carried by the edges outgoinj^ from v. Thus we have a set I :— We denote 

the set of N edges outgoing from v by E(v), and the set of spins incident at the vertex v by jy := {ji}\ei(^E{v)- By 
definition (|2.25p . (|2.26p . the action of the volume operator V can be separated to the single vertices v e V{'j). We 
can thus restrict our analysis to V acting on a single but arbitrary vertex v, since we can determine the action of V 
on arbitrary graphs from that. 

V is sensitive to the sign factors e{IJK) := sgn[det(e/, ej, ej<-)] which are the sign of the determinant of the 
tangents ei,e,j,eK of each ordered triple {ei,ej,eK), [I < J < K) oi edges outgoing from v. A particular sign 

configuration e := {e{IJKy}\ej,ej,e,^eE{^) is fixed by the chosen graph 7 (as a representative of its orbit under the 

i<j<K — 
diffeomorphism group). Note again that, upon imposing gauge invariance due to (|4.12p . there is a unique map 

S : e a where a := {a{IJK)}\cj,cj.ciieE(v) . Note that S is not invertible, in particular it will happen that for 

i<j<K 

different e, e '-configurations we have 5(e) = S{e ') = a. Let us denote the number of e-configurations giving the 
same (T-configuration by Xff- We call xs the redundancy of a. In the end we will compute the volume spectrum 
for diffeomorphism equivalence classes of vertex embeddings, each of which has a corresponding sign configuration 
e. Since only the (i-configurations enter into the eigenvalue computation, we must simply keep in mind that each 
(T-configuration comes with degeneracy xs- 

Given a fixed A'^-valent vertex v with fixed spin configuration jy — {ji < < • • ■ < jw} and fixed sigma 
configuration a, we compute the resulting matrix Q given in (|2.26p . As mentioned earlier we will employ Planck 
units, with £p, and the Z appearing in (|2.25p (which is a function of the Immirzi parameter /3 and the regularization 
procedure), set to 1. We then compute the eigenvalues Ag. The square root of the modulus of the Xq gives the 
eigenvalues of Xy of V. Note that the eigenvalues Ay = •\/|±Aq| come in pairs or are identical to 0. We store one 
number from each of the eigenvalues pairs, along with the corresponding spin and sigma configurations, in binary data 
filesS 

Ideally we would like to compute the full spectrum of the volume operator. Since its action on spin networks can be 
reduced to that on a single vertex, we consider only the latter. While in principle one must compute the eigenvalues 
above for all valences N and spin configurations j, it will be necessary, due to finite computational resources, to 
introduce a cutoff jmax on the maximum allowed spin in the spin configurations, and restrict ourselves to valences 
N < 7. We thus compute, for each valence N < 7, for every spin configuration jy = {ji < j 2 < • ■ • < Jn}, for every 
sigma configuration a described in section the set of eigenvalues Ag of the resulting matrix Q as given in (|2.26p . 
From that one arrives at the set of individual eigenvalues Xy for V. Each eigenvalue thus obtained is regarded as 
coming with a redundancy of D{j) ■ xs- The first factor arises from the degeneracy of the spin configuration (j2.40p as 
described in section [2.51 and the second from the degeneracy of the sigma configuration mentioned above. 

In the following we display the collection of eigenvalues we have computed in a number of graphs. Most of these 
take the form of histograms: the real line between and the maximum occurring eigenvalue A^'^'^^ is divided into A'bins 
equal width bins. The number of eigenvalues falling into each bin, taking into account the degeneracies mentioned 
above, is plotted. We ignore zero eigenvalues in all histograms. (In general there are many more zero eigenvalues than 
the number of eigenvalues in any histogram bin that we display.) We have attempted to choose the number of bins 
for the histograms such that we capture as much fine structure of this 'spectral density function' as possible, without 
using a too fine bin width such that the histogram appears 'noisy'. Figure d] illustrates the effect of bin width on the 
histogram. The histogram is like those described in section 16.3.21 showing the density of eigenvalues for a 5-vertex 
with ^-configuration^ (0 -2 -4 -4) and all spins up to jma.x — We see that decreasing the number of bins reduces 

^^Instructions for downloading the code will appear shortly. 

^"Recall that we can always redirect edges such that at every v £ V{-f) with valence greater or equal 3 we only have outgoing edges. 
These are stored in a native binary format. The results described in this paper come from 188GB of data. 
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Figure 4: Illustration of effect of bin width on histograms. 



the apparent noise in the curve. For a very small number of bins some fine detail of the density curve is washed out. 
Here, when increasing the bin width by a factor of 2, we also multiply the number of eigenvalues in that bin by ^, so 
that the curves line up with the same vertical scale. Throughout the rest of the paper, whenever two histograms of 
varying bin width appear in the same plot, we employ this same normalization technique, so that the histograms can 
be displayed with the same vertical scale. 

Additionally we will present various fits, which measure characteristic parameters of the eigenvalue set we have 
computed. For the number of eigenvalues as a function of j,„ax we employ a polynomial fit according to the estimate 
in section 14.1.11 Unless otherwise stated we have used for all other fits linear regression for logarithmic (in case of 
exponential dependencies) or double logarithmic (in case of possible polynomial dependencies) plots. This results from 
the fact that there is currently no applicable analytic model for the spectrum at hand, which could give us at least the 
functional dependencies of different characteristic spectral parameters. Therefore the fits we present here have to be 
seen rather as a suggestion to a possible behavior than as a proof. Also the error ranges we give should be considered 
within this context. 



6.2 Gauge Invariant 4- Vertex 

The 4- vertex has been studied extensively in [53] . New analytic results concerning the degeneracy of the spectrum and 
on the smallest non-zero eigenvalue are presented in section [T] Here we present only a histogram of the eigenvalues in 
figure m for completeness, as well as a fit of the total number of eigenvalues. 

Note that the 4- vertex has five realizable (T-configurations, with (t(123) able to take any of its possible values -4, -2, 
0, 2, or 4. The sign has no effect, so the only role of the a-configuration, in this case, is to scale the overall spectrum by 0, 
2, or 4. For consistency with the higher valences we include all eigenvalues contributed from these five (T-configurations 
in our analysis, whereas in [24] only the bare matrix 5123 was analyzed. Our numerical analysis extends the criterion 
for the smallest non-zero eigenvalue A^'"'' given in (24' to a higher maximum spin of jmax = = 63. This sequence 
is contributed by spin configurations ji — j2 — \ and = = jmax ■ As shown in [53] these configurations give rise to 

an eigenvalue sequence \^™^'^\imax) oc \J + !)■ This eigenvalue sequence as well as its spin configuration 
are in good agreement to the analytic lower bound (|7.55p for the eigenvalue^ Ag , as Ag > jmax ■ 

6.2.1 Histograms 

Our results for the 4-vertex are depicted in figure [3 which was generated with A^bins — 2048. Each line represents 
the histogram of the eigenvalues for all (T-configurations and all values of edge spins up to a maximum spin jmax (and 
all five (T-configurations). Thus the histogram for a given jmax includes all eigenvalues which were used to generate 
the histogram for jmax — 5- We will refer to this diagram as an 'overall' histogram, for the given valence. Here the 

22 Recall that = yj\ ± Aq|. 
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Figure 5: Overall 2048 bin histograms at the gauge invariant 4- vertex up to j„ 
eigenvalues in all, of which 17,364,833,136 are zero. 
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topmost line (gray) is the histogram for jmax = -f^, the red line beneath it is the histogram for j„ 
forth. There are 17,364,833,136 zero eigenvalues not included in the histogram. 



and so 



6.2.2 Number of Eigenvalues 



As a consistency check with the results in [21], we present, according to section 14.1.11 the total number number of 
eigenvalues (including zero eigenvalues) for all vertex configurations ji < j2 < js < ji — jmax- Fitting to a polynomial 
of degree four in jmax we obtain 



-'^evals Umax; 

with coefficients (and their 95% confidence interval) 



4 

E 

k=Q 



rk ■ (jmax)* 



(6.1) 



/ro = 54.4937 ±5.66522 \ 
ri = -117.403 ±1.22582 
ra = -10.6641 ±0.07806 
ra = 213.333 ±0.00184 

\r4 = 106.667 ±0.000014 / 



^ovais(jmax) " ^ovais (jmax)) takcs thc numerical value = 4533.14 for the 



For this fit thc quantity v^'ovaisumax; - ^Vvais 

set of 126 data points. 

As the number of eigenvalues is governed by the construction rules for the recoupling scheme basis, one can also give 
a model having a minimal number of parameters in order to get the behavior of iVovais in the leading order in jmax as 
^cvais (imax) = ^ ■ (jmax)^- In this casc wc obtain (using the non linear fit routine of Mathematica) r = 128.981 ± 0.428 
and s = 3.96158 ± 0.00082, = 1-18 • 10^^ which approximates the data quite weU. 

Note that in [23 the cumulative number of eigenvalues was analyzed, that is all eigenvalues, including zero eigen- 
values, contributed from spin configurations < ji < j2 < ja < < Jmax- In the present paper we have excluded 
all spin configurations with — ji, because setting the spin of an edge to corresponds to deleting that edge from 
the graph it belongs to. Moreover the present analysis considers the redundancy of the cri23 = 0(6), ±2(4), ±4(1) 
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Figure 6: Number of eigenvalues at the 4- vertex. The 
sohd hue is a fit to a fourth order polynomial as given in 
the text. 



Figure 7: Number of eigenvalues at the 4-vertex in a 
double logarithmic plot using a two-parameter model. 
Again the computed data is represented by dots, the fit 
is a solid line. 



prefactor: In brackets we have given the number Xff of e* = (e(123), e(124), e(134), e(234)) configurations of linearly 
independent edge tangents which give rise to the according 17123 prefactor. As the sum over all X3 is 16, we include 
this additional factor in our number of eigenvalue counting, as opposed to |24j which does not consider this factor. By 
taking into account this difference we exactly reproduce equation (95) in |24j . 

6.3 Gauge Invariant 5- Vertex 

6.3.1 Histograms 

According to section B:.2.2l we have 171 tr-configurations. Due to the modulus taken in the definition (|4.12p . configura- 
tions which differ by an overall sign can be identified. Therefore (excluding the a{IJK) = V I < J < K configura- 
tion) we are effectively left with 85 nontrivial (T-configurations. (There is clearly a symmetry in the a-configurations 
whereby if one cj-configuration is realizable then its pair with all opposite signs is also realizable.) 

We present here the histograms for the 5-vertex, for each maximum spin jmax ^ These include all 85 a- 
configurations. It turns out that, as opposed to the gauge invariant 4-vertex, the smallest eigenvalue does not always 
increase with jmax ■ Rather this property is dictated by the particular a-configuration, as can be seen in the discussion 
of the smallest eigenvalues below. Here this manifests itself in the growth of spectral density with jmax at zero, as can 
been seen in figure [9l Figure [TOl demonstrates that, upon dividing the number of eigenvalues A^evais in each histogram 
bin by the total number of eigenvalues A'cvais(jmax) (as given in section [B.3.5p . one finds that this growth of the spectral 
density becomes more distinct at larger jmax ■ 

Note that this cannot be traced back to numerical errors (see section 16.61 for a detailed discussion), but is a 
property of the volume operator itself: Eigenvalues are accumulated close to zero, resulting in a non-vanishing and 
in the limit jmax 00 possibly infinite spectral density there. As we go from Ay = to larger Ay the eigenvalue 
density approaches a minimum and then goes over into a rising edge. This edge then finally decays into a long tail of 
the eigenvalue distribution, which can be traced back to the finite cutoff in jmax- 

Figure [TT] portrays 16 histograms, one for each value of jmax = • ■ • ; Here each bin occupation number iVbm 
is divided by the maximum bin occupation number iVj^^'^(jmax) of eigenvalues for that jmax, and in addition each 
eigenvalue is divided by the maximum eigenvalue at that jmax, A™"^(jmax)- Thus each of the histograms is normalized, 
such that we may see if they possess a universal shape which is independent of jmax- Here one can see that the 16 
normalized histograms each have a quite similar shape. 

6.3.2 Effect of Different (T-configurations 

To give a feel for the effect of the tr-configuration, we show in figure a sampling of histograms as in figure [HI but for 
various fixed (T-configurations. The first three a-configurations shown (which have at least one a = ±4) give a smallest 
non-zero eigenvalue which decreases with jmax (cf. section I6.3.4|) . The next two a-configurations give a smallest 
non-zero eigenvalue which increases with jmax- The last histogram is for the a-configuration which gives a constant 
smallest non-zero eigenvalue. Some details regarding these histograms are given in table[3] The Xff is defined in section 
16.11 Nzems is the total number of zero eigenvalues which occur for each <T-configuration. The last three histograms 
look sharper with fewer bins because there are more e-configurations corresponding to these (T-configurations, which 
causes more jitter in the spectrum. 

Note how the a-configuration affects many aspects of the spectrum, including the accumulation at zero, a dip in 
the spectral density at small eigenvalues around 10, the presence of a secondary peak at 30 or so, the presence of 
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Figure 8: Overall 2048 bin histograms at the gauge invariant 5-vertex up to jmax = There are 4,680,292,624,128 
eigenvalues in all, of which 450,641,502,288 are regarded as zero. 
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Figure 9: Zoom of eigenvalues close to at the gauge invariant 5-vertex up to j,, 
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Figure 10: Zoom of log-plot of eigenvalues close to at the gauge invariant 5-vertex up to jmax = We display 
the normalized histogram here: The number of eigenvalues per bin iVbin is divided by the total number of eigenvalues 
-^Gvais(jmax) (as given in section fG.S.Sp . such that the sum of all histogram bin occupation numbers gives 1. For larger 
jmax we clearly see a drastically increasing eigenvalue density close to zero. 
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Figure 11: Normalized fixed- jmax histograms for the 5-vertex for j„ 
superposed are shown in the left figure. 
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(right figure). All histograms 
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Table 3: Details of some fixed a-configuration histograms for the 5-vertex. 
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jagged edges for larger eigenvalues, and the magnitude of the largest eigenvalues which arise. 
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Figure 12: Histograms for individual (T-configurations of the 5-vertex. 

In figure [l2] we show the histograms for each (j-configuration, without imposing any constraints on jmax (save 
that it is < ^ ) . Here we have also ignored the a-configuration redundancy {xa of section I6.ip , and weighted each 
individual (T-configuration equally. Each histogram has 2048 bins, save those of the a-configurations corresponding to 
the (7-indices 30, 32-34, 36, and 40 from table |4] of section [6.3.4( which have only 512. These later ten histograms 
have their A'^cvais divided by 4, so they appear on the same scale as the other 75 histograms of figure [T31 This 'coarse 
graining' improves the clarity of the figure, which is otherwise obscured by the rapidly fluctuating A^cvais for these ten 
histograms. 



6.3.3 Fitting of Histograms 

The overall histogram for the 5-vertex, shown in figure [SI seems to possess an exponentially rising edge for eigenvalues 
around Ay = 10 to 30 or so. Note that as we consider larger and larger jmax, this edge does not change much as 
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Figure 13: Histogram for each of the eighty-five (j-eonfigurations of the 5-vertex. 

compared to the increase in the remainder of the histogram, save being extended to larger eigenvalues. Thus this 
exponentially rising edge may be a feature of the spectrum in the limit of jmax oo. Figure [M] shows an exponential 
fit to this rising edge. The range of data points to fit, as well as the form of the fitting function, are chosen arbitrarily, 
as there is no directly applicable analytic model for the eigenvalue distribution available at the moment (though in 
[T5] an expression for the 4- vertex was derived, see comments in Thus we would like to emphasize that the errors 

given for fits in the spectrum throughout this paper, although meaningful within the chosen set of data points, have 
to be taken in the context of the choice of fitting model and data point set we have made by hand. Our point here is 
to demonstrate that the spectrum in principle allows for an exponentially increasing number of eigenvalues. The fit 
gives the logarithm of the number of eigenvalues A^ovais(Ay) as a function of the eigenvalue by: 



ln[7Vevals(A^)] 

20 data points 



^f^t) _ 



(14.87 ± 0.20) + (0.43 ± 0.02)A^ 



2.87 • 10^ • e°-^^ 



(6.2) 



The quantity ■= EfcLi (^ivaii(-^i/^) " ^cviis(-^y))^ ^as the numerical value = 1-06 • lO^^ on the chosen set of 



6.3.4 Extremal Eigenvalues 

We have seen in section 16.3.21 that the a-configuration has a considerable effect on the shape of the distribution of 
eigenvalues. Here we examine the largest and smallest eigenvalues which arise in the spectrum, as a function of jmax 
and (7-configuration. Figures (TK] and show the minimum and maximum non-zero eigenvalues of the spectrum, for 
a given (T-configuration and jmax- 

As it turns out from the numerical analysis many of the a-configurations give identical minimum eigenvalues A^'"' 
for each value of jmax, and we have separated them into equivalence classes, such that any two (i-configurations which 
have the same A^™''(jniax) will lie in the same equivalence class. The decision whether two (j-configurations belong to 
the same equivalence class is made as follows: Along with the smallest eigenvalue data we have recorded the numerical 
error bound of each eigenvalue according to section 15.31 We take the maximum numerical error which occurs among 
all the computed smallest eigenvalues as a threshold and look at the absolute differences in an elementwise subtraction 
of two eigenvalue lists. If the maximum numerical value of all the differences does not exceed the threshold then 
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Figure 14: Logarithmic plot of the fit to the rising edge of the spectral density for the gauge invariant 5-vertex. The 
points shown in red were used in the fitting process, the blue solid line is the fitted function. 



we consider the two (j-configurations as elements of the same equivalence class. These equivalence classes, which are 
detailed in tables [4] and [H are used in the figures [TH and [TH 




Figure 15: Smallest non-zero eigenvalues at the Figure 16: Largest Eigenvalues at the gauge invari- 
gauge invariant 5-vertex up to jmax = ^ ■ The a- ant 5-vertex up to jmax = ^ • The (?- index is defined 
index is defined table HI in table El 

By inspection of figure [TSl one can clearly see that some tr-configurations give rise to a minimum eigenvalue A^'"^ 
which decreases with jmax, some give rise to a A^'"'' which increases with jmax, and there is one in the middle for which 
^(inin) ^Qgg vary with j'max- These three behaviors of the minimum eigenvalues serve to classify the cr-configurations 
into three groups. In the following table we count the number of (j-configurations in each of these groups, and also 
the number of inequivalent sequences A^™'' (jmax) which arise from the a-configurations in each group. (We will often 

refer to ordered sequences of eigenvalues, such as A^""''(jmax), A^^'^'^'' (jmax), as jmax ^ oo, as eigenvalue sequences.) 
We findEl 



increasing u-conf. 


equiv. classes 


10 




total 


24 


decreasing cr-conf. 


equiv. classes 


29 




total 


60 


constant a-conf. 


equiv. classes 


1 




total 


1 



Note that the Ust of equivalence classes given in tables ^ and [S] may be regarded as preliminary, because the eigenvalue sequences 
max) for two different (X-configurations may match only up to some finite jmax- As observed already in |24l . the structure of the 
volume spectrum becomes richer as we increase the cutoff jmax, in particular new eigenvalue sequences may show up. Thus the number of 
equivalence classes could in principle increase at yet larger jmax. 
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The analysis of the data for the smallest non-zero eigenvalues reveals that increasing smallest eigenvalues are con- 
tributed from (T-configurations containing only (0, ±2) but not ±4, whereas decreasing eigenvalues are contributed by 
(T-configurations containing (0,±2,±4). Table [H lists the equivalence classes of a-configurations showing a particular 
behavior for A^"'"^(jmax) in detail. Recall the definition a :— {a-(123), (t(124), (t(134), (j(234)} where the arguments of 
a{IJK) are the edge labels for the particular edge triple {ei,ej,eK)- We have assigned an index '(T-index' to every 
equivalence class of a-configurations which relates it to the according curve in figure [TH 
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Table 4: Equivalence classes of a-configurations with regard to A^, 
accordance with the behavior of the smallest eigenvalues. 



(mill) 



(jmax)- The (T-configurations are grouped in 



The fact that there is a a-configuration giving constant smallest eigenvalues, independent of the value of jmax, 
seems to be puzzling, however this can be easily understood when we realize that this is the configuration a := 
{(t(123), (t(124), (t(134), (7(234)} = {2, 0, 0, 0}, as follows. If we look at equation (|3.9p . which gives the matrix elements 
for 9123; then the corresponding matrix only depends on the spins Ji,j2,j3 and the intermediate recoupling 02 and 
a^. So it is effectively the matrix of the gauge invariant 4- vertex. Now even if becomes large, there still exist spin 
configurations with a large J4. Recall that the intermediate recouplings in a recoupling scheme at the gauge invariant 

5-vertex are defined as: 02(^1 ^2), ^3(02 js), 04(03 ji) = jb where one has to apply the rules for the recoupling according 
to section 14.11 Now the smallest absolute numerical value for eigenvalues of the gauge invariant 4- vertex was found in 
[21] to be given by a configuration with ji ^ ~ 1 and J3 ~ 03 = J4 ~ jmax 13 But now we have this case embedded 
in a 5-vertex. This means that if js is tiny then 03 can be tiny as well (since ji < j2 < is < ^4 < is) whereas it might 
still happen that -I- ja ^ js- This case can occur for every js = jmax, since we can still freely choose ^4 J5 (because 
we consider all ji, . . . , j4-combinations up to j^). So we always get the same value. In the case of the 5-vertex, since 
the number of spins is odd, there is a sort of 'struggling' going on in the spin configuration which gives rise to the 
smallest eigenvalues, due to the condition that the sum of the 5 spins must be an integer (because otherwise we cannot 
recouple them to resulting zero spin). Thus different spin combinations can contribute the smallest eigenvalue: At the 
gauge invariant 4-vertex we saw that the smallest eigenvalue is contributed by spin configurations ji = J2 = 5 and 
j3 = j'4 = jmax- Now assume a similar statement is true for the gauge invariant 5-vertex, for exampleo something 
like ji = ^2 = 5 and J3 = J4 = J5 = jmax ■ Then it is easy to see that while in case of the 4-vertex it does not matter 
if jynax is integer or half integer, this does matter at the 5-vertex. The combination here can only be realized if jmax 
is integer, otherwise there is no recoupling to resulting spin possible. Thus for jmax half integer, a different spin 
combination will contribute the smallest non-zero eigenvalues. This is nicely demonstrated by the oscillating curves 

^^This result is proved below in section [7.4.21 

^^This is just a fictive example to illustrate our point. 
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in figure [TOl 



Overall smallest eigenvalues As can be seen from tlie analysis of the data, the overall smallest non-zero eigenvalue, 
for each value of j^ax, is always contributed by a := {cr(123), cr(124), cr(134), cr(234)} = {4,2,2,0} = {4,-2,-2,0} 
=: (^min- It is contributed by the spin configurations listed in table O (here js = jmax)- 
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Table 5: The spin configurations contributing the smallest eigenvalues for <7min- 

These minimum eigenvalues are plotted in figure fT71 with the eigenvalue axis in log scale. Note that the eigenvalues 
become comparable to the numerical noise at around jmax — ^ ■ This can also be seen in the breaking of the spin- 
pattern in table [5] for ^5 = jmax > By inspection of figure [IT] one can extract information on the behavior of A^'""* 
depending on jmax- A linear fit for a chosen set of data (see caption of figure \mi reveals that aJ™"')^^-^**) (2 • j^ax) ~ 
^(3.395±o.073) ^-(o.394±o.oo3)2.j_ ^ f^^j that EHi (A^™"^'^^'*^ (jf 2x) " x'f"'\&Y = 7.7 • IQ-io for the 

chosen set of 11 data points. 

In order to show the generic behavior of the smallest eigenvalue sequences we have included figure[Tni Note that the 
eigenvalue sequences are not always clearly separated from each other numerically, even when the eigenvalues are large 
compared to the numerical errors. They are also not always monotonically decreasing. This can be traced back to the 
following two reasons: Firstly certain eigenvalue subsequences show up only above a certain value of jmax, therefore 
we can expect a 'change of the smallest eigenvalue sequence' as jmax increases. Secondly, gauge invariance demands 
the sum of the five spins be an integer. Hence the spin configuration which gives rise to the matrix contributing 
the smallest eigenvalue can be forced to change, depending on the parity of jmax^ we may 'jump' between different 
eigenvalue sequences which contribute the smallest eigenvalue. These effects give rise to a sort of 'oscillating' curve, 
whose enveloping curve nevertheless still decreases. As the numerical noise is approached at jmax > ^ ^ clean 
distinction between the smallest two eigenvalue sequences is lost. Comparison with the numerical errors shown in 
figure [T71 and figure HSl indicates that this is caused by numerical noise. 
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10 20. 30 40 



Figure 17: Smallest eigenvalue sequence at the 
gauge invariant 5-vertex, contributed by (Tmin ■— 
{cr(123), cr(124), cr(134), cr(234)} = {4, 2, 2, 0}. This is a- 
index number 1 in figure [TSl 
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Figure 19: The twelve smallest eigenvalue sequences at the gauge invariant 5-vertex. The red line indicates a, 



Figure 18: Fit of the logarithm ln[A^^'""'^] of the small- 
est eigenvalue sequence of figure [Tjl We have used 
A^^^Omax) for jinax = for the fitting (blue 

dots). The red line indicates the eigenvalue sequence of 
f?min, the black line the fit function. 
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Table 6: Equivalence classes of a-configurations with regard to A^^'''(j„ 



Overall largest eigenvalues In the same way as for the minimum eigenvalues, we place the (T-configurations into 



equivalence classes according to the maximum eigenvalue sequence 



(max) 



Omax)- In this case there are only seven 



equivalence classes, out of the 85 ir-configurations. They are listed in table [H The tr-indices listed there are used to 
plot the maximum eigenvalue sequences A^^'^^ (jmax) in figure [TBI 

The overall largest eigenvalue sequence is contributed by those a-configurations which arc contained in the equiva- 
lence class of (T-index 7. These maximum eigenvalues are plotted in figurejini with the eigenvalue axis in log scale. It is 
interesting to note that the tr-configuration which gives rise to the smallest eigenvalues, a — {4, 2, 2, 0} = {4, —2, —2, 0}, 
also lies in the equivalence class of (i-configurations which gives rise to the largest eigenvalues. Figure [21] shows a fit to 

\(max) , (max),(/it) /„ • n _ ^0.067±0.043 Icy ■ \1.484±0.013 ,2 Y^29 f .{max) ,{fit) , .(k) n >(max)/.(fe) s\2 _ 

as [^■Jmaxj—'S ■[^■Jma.x) , Wltn X 2^k=l [■^y {.]raa.x)~A^ Umaxjj — 

57.3 for the chosen set of 29 data points. 



2j„ 



Figure 20: Largest eigenvalue sequence at the 
gauge invariant 5 vertex contributed by 
{cr(123),cr(124), (7(134), cr(234)} contained in the 
equivalence class with index number 7 in figure 1161 
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Figure 21: Fit of the logarithm ln[A " ^] of 
the largest eigenvalue sequence of figure [zDj versus 
ln[2ji„ax]- We have used a|™'"'-' (jmax) for j,„ax = 



V 

. . . ^ for the fitting (blue dots). The red line in- 
dicates the eigenvalue sequence, the black line the fit 
function given in the text. 
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6.3.5 Number of Eigenvalues 

According to (|4.7p the number of eigenvalues A^evais(jmax) for configurations with ji < j2 < ja < j4 < js = jmax will 
be given by a polynomial of the order of [jmax]'^^~'^ = [jmaxf- 

Nii:ll{j^.^)=j2rk ■ (Jmax)'^ (6.3) 
fc=0 

with coefficients (and their 95% confidence interval) 





= 106.59 


±22515.70 \ 


ri 


= 1523.19 


±26095.85 


r2 


= -164.42 


±9685.98 


rs 


= -8467.44 ±1581.95 




= 9014.52 


±126.55 




= 15678.82 


±4.86 




= 5226.69 


±0.07 / 



For this fit the quantity ■= S -i (^ovais(imax) — -^o^^^^j^ (jmax)) takes the numerical value — 2.26 • 10^ for 
the set of 44 data points. 




5 10 15 20 1 1.5 2 3 5 7 10 15 20 

Jmax Jmax 



Figure 22: Number of eigenvalues (dots) at the 5- Figure 23: Number of eigenvalues at the 5-vertex 
vertex, fitted by a sixth order polynomial (solid in a double logarithmic plot for the two parameter 
curve). minimum model. 

As in the case for the 4- vertex, one can also try a non linear two parameter fit A'^^vai's (jmax) — r ■ (jmax)^, in 
order to measure the behavior of iVovais in the leading order in jmax. In this case we obtain r — 9296.6 ± 96.02 and 
s = 5.85582 ± 0.00341, := EJ^^^i (^cvak(jmax) - Afevafi(imax))' = 1-34 • lO^^ for the set of 44 data points (see 
figure [ 
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6.4 Gauge Invariant 6- Vertex 



There are 16413 tr-configurations for the 6- vertex, according to section [4.2.21 As mentioned in section [B. 3. 41 configu- 
rations which differ by an overaU sign can be identified. Therefore (exciuding a{IJK) — V J < J < K) we are ieft 
with 8206 non-trivial tj-configurations. 

6.4.1 Histograms 




2e+09 



Figure 24: Overall histograms at the gauge invariant 6- vertex up to jr, 
in aU, of which 212,946,944,688 are zero. Nuns = 2048. 



120 



= ^. There are 7,349,844,794,112 eigenvalues 



As one can see in figure !^ the overall histograms for the gauge invariant 6- vertex show a similar behavior as those 
of the 4- and 5-vertex. However, there is no 'lip' in the spectral density close to zero visible (and quadrupling the 
number of bins does not alter this). Nevertheless the number of eigenvalues near zero increases with jmax and, as can 
be seen from section [6.4.31 below, the behavior of the smallest eigenvalue of a particular vertex configuration depends 
on its (7-configuration. Thus it seems that upon taking all eigenvalues from all a-configurations the eigenvalues at the 
6-valent vertex are rather equally spaced from on, whereas the spectrum of the 5-vertex has an accumulation point 
at 0. Nevertheless there are also single a-configurations which show a lip in their spectral density, as demonstrated in 
figure [m 

In figure [25] we show the normalized overall histogram for the 6- vertex, analogous to figure dni of section 16.3.11 in 
which each bin occupation number A'bin has been divided by the total number of eigenvalues A"cvais(jmax) at that jmax- 
There we see that, as opposed to the situation with the 5-vertex, the normalized spectral density does decrease close 
to zero. 

Figurel^ shows two histograms for fixed a-configuration a = (ct(123), (t(124), (t(125), (7(134), (t(135), a(145), cr(234), 
(t(235), (t(245), (t(345)). The (i-configuration on the left yields a minimum eigenvalue which increases with jmax, while 
the fT-configuration on the right yields an decreasing minimum eigenvalue. Some details regarding these histograms 
are listed in table [71 Figure [27l shows the portion of the histograms on the right side of figure [26] for Ay < 10. There 
we see the lip quite clearly, which indicates that it is the 'averaging' over all (T-configurations that causing the lip to 
flatten out in figure [Ml 
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Figure 25: Normalized overall histogram for the 6-vertex. The vertical axis is in logscale. The histograms for the 
smallest four values of j,nax have been removed to give a clearer picture of the behavior at larger j,nax- -^bins = 2048. 





Figure 26: Two sample histograms for fixed (j-configuration. On the left <j = (0, 0, 0, 2, 2, 2, 2, 2, 2, 2) and on the 
right a = (2, -4, -4, 0, 0, 0, 4, 2, 2, 0). 



(7-configuration 


Xa ^zcros A^bins bchavior 


02222222 
2 -4 -4 0004220 


8 169,083,736 256 increasing Aj^""^ 
2 2,998,262 512 decreasing a|.™"^ 



Table 7: Details of two fixed tr-configuration histograms for the 6-vertex. 
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Figure 27: Zoom of the lip of the decreasing A^'" histogram of figure [211 Here we use iVbins ~ 1024 to show more 
detail. 

6.4.2 Fitting of Histograms 

As in section 16.3.31 figure seems to possess a spectral density which rises exponentially with jmax- A fit to this rising 
edge is shown in figure 1281 in which the logarithm of the number A'cvais of eigenvalues as a function of the eigenvalue 
Ay is given by: 

ln[iVevais(Ay)](^**) = (19.081 ±0.007) + (0.549 ±0.008)Ay — > Af,^^^\*](Ay) - 0.19 • 10^ • e^-^^ (6.4) 

The quantity ^ := ^^Li (^evaii(^t>^) - ^fl\jK^v)t has the numerical value ^ = 7.257 • lO^"* on the chosen set of 
71 data points. 

6.4.3 Extremal Eigenvalues 
Smallest Eigenvalue 

Figure [29] shows how the smallest eigenvalue varies with jmax for a sample of 21 a-configurations. As is the case 
for the 5-vertex, some a-configurations lead to rising minimum eigenvalues, and some to falling. Figure [30] shows the 
same for a collection of ci-configurations which only yield a decreasing minimum eigenvalues. 

As it turns out, it is not possible to identify a genuine smallest eigenvalue sequence for the 6- vertex, at least 
not in the computed spin parameter range. This is due to the fact that the different eigenvalue sequences intersect 
each other and can hardly be separated. Figure [BT] illustrates this property. Nevertheless we can define a smallest 
eigenvalue sequence by taking into account the ordering of the smallest eigenvalue sequences at a fixed jmax- The 
chosen sequence is displayed as red curve in figure [311 it is obtained by comparing all eigenvalue sequences at jmax = ^■ 
It is contributed by the spin configurations shown in table [SI (here je — jmax)- 

We can then (as in case if the 5-vertex) apply a fitting of ln[A^'"^] versus 2jniax- A linear fit for a chosen 
set of data (see caption of figure [321) reveals that A^""^(jmax) = (B(3-723±i.ii9)(g-(i.2i5±o.i30)2-j„,,^ Yi&xQ — 
X^fcLi ('^^"™''''^**''(imax) — A^'"''(jm2x))^ — 0.32 for the chosen set of 11 spin configurations. 

Largest Eigenvalue 

In figure[nSlwe show the maximum eigenvalues for the 30 equivalence classes (according to the maximum eigenvalues) 
of (7-configurations. 
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Figure 28: Logarithmic plot of the fit to the rising edge at the gauge invariant 6- vertex. The points that were used in 
the fitting process are indicated in red, the fitted function is the solid blue line. 
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Table 8: The spin configurations contributing the smallest eigenvalues for amin ~ (2, 0, 0, 0, 0, 0, 0, 0, -2, 0), (2, 0, 
0, 0, 0, 0, 0, 0, 2, 0). 




Figure 29: Minimum eigenvalue vs. jmax for an arbitrary 
sampling of 21 a-configurations. (In all plots for the 6- 
vertex, the integers labeling the (T-configurations have no 
intrinsic meaning.) 



Figure 30: Minimum eigenvalue vs. jmax for a sampling 
of 21 tj-configurations for which all A^'"-* decrease with 

Jmax- 
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Figure 32: A curve for for the smallest eigen- 
value sequence of the 6-vertex, which occurs for 
t7™„ = (2, 0,0, 0,0, 0,0, 0,-2,0), (2, 0,0, 0,0, 0,0, 0,2,0). 

We have used A^""^(jinax) for jmax = |, ■ ■ • ^ ^"'^ 
fitting (blue dots). The red line indicates the eigenvalue 
sequence of ffmim the black line the fit function. The 
error bars for these data points appear in figure 1511 







A (max) 


increasing cr-conf. 


equiv. classes 


58 


30 


total 


323 


8206 


decreasing cr-conf. 


equiv. classes 


3645 




total 


7856 




constant cr-conf. 


equiv. classes 


5 




total 


27 





Table 9: Numbers of equivalence classes of (T-configurations for increasing and decreasing minimum eigenvalue, and 
maximum eigenvalue, for the 6-vertex. The information whether a configuration is increasing, decreasing or constant 
is obtained by comparing the (minimum) eigenvalue for jmax — 2 with that from jmax = 



Figure 31: The 50 smallest eigenvalue sequences (sorted 
with respect to the smallest eigenvalue at jmax = -y) at 
the gauge invariant 6-vertex. The red line indicates the 
smallest eigenvalues from <Jmin in this ordering. 



Figure [311 displays a fit to the largest of the maximum eigenvalue sequences. The fitting function is A^™'''^ (jmax) = 
^o.270±o.026 . (2 . ^■^^^)i.4io8±o.oi2_ Moreover ■= EZi (A^"°"^'^^'*^ (jil) - aJ""'^^(jSx))' = 0.23 for the chosen set 

of 10 values of jmax- 

Table[n]shows the number of equivalence classes of a-configurations for both the decreasing and increasing minimum 
eigenvalues, as in section 1^.3.41 



6.4.4 Number of Eigenvalues 

According to (|4.7p the number of eigenvalues Ncvaisijmax) for configurations with ji < j2 < js < ji < js < je — jmax 
will be given by a polynomial of the order of [jmax]'^^^^ = [jmax]^'- 

N^J:lliM..)=j2^k ■ (Jmax)'^ (6.5) 
fe=0 

with coefficients (and their 95% confidence interval) 
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For this fit the quantity := (^cvais(jmax) - ^evaL (jmax)) takes the numerical value = 3.64 - 10". 



CT-connguraii 



Figure 33: Maximum eigenvalue vs. jmax for the 30 
(T-configuration equivalence classes at the 6- vertex. 



Figure 34: Fit of the logarithm ln[A^'^'^''] of the largest 



((T-configuration index 
(jmax) for 



We have used x[!:^'^^^ 



V 



eigenvalue sequence of figure] 
30) versus ln[2j 

ma 

Jmax — |i ■ • ■ for the fitting (blue dots). The red 
line indicates the eigenvalue sequence, the black line 
the fit function given in the text. 
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Figure 35: Number of eigenvalues at the 6- vertex fitted Figure 36: Log-log plot of the number of eigenvalues 
by a polynomial of degree eight (solid line). at the 6- vertex. The solid line shows a two parameter 

fit to the data. 



As in the case for the 4 and 5-vertex, one can also do a non linear fit to a minimal model N^ilisUma^) = t ■ (jmax)*, 
in order to get the behavior of iVevais in the leading order in jmax- In this case we obtain r = 2.94389 • 10^ ± 85970 
and s = 7.47748 ± 0.01586, (^^evais(jmax) - N^,i2{j^^^)f = 3-49 • lO^^ for the set of 13 data points 

(see figure IM)). 

6.4.5 Cubic 6- Vertex 

Here we include results for a 6- vertex which has linearly dependent edge tangent vectors. The cubic 6- vertex arises 
due to its convenient topology and geometry, and is frequently used in order to construct arbitrarily large periodic 
graphs with definite vertex structure. It is used in particular when working with coherent states, as in [5], as well in 
the recently initiated program [lOj. As it turns out cubic 6- vertex configurations provide a volume gap: the smallest 
non-zero eigenvalue grows with the maximum spin. 

Histogram 

The cumulative histogram in figure [3S] nicely illustrates two main parts of the spectrum: a rising edge for small 
\y ~ 10, which at larger is dominated by the cutoff imposed by finite jmax- 
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Figure 37: Cubic 6- vertex: setup and sign configuration. 
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Figure 38: Overall histograms at the gauge invariant cubic 6- vertex up to j„ 
eigenvalues in all, of which 535,933,540 are zero. iVbins = 2048. 



There are 11,121,868,100 



Extremal Eigenvalues 

As one can see from the plot in figure 1391 the cubic gauge invariant 6- vertex belongs to the class of (T-configurations 
having an increasing smallest non-zero eigenvalue A^™^ as the maximal spin jmax is increased. The oscillatory behavior 
of the smallest non-zero eigenvalue can be traced back to different spin configurations which contribute the smallest 
eigenvalue for fixed je = jmax- We find (except for ji = j2 = js = J4 = js = je = jmax = 5): 
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This nicely illustrates the effect of gauge invariance, namely that the sum of all spins must be an integer in order to 
have the chance to recouple them to resulting angular momentum. 

The plot for the largest eigenvalues is given in figure |40l 

In figures UT] and 221 we show fits of these data. The fitting functions are: 

Here := X^lLi ('^[>'""''''''^**' (jmax) — A^'"^(jmax)^ = 0.02 for the chosen set of 15 data points. For the largest 
eigenvalue sequence we obtain 

^(max)Xfit) ^ _ ^0.859±0.0157j^2j -)1.44954±0.0064 [G 7) 
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Figure 39: Smallest eigenvalues at the gauge invari- Figure 40: Largest eigenvalues at the gauge invariant 



ant cubic 6-vertex up to jmax = ^ 
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Figure 41: Fit of the smallest eigenvalue sequence at 
the gauge invariant cubic 6-vertex up to jmax = ^• 

We have used A^'"^(j,nax) for j,„ax = f, • ■ • ^ for 
the fitting (blue dots). The red line indicates the 
minimum eigenvalue sequence, the black line the fit 
function. 



Figure 42: Fit of the largest eigenvalues at the gauge 
invariant cubic 6-vertex up to j,na> " 



used A^'^''^(j,nax) for jmax = | , . . . for thc fitting 
(blue dots). The red line indicates the maximum 
eigenvalue sequence, the black line the fit function. 



f. We have 



Here := EI-Li (A^"""^'^^'*^(ji?ix) - aJ^"""^ (jH))' = 4.04 for the chosen set of 17 data points. 
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6.5 Gauge Invariant 7- Vertex 

The 7-vertex has 3,079,875 (T-configurations, of which 1,912,373 are 'non-trivial' (not all zero, and not equivalent to 
another up to an overall sign). This, along with the septic number of spin configurations for a given jmaxj makes 
it extremely expensive computationally. Thus we have 7-vertex data only for much smaller spins than for the other 
valences (though simply in terms of number of eigenvalues our 7-vertex data vastly overwhelms that for smaller 
valences!). 

6.5.1 Histograms 
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Figure 43: Overall histograms for the gauge invariant 7-vertex up to jmax = §■ There are 9,475,991,439,360 eigenvalues 
in aU, of which 146,895,794,400 are zero. iVbi„s = 2048. 

In figure[43]we again see no lip in the spectral density close to zero, and furthermore notice that the density appears 
to be linearly increasing for eigenvalues < 3. Both of these features are likely due to the fact that we have computed 
these spectra only up to the relatively small spin jmax = §• At the 5- vertex, for small values of jmax, we observe that 
there is no lip present (e.g. from figure [TU)) . Additionally, at the 6- vertex, for extremely small values of jmax, we also 
fail to observe an exponentially rising edge in the spectral density, but rather see linear growth as above. It seems 
likely, then, that both the exponentially increasing 'edge', and the lip at zero, will arise with larger jmax, as it does at 
smaller valences. 

6.5.2 Fitting of Histograms 

As mentioned above, figure [13] seems to possess a spectral density which rises linearly rather than exponentially with 
Xy. A fit to this rising edge is shown in figure HH in which the logarithm of the number iVgyais of eigenvalues as a 
function of the eigenvalue Xy is given by: 

ln[iVevais(Av.)](^'*^ = (19.72 ± 0.084) + (0.7943 ± 0.067)A^> 

The quantity EfcLi (^ivaiiC^t^^) " ^ovais(A^*'^))^ has the numerical value = 1-25 • lO^^ on the chosen set of 

51 data points. We find that, although the overall shape of the rising edge does not resemble a line in the logarithmic 
plot, the exponential fit is better than a linear fit, at least for the chosen set of data points. 

Figure [33] seems to indicate that we have not yet entered the regime of large jmax- As the exponential growth of 
the spectral density is expected to be a property of the spectrum at large jmax, the exponential fitting is not obviously 
applicable for the computed range of spins. 

6.5.3 Extremal Eigenvalues 

As for the 5 and 6- vertices, by inspection of the eigenvalue data, we find smallest eigenvalue sequences which are 
decreasing, constant, or increasing as jmax is increased. 




A^i'l'(Ay)-0.37.109.e°-™^*- (6.8) 
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Figure 44: Logarithmic plot of the fit to the rising edge at the gauge invariant 7-vertex. The points that were used in 
the fitting process are colored red. 

Since we can go only to a maximum spin of jmax = |j we only have 4 data points (the case where all spins are 
equal to ^ is excluded by the odd valence and the requirement that the sum of spins has to be integer in order to 
obtain gauge invariance), which is too few to identify an overall smallest eigenvalue sequence. However, we can take 
the overall smallest eigenvalue in our numerical data, which arises e.g. from (j-configuration (2 2 2 -4 -2 -2 
224-2 4400 -2)1^ and take this as the smallest eigenvalue g-configuration. The smallest eigenvalues of this 
sequence are contributed by the following spin configurations (here jj — jmax): 
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Table 10: The spin configurations contributing the overall smallest eigenvalue for 
= (2, 2, 2, -4, 0, 0, 0, 0, -2, -2, 2, 2, 4, -2, 4, 4, 0, 0, 0, -2). 

The overall smallest eigenvalue is contributed from the spin config with jj = jmax = |, its numerical value is found 
to be A^'''^|j_^^^^| = 0.000597922. By inspection of the data we find that this eigenvalue is in fact obtained from 
120 different (T-configurations, each with the spin configuration in the last line of table Uni This already indicates 
the importance of understanding the orbits of the (T-configurations under permutations of the edges at a vertex [26] : 
We have partially fixed this symmetry by demanding ordered spins (|2.39p . however as 5 spins are equal in the spin 

configuration leading to ^^^'"^ lj,„^^=| ; we are free to permute 5 out of 7 edges, that is we have 5! = 120 possibilities, 
each giving a different tr-configuration. 



°The sigmas are again ordered in the same pattern as for the smaller valences, i.e. 123 124 125 126 134 135 136 145 146 156 234 235 
236 245 246 256 345 346 356 456. 
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6.6 Numerical Errors 



An interesting feature of our data, which is discussed in section [B. 3. 1[ is a 'lip' in the spectral density of the 5- vertex at 
eigenvalues near zero, as manifested in figure [H One may be concerned whether this is a genuine feature of the theory, 
or if it arises from numerical noise. A look at figure [15] causes alarm, in that it shows the smallest non-zero eigenvalue 
decreasing monotonically to values arbitrarily close to zero. Since the numerical errors are obviously bounded away 
from zero by machine epsilon, for some finite jmax the minimum eigenvalues must inevitably become comparable with 
the numerical noise. 

6.6.1 Error bound provided by LAPACK 

In presenting our data we regard as zero any singular value 

Si < TSsi (6.9) 

where Ssi is given by eqn. (|5.4p and T is a constant numerical threshold. As we reduce T, we will regard smaller and 
smaller eigenvalues as non-zero. For the data presented up to this point T — 1. From section we recall that the 
eigenvalues are related to the singular values by Ay = so in terms of the eigenvalues, our error bound is 

^^^^ 

The Ssi are proportional to the largest singular value of the matrix which gives rise to Si. This is bounded by the 
square of the largest eigenvalues which occur. A more significant effect is that, because of the square root, the errors 
in the eigenvalues are inversely proportional to the eigenvalues themselves. Thus at large jmax, as Ay converges to 
zero, our errors grow without bound. Because of this the computation of volume eigenvalues for valences TV > 4 is 
bounded by machine precision, rather than compute time or system memory. 

The error bars in figure fl31 reveal that, for the 5-vertex, the smallest non-zero eigenvalues are comparable with our 
error bound for jmax ^ ^ ■ Figure US] portrays the closeness of the remainder of our 5-vertex dataset to the LAPACK 




Figure 45: Smallest non-zero eigenvalues of 5-vertex (top, red), along with their LAPACK error bound (bottom, blue). 
The vertical axis is in logscale. 

error bound. The top (red) surface is the minimum eigenvalues for each tj-configuration and jmax- (This is equivalent 
to figure [TSFI ) The bottom (blue) surface is the LAPACK error bound for these minimum eigenvalues. The vertical 

•^'^The (T-configurations of figure l45l are labeled by a sequential index (starting from zero) over the (T-configurations shown in table |4l 
rather than the (T-indicos of figure [TSl 
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axis is in logscale. Note that the errors grow in inverse proportion to the eigenvalues, as expected by (|6.10p . They 
become comparable only for a tiny 'corner' of the dataset, with jmax ^ and ir-indices (from tabled) 1-4. (The 
smallest eigenvalues are always greater than the error bound by construction, c.f. (|6.9p with T = 1.) Zooming into 




Figure 46: Smallest eigenvalues of 5-vertex (top, red), along with their LAPACK error bound (bottom, blue), for 
'large jmax small (i-configuration corner' of dataset. The vertical axis is in logscale. 

this corner, as shown in figure \^ illustrates the effect of the error bound. Consider the (T-configurations labeled 6 
and 7. (These are a-'mdex 3 in tabled) The smallest non-zero eigenvalue decreases steadily until jmax = at which 
point it rises suddenly. For a-configurations 2-4 we see similar behavior, with the smallest non-zero eigenvalue rising 
at jmax = and a-configurations and 1 see a rise at both jmax = ^ and This is likely due to the smallest 
non-zero eigenvalue dipping below the error threshold, such that it is regarded as zero from inequality (j6.9p . Then the 
next-smallest eigenvalue is regarded as zero, resulting in an apparent jump in figure |3S1 Thus we can see that some 
small number of our non-zero eigenvalues, for jmax ^ are descending beneath our error bound. In scction r6.6.3l we 
will estimate the number of eigenvalues that are thus 'lost'. 

6.6.2 Characterization of true numerical noise 

How well does the error bound provided by LAPACK describe our numerical noise? To explore this question, we can 
vary T in (j6.9p . and watch what happens to our data. 

The simplest measure of the effect of varying the error threshold is to count the number of zero eigenvalues. 
Figure HT] plots the overall number of zero eigenvalues as a function of T. The left figure shows that the number 
of zeros is essentially constant for any T > 0.1. Below that value the number of zeros drops off rapidly, indicating 
that a significant number of zero eigenvalues are being regarded as non-zero due to the threshold being less than the 
numerical noise. The plot on the right shows a relatively slow increase in the number of zeros for larger thresholds, 
which is expected. Here genuine non-zero eigenvalues are being regarded as zero due to the large error threshold. 

Given the results of the above paragraph, we may be tempted to use T = 0.1, so as to catch as many genuine 
small but non-zero eigenvalues as possible. Figure l48l shows what happens to figure [45] when T is reduced to 0.1. 
Two features are of note. In our 'large jmax small tr-configuration corner' we see that with T < 1 the error now can 
exceed the smallest non-zero eigenvalues, which it does. We also note a huge spike for tr-configuration 72 (tr-index 36), 
jmax = for which the smallest non-zero eigenvalue and its error jump to meet (in fact cross) at around 10~^. 

This occurs because numeric noise has exceeded the error threshold, and is now being regarded as a genuine smallest 
non-zero eigenvalue. Note that this essentially random event can occur at any point in our dataset. We thus see that 
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Figure 47: Number of zero eigenvalues for 5-vertex, as a function of error threshold T. The green line at the bottom 
of the left plot indicates the number of 'obvious' zeros, with T = 0. The plot on the right zooms in on the data points 
for T > 0.08. 




Figure 48: Smallest eigenvalues of 5-vertex (top, red), along with their LAPACK error bound (bottom, blue), for 
T = 0.1. The vertical axis is in logscale. 
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our true numerical error can exceed one tenth of LAPACK's bound, and thus stick with their reported bound with 
T = 1. The fact that this problem occurs just below the reported bound testifies to it being an accurate and tight 
estimate of the true numerical error. 

At T = 0.01 the analog of figure H51 becomes mostly noise: The smallest non-zero eigenvalues are a constant ~ 10~^ 
for most (T-configurations and values of j'max- 

6.6.3 Effect of noise on results 

Reducing T allows us to probe the 'problematic corner' a bit closer to the numerical noise. Figure H^l portrays the 
effect of smaller T on the smallest non-zero eigenvalues. Clearly some valid non-zero eigenvalues are being washed out 
due to the noise. 
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Figure 49: Smallest eigenvalues of 5-vertex, for T = 1 (top surface: upper side blue, bottom side magenta), and 
T = 0.1 (bottom surface: upper side red, bottom side green). The vertical axis is in logscale. 

The number of eigenvalues 'between the two surfaces' of figure HH is the difference between the two leftmost data 
points on the right hand plot of figure Wl\ which is 413,160. The largest difference between data points on the right 
hand plot is 8,230,920. It is impossible to determine how many more non-zero eigenvalues there are below the bottom 
surface, but we expect it to be approximately of this larger order of magnitude. Now the bin width of the histogram 
in figure[n]is ^^''"2048^^" ~ 0.242785. Clearly this is much larger than any of the eigenvalues close to the error threshold 
that we are considering, so the only effect of noise on the 5-vertex histogram is in the bin for smallest Ay. The 
population of that bin, for T = 1, is 154,430,880. The true value is slightly larger than this. An addition of 8,230,920 
to this number will have an indiscernible effect. In any event the lip is a genuine feature of the volume operator, and 
not a result of numerical error. 

6.6.4 Valences N 

Figure [50] depicts the distance between the error bound and smallest non-zero eigenvalues for the 6- vertex. The minimal 
eigenvalue is closest to the error bound for the a-configurations 7751 & 7752 shown, out of 8207 a-configurations in 
total. Figure ISTl shows the same numbers for (j-configuration 7751, with error bars. This curve also appears as one of 
the fifty shown in figure [311 It appears that we are just at the threshold of the numerical noise for the 6- vertex. Given 
the analysis of the 5-vertex errors above, we expect the errors to have negligible effect on the 6-vertex results. 

(One of) the smallest eigenvalue(s) for our 7- vertex data, as described in section [B. 5. 31 arises from a-configuration 
(2 2 2 -4 -2 -2 2 2 4 -2 4 4 -2), with spins (2 2 2 2 2 3 5). It's value is 0.000597922, with error bound 
2.54519 • 10^^. Clearly numerical error is not a significant issue for the 7-vertex. 
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Figure 50: Some of the smallest non-zero eigenvalues of the 6- vertex (top, red), along with their LAPACK error bound 
(bottom, blue). The vertical axis is in logscale. 
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Figure 51: Minimal eigenvalues for (T-configuration 7751 = (0, 2, -2, 0, 0, 0, -2, 0, 2, 0), which yields the overall 
smallest non-zero eigenvalue at the 6-vertex up to jmax = ^ ■ ^ = 1 . 
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Since the smallest non-zero eigenvalues for the 4- vertex rise with jmax, we do not expect any effect from numerical 
noise. The relative error for every eigenvalue at the 4-vertex will be negligible. The same situation holds for the cubic 
6- vertex of section 16.4.51 

6.6.5 Symmetric Eigenproblem 

Earlier we computed eigenvalues more directly using subroutines LAPACK provides for the 'symmetric eigenproblem'. 
To cast the matrices into symmetric form we computed eigenvalues for QQ^ rather than Q. In this case numerical 
error is much more prominent, for a number of reasons. 

For the symmetric eigenproblem, LAPACK employs the 'relatively robust representation algorithm', which involves 
an L£)L"^ factorization and a number of 'translates', to return the matrix eigenvalues with a fixed absolute error 

^e, = e||QQt||^ (6.11) 

where HQQ^Hi is the 1-norm (maximum column sum) of QQ^ . This 1-norm is bounded from below by the largest 
eigenvalue. For singular value decomposition (SVD) the relative error is fixed and thus proportional to the largest 
eigenvalue. 

We see that the relative error resulting from any algorithm is bounded below by the machine precision. Thus the 
error in the smallest eigenvalues can never be less than the maximum eigenvalue times machine epsilon, so that, when 
considering numerical error, what is important is the ratio between the largest and smallest eigenvalues which arise 
from a given matrix. When we square the matrix, we double this ratio, which in effect halves the relative precision of 
the computation. 

Note that, as manifested in figures HSl and [50 1 the error bound in the smallest eigenvalues appears to be growing 
exponentially with jmax- This occurs simply because we are taking the square root of the singular values to get 
the eigenvalues, and the minimum eigenvalues decay exponentially toward zero (as one can see from figure fT7|). so, 
according to can. 16.101 the error in the square root must grow exponentially. In the case of the symmetric eigenproblem, 
since we square Q before computing eigenvalues, we must take a 4th root to get the volume eigenvalues. In this case 
the analog of (j6.10p has the volume eigenvalue cubed in the denominator. Thus the log of the errors in the volume 
eigenvalues will grow three times as fast as for SVD. This effect, along with the effect of squaring Q discussed in the 
previous paragraph, represents a considerable increase in numerical error. We are thus able to explore the spectrum 
to much higher values of jmax with the singular value decomposition. 

6.7 Histograms of Entire Data Set 



5e+11 



4.5e+11 



"T" 



"T" 



4- vortex 
4, 5- vertices 
4, 5, 6- vertices 
4, 5, 6, 7- vertices 



4e+11 



3.5e+11 



3e+11 



2.5e+11 



2e+11 



1.5e+11 
1e+11 



5e+10 



200 



400 



600 



800 



1000 



1200 



1400 



1600 



1800 



Figure 52: Histogram of entire dataset. A^bins — 4096. 



Figure [52] presents four histograms of our entire data set. Each represents all eigenvalues we have computed up to 
a given valence. (Thus the blue curve, for example, includes all 4-, 5-, and 6-vertex eigenvalues we have computed. 
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Note that only three of these are new. The histogram for the 4- vertex already appears as the top curve in figure [51) 
There are 21,551,022,251,376 eigenvalues in all, of which 827,849,074,512 are (regarded as, c.f. discussion on errors in 
section [6.61 above) zero. 
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Figure 53: 'Central region' of histogram of entire dataset. 

Figure [S3] zooms in on the 'central' region of figure We see that the lip of the 5-vertex histogram gets washed 
out by the higher valence data. This is likely an artifact of the fact that only quite small values of spin are considered 
for the larger valences. We do note that in any event the spectral density at zero increases with valence. 

Figure [SD displays the tail of figure [221 Here all four curves overlap, indicating that these large eigenvalues come 
from the 4- vertex, for large spins. 
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Figure 54: Tail of overall histogram. 

7 Analytical Results at the Gauge Invariant 4- Vert ex 

In this section we present three analytical results on the spectral properties of the volume operator at the gauge 
invariant 4-valent vertex. We first prove the spectral simplicity and second give a general formula for the eigenvectors 
of the volume operator at the 4- vertex depending on the eigenvalues, which agrees with the earlier publication [17j . 
Third we analyze how the smallest non-zero eigenvalue scales with the maximal spin j,nax at the vertex. 

For the methods we are using here the particular shape of the operator, when represented as a matrix, will be 
crucial. In fact many theorems developed for Jacobi matricesEl can be carried over to the particular case under 
investigation. Here [27] provides the crucial theorems and proofs. 

7.1 Setup: Explicit Expression for the Volume Operator as D-dimensional Matrix Qd 

The general case of the gauge invariant 4- vertex deals with 4 edges ei, . . . , 64 outgoing from the vertex v. Each edge ex 
{K = 1 ... 4) carries a 2jK + 1 dimensional representation of SU{2). Due to gauge invariance the standard recoupling 
basis states (|B.17P are given in this case as 



a2 



(J1J2) a3(a2 j3)=j4 J(a3j4) = M = ) 



(7.1) 



where the intermediate recoupling 03(02^3) has to equal J4 due to the Clebsch-Gordan Theorem and as a result the 
intermediate recoupling 02(^1 ^2) is the only degree of freedom for fixed spins ji, . . . Therefore expression (j4.12p 
for the matrix representation of the volume operator simplifies dramatically to give 



Z-a{123) 9123I 



Z|2. [a(123)]2.QtQ 



(7.2) 



where we have introduced the shorthand Q := 5123. Note that (t(123) = e(123) - e(124) + e(134) - e(234) = 0,±2,±4 
gives a constant numerical prefactor depending on the relative orientations of the 4 edges only. In the following we 
will assume that cr(123) 7^ but leave its numerical value unspecified, that is we will drop it from our formulae and 
reinsert it at the end. We set the regularization prefactor Z = i, c.f. section Using p.9p together with the 

definitions of section 13.3.11 we obtain for the matrix Q = Qd [D denotes the dimension) the following shape [131 [M] '■ 



^^That is symmetric tridiagonal matrices. 

^® Note that due to the definition of the eigenvalue of V as := ^|Aq|, this docs not change the volume spectrum. Because Q is 
an antisymmetric matrix, setting its matrix elements to be real or imaginary only rotates its spectrum in the complex plane either to the 
complex axis or the real axis. 
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■•• 


-qo-b 





























qo-b 





-9-D-4 





























9-D-4 





-qD-3 





























qD-3 





-qD-2 





























qD-2 





-qo 


V 


























qo-i 






(7.3) 



where the matrix elements are given by 

9fc = 



^(2a2 + l)(2a2-l) 



(jl + j2 + 02 + l)(-jl + j2 + a2)(jl - j2 + 02) (jl + j2 
X (j3 + 34 + 02 + l)(-j3 + j4 + a2)(j3 - ^4 + a2)(j3 + J4 



02 
0,2 



1)X 
1) 



and the intermediate recoupHng step 02 = 02(^1 ^2) is given as 02 = 



(min) 



k, that is 



fl ■ • I I • ■ ll (min) ^ ^ (max) . r . . . . ] 

max [|j2 - b4 - JslJ = a2 < 02 < - mm [ji + j2, J3 + J4J 



(max) 



(7.4) 



(7.5) 



and D = dim Qi 



n ^ 



a. 



{min) 



7.1.1 Known Results 

1. Spectrum: By choosing Z ^ i, (j7.3p is an antisymmetric matrix with purely imaginary matrix elements. Its 



eigenvalues A are thus real and come in pairs A = ±|A|. Moreover its eigenvectors '^^ = (^'^, 
orthogonal and there exists a unitary matrix 



Ud = 



*2^ *2 



*^)^ are 



(7.6) 



such that Ud^^d = QdUd, = U^'^ Qd Ud = uj^ Qd Ud, = diag(Ai, . . . , A^) 

2. 0-eigenvalues: It has been found in 24J that as a single eigenvalue is contained only in the spectrum of odd 
dimensional matrices Qd- 

7.2 Simplicity of the Spectrum of the Matrix Qfj 

By inspection an eigenvector — {'i'^, ^^''^ '^'^ 



for the eigenvalue A of with {Qd - Al)^'^ = wiU fulfiU: 



(I) 




= 


(11) 




= 


(III) 




= 



Kk <D 



(7.7) 



Now suppose A 7^ Ocj Setting ~ leads, together with (I) and the fact that 7^ for fc = 1, ... TV — 1, to 
^2 = 0, which then in turn implies together with (II) that = Vfc. Therefore in order to get a non-vanishing 

! 

eigenvector ^f-^, A 7^ we must have 7^ 0. 

Now suppose there exists an eigenvalue A 7^ of Q^i that has multiplicity > 1 : Let vj/-^^ vj;-^ be two eigenvectors for 
A. Then we can find two constants a,a G (D, a,a ^ 0, such that we can construct an eigenvector — a^'^ + 
with the property <I>^ = 0. As we have seen this implies $^ ^ Vfc, and thus = a'^^ + ci^'^, that is ^f^, are 
linearly dependent and therefore there exists only one eigenvector ^''^ for each eigenvalue X ^ oi Qd- We have thus 
proved the following 



^'^The case A = is already discussed in | 24| . 
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Theorem 7.1 Spectral Simplicity of Qd 

(i) The spectrum spec^Qo) 7^ is simple (consists of D distinct real numbers) 

(a) If^^ = {^^, . . . , *^)^ is an eigenvector of Qd (Qd'^^ = A^*^, A 7^ OJ then 7^ 0, ^ 

7.3 Eigenvectors of the Matrix Qd 

Starting from (I) in (j7.7p by setting = x = const one can explicitly construct the components One gets, in 
agreement with [iTj (using the integer number L> 1 with L < ^ for even D and L < for odd D): 



2L-1 

n gfc 



L-1 2L-2 fcAf-2 fei-1-2 fc3-2 fe2-2 

M=l fejv/=2i>-f-l *:m-1=2(A/-1)-1 fc,=2(Af-!)-l '=2=3 fci=l 



'2L+1 - 2L 

n Qk 



+ E ^ 



L 2L-1 
2{L-M) 



fcA^-2 fci-l-2 fc3-2 fc2-2 

E E ■■• 

fcM^2M-l kM-i^2(M-l)-l ki=2{M-l)-l ^2^3 fci^l 



Ev V 2 2 2 2 2 



(7.8) 



where L is a positive integer < ^ . When M — \ the upper bound of the fci sum is that of the /cm sum. One may 
explicitly check that these states fulfill (II) in ()7.7p . In order to fulfill (III) we find the conditions 



D = 2L, even 



2L-1 

n 

fc=i 



fc3-2 fc2-2 



L 2L-1 

,2.^^,2(.-M) ^ ^ ^ E«.-. 

M=l kM=2M-l fcM-i=2(Af-l)-l fe2=3 fci = l 



■9fe2 9fci 



(7.9) 



D = 2i + 1, odd 


„ 1 a;A 


A2^ + 




n Ik 





L 2L 



Af=l 



X! X! 'ikMikM-i 

kM='2M-l fcA_f_i=2(Af-l)-l fc2=3 fci = l 



■9fe2 9fei 



Thus, in order for (j7.9p to be satisfied, A has to be a root of the characteristic polynomial 

L 



7r2L(A) 



M=l 

^2L+i(A) = A[A2^+E&i/"'A^^^-''^ 
with coefficients 



M=l 



fcn-2 



fcM-2 



r,(feo) 



E E 



D = 2L even fco = 2L + 1 = D + 1 

D = 2L + 1 odd fco = 2L + 2 = D + 1 



fci-l-2 A;3-2fc2-2 

E\ ^ \ " 2 2 2 2 2 



(7.10) 



(7.11) 



kM=2M-l feM_i=2(Af-l)-l fci=2(J\/-i)-l fe2=3fci = l 



Since the bracketed terms in the characteristic polynomials (|7.10p contain only even powers of A we may replace 
^2fc _ ^k Tfj-^iY^ ~ |Ap in order to arrive at a reduced purely real notation for them: 

= A^^ + E h^M^X^^"^-^'^ 

M=l 



(7.12) 



Af=l 



where we have pulled out a prefactor (—1)^^ coming from the purely imaginary nature of the matrix elements (17. 4p . 
7.3.1 Traces of (Q|,)" 

Knowing the coefficients ()7.1ip of the characteristic polynomial, together with Newton's equation valid for an arbitrary 
matrix A with tr[A"] = s„ 



n-l 

(ko) , j,(fco) 



Ec; SI 



1=1 



n> 1 



(7.13) 
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^ D-1 

we are now enabled to successively calculate the traces of all powers of starting from si — —2 ^ q^. This 

fc=i 

is of special relevance because it allows us to reconstruct the spectral density function p(A) of eigenvalues A of the 
hermitian operator Qo in question (see e.g. |33| . chapter 2). 
p(A) fulfills the properties of a probability distribution 

(i) > V A positivity 



+ 00 

(ii) J p(A) = 1 normalization 



— oo 

D ^ 

(iii) p{X) = J2 ^{^ ^ ^k) for a discrete spectrum, as is the case for Qu 

k=l 

Then 



II- n 

tr[(Oi,)"] = ^(Afc)" = TO„ = / A" p{\) dX (7.14) 

k=l •{„ 



is called the n*'^ moment of the distribution p(A). The Fourier transform of p(A) 



X{r) = J e-'^^ p(A) dX (7.15) 

— OO 

is called the characteristic function. By inversion of the Fourier transform p{X) can be reconstructed from x('^)- 



+ 00 

1 



P(A) = ^ y e'^^ xir) dr (7.16) 

— oo 

Now by expanding the exponent in (j7.15p and inserting (|7.14D . the Fourier transform xi^) of the eigenvalue distribution 
p(A) can be expressed in terms of the momenta (|7.14p as 



i-ir 



n=Q 



A" p(A) dX 



Y^^-r-m^ (7.17) 



n 

n=0 



Inserting this back into (|7.14p one finally obtains p(A) in terms of the momenta. However this formal way of obtaining 
p(A) will not be directly accessible in most cases, since the series (|7.17p is not convergent, and one needs a closed 
expression for x{r)j as otherwise inserting the formal power series (|7.17p into ()7.14p will result in ill-defined divergent 
expressions. One then has to put in by hand an assumption on the kind of distribution one is working with (see e.g. 
|19j). A note on the results obtained in [TH] is given in the companion paper There are also numerical procedures 
to reconstruct a statistical distribution from a finite number of its moments using spline interpolation [20, . 



7.4 Upper and Lower Bounds on the Eigenvalues of 
7.4.1 Upper Bound 

One can find upper bounds for the eigenvalues by applying the theorem of Gersgorin (see e.g. [31] . [34]): 

Theorem 7.2 Gersgorin 

Every eigenvalue X of an (n x n)-matrix Q lies at least in one of the discs 

n 

|fe-A| <Y\q^j\ i = l,...,n 

3 = 1 

That is every eigenvalue lies in a disc centered at the diagonal element qu with radius the sum of moduli of the 
off-diagonal-elements g^j, i ^ j oi the i"^ row or column. In the case of the gauge invariant 4- vertex this theorem 
simplifies, due to the banded matrix structure of (|7.3p and the fact that qu = 0, to 

|A| <^|gy| = |g, ,_i| + |(7. ,+i| (7.18) 

A general upper bound was found already in [24 by directly employing theorem 17.21 
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7.4.2 Lower Bound 



Finding a lower bound on the modulus of the smallest non-zero eigenvalue turns out to be a much harder task in [21] . 
The idea is to invert the matrix Qd and find an upper bound on the eigenvalues of the inverse Q^ , which then in 
turn would serve as a lower bound on the spectrum oi Qu- The problem is that in the case of odd dimensional Qd 
cannot directly invert the matrix but have to project out the null space first. That results in a messy matrix which 
has lost some of the banded structure of the original Qd given in (|7.3p . 

We will show now that this procedure is not necessary. In order to see this consider the D) component of the 
resolvent of Q^: 



(7.19) 



1 




Qd- AId 


DD dct[QD-i - A1d-i] 7rD-i(A) 




DD det 


Qd - AId 


det 


Qd - AId 


7rD(A) 



where we have used the general definition of the matrix element of the inverse matrix [Q — Qkj obtained 

by deleting row k and column j from Q and taking the determinant of the remaining submatrix. Moreover 7rD(A) 
denotes the characteristic polynomial of Qd, 7rD-i(A) is the characteristic polynomial of the D — 1 dimensional sub- 
matrix Qd-1 one obtains by deleting row D and column D of Qd- In what follows we will drop the index D from 
Qd, Ud etc. for clarity. 

From (|7.6p we know that Q = UAU\ A — diag(Ai, . . . , Ad), and U^U — UU^ = 1, U being a unitary matrix. 
Therefore we know that 



elm 



= J2 Ulr Ks C/L 

r.s 

= Ulr [Xr ■ Srs] 

r.s 

= J2XrUlrUL 



(7.20) 



and, since U diagonalizcs (Q — Al) ^ as well: 



Y^Uir [(Q-Al)-!]^^ Ul 

r,s 

Y^Uir [(A-Al)-i]^^[/L 

r,s 
r,s 



(7.21) 



Thus we find using (|7.6p : 

[iQ-Xl)-']DD = 



7rD-i(A) 



7rD(A) 

where the overline denotes complex conjugation. 



Xr-X 



E 



A,. - A 



E 



At- a 



(7.22) 



Now from theorem 17. II we know that we must have J 7^ Vr = 1 . . . Z?. Therefore we know that 



V{X) = 



7rp-i(A) 
7rp(A) 



\ " I ^p I 

^ Xr-X 



(7.23) 



is a real valued function and its derivative 



|^nA)^E^^>OVA 



(A. - xy 



(7.24) 



That is, V{X) is monotonic increasing V A S R, has poles at every X = Xr, r — 1 . . . D, and is finite V A 7^ A,.. So 
we conclude that T^(A) has the form illustrated in figure [m The function V{X) has precisely one root in the interval 
[Afc,AA;+i]. This means that the roots of 7rp_i(A), 7rp(A) interlace: in between two eigenvalues Xk^Xk+i of Qp there 



70 



V{X) 



^feyf 1 



D odd 



Figure 55: The function V{X). 



D even 



Qd-1 

Qd~-2 



Qd 

Qd-1 

Qd-2 



Figure 56: Interlacing of the spectra for Figure 57: Interlacing of the spectra for 
odd dimension D. even dimension D. 



is precisely one eigenvalue of its submatrix Qu^i. From this we conclude that the structure for the spectra oi Qd 
and its submatrices Qd-i,Qd-2, ... is as illustrated in figures 15^ and 1571 The D odd case in figure EH shows that a 
lower bound for the smallest non-zero eigenvalues of Qd is given by the smallest eigenvalue of the even dimensional 
submatrix Qd-i- However, as mentioned above, one can show [24j that even dimensional matrices of type of Qd have 
no zero eigenvalue in their spectrum and thus it seems feasible to give a lower bound for their eigenvalues by finding 
an upper bound for the eigenvalues of its inverse matrix Q^ . 

N-l 



Set L(*'^^^) = 



1 



n 



q2i 



for M < N and i^^^'*^) = 



1 



92A/-1 ^-^Jj 121+1 



<l2M-l 



such that e.g. = L^^'^) ^ 



or 



^(2^2) ^ J_ ^(2,3) ^ _24_ ^(2,4) 

' 9395 ' 



Then if is as given in (|7.3|) and D is even we have for its inverse (using 



n 



Q 



D 



( 





L(i'2) 








l^(l,n-2) 





^(l,n-l) 

















• 



























L(2.2) 


L(2.3). 





2,(2,«-2) 





2,(2,n-l) 





_L(2,n) 




_L(i^2) 





_L(2,2) 


• 






























L(3,3). 





2,(3,«-2) 





i(3,n-l) 





i(3,«) 




_i(1.3) 





_i(2.3) 


_i(3.3) . 






























• • 





2^(n-2,n-2) 





2^(«-2,«-l) 





j^(n-2,n) 




.L(1,"-2) 


- 


_L(2.«-2) _ 


.^(3,«-2) . . 


_7^("-2,«- 


2) 

























• • 











^(n-l,ri-l) 





^(n-l,ri) 




,2,(1."-!) 


- 


-i(2,n-l) - 




j^{n — 2,n — 


1) - 




1) 



















• • 



















V 







-L(2,n) 


_L(3,n) • • 


j~^{n— 2.71 


) 




) - 


L(" 


n) / 



(7.25) 
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We are now left with the task of finding an upper bound for the row/column sums of Q in order to apply theorem 
17.21 For this we have to use the explicit form of the matrix elements (|7.4p . Let us introduce the following shorthands: 



(m.n) _ 
fc 

(m,n) 



a. 



(min) 
2 

(min) 



k+l 



iim + in + a^"™'' + k+ l)(-jm + jn + O^"™^ + k){j 



, (min) 
■ Jn + a2 



k) 



(7.26) 



which enable us to write for 

N-l 

92M-1 



M < N 



1 



q2i 



l=M 



N-1 



= [2(0^"^ + 2M - 1) - 1] [2(ar"' + 2M - 1) + l] H 



[2( 



{min) 



+ 2/ + 1) - 1] [2(4™") + 2^ + 1) + 1] 



[2(4™^") + 21) ~ 1] [2{a: 



(min) 
2 



N-1 (1,2) ..(3,4) 



1 



..(1,2) ..(3,4) 11 ..(1,2) ..(3,4) 
'-'2M-1 ^^2^1-1 l=M '■^21+1 "^2i + l 



7?' 

rC. 



(1,2) p(3,4) 



2M- 



W-l p(l,2) p(3,4) 

n ^2l ^21 
n(l,2) p(3,4) 



1 -^2M-l i=A/ -^2; + l -^21+1 



21) + 1] 



2a 



(mm) 



4M - 3) (2a 



(min) 
2 



4iV- 1 



..(1,2) ..(3,4) n(l,2) n(3,4) 
'^2M-1^2M-l-^2Af-l-^2Af-l 



where in the last line we have introduced 



N-l p(l,2)p(3,4)..(l,2)..(3,4) 

n ^2l ^21 ^2i ^2i 
p(l,2)p(3,4)..(l,2)..(3,4) 
l=M ^21-1^21-1'-^ 21+1 "-^ 21+1 



x{2l) 



(7.27) 



x{2l) 



n(l,2)n(3,4)..(l,2)..(3,4) 

^2i ^2; ^2; ^2i 

p(l,2)p(3,4)..(l,2)..(3,4) 
-n-2i-1^2/-l'^2; + l ^2i+l 



(7.28) 



Obviously, due to positivity of all terms contained in [/^"'-"^ -v^e always have 

Let us take a closer look at the product terms contained in (|7.27p . 

Using x^™'") = jm + Jn — 02™*") + 1 we may write 



X <l 



N- 



n 



1 7-,(m,n) 
^21 



J. J. n(m,n) 



-'^ 1 ^(m,n) 21 



n 

l=M 

-N + 1 

n 

l = -M 



r^(m.n) -21 + 1 

2,(m,n) 2Z 
,(m,n) +2^+1 



^ *^ - 27V + 2r 



n 



[;(m,n) - 27V+ 2r + 1 



introduce: r = I + N 



r=l 

r(i + - M) r(f + ^^~ TV) 



r(i + 2±i2> _ AT) r(| + ^ 



M) 



(7.29) 
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Using 



(m,n) 



I ■ I (min) -I 



(mi7i) 



(7n,n) 



Jm Jn ^ 



(min) 



we can write 



N-l jj{m,n) N~l 

n 



l=M ^21+1 



n 

l=M 
3 

3 

= n 

3 

= n 

3 r 

3 



x^ + 2l + l 



n 

n . 
n 

N-A 

n 



2/ 



introduce: s = / — M + 1 



Xf, + 2M + 2{s - 1) 
f 2M + 2(s- 1) + 1 

+ 2M + 2s 



s=0 

m-M 



2M + 2S + 1 



Xf, + 2M + 2s 



+ 2jV + 1 



Xf, + 2M 



Xn 



2M + 2s + 1 



. s=i 

T{l + ^ + N) r(| + ^ + M) 



Xf, + 2N 



2M+1 



2N +1 



2M 



N) 



x^ + 2N Xf, + 2M+l 



r(i 



M) r(j 



+ iV) 



(7.30) 



since ^|^^_^^jv^ x^+2A4+i — ^ because M < N and the function f{y) — is monotonically increasing. Using the 
expansion of the F-function 28J one may check that 



'x'^-^V^i^ < r(x) < 2-<B-''x'=-iV2^ 



moreover we have 



< (E < 2 • 1 



V a; e R a; > 1 



V a; e R x> I 



and we can estimate 



nx) 



T{x- 



< 



2V2 



and 



r(a;- 



r{x) 



< 2y/^ 



and finally obtain from (j7.29p 



p(m,ri) 
l=M ^21-1 



r(i 



M) r(| + 



N) 



For ((73011 we find 

n 



N-l jj{m,n) 
^21 
(m.n) 



< 



l=M ^21+1 

Now recall that 



n 



r(f + 2;^-Af) r(| + ^2^-M) 



< iV2 



1 + 



N 



M 



(7.31) 



(7.32) 



(7.33) 



(7.34) 



r(i 



M) r(| + ^iv^ + iv) 



< (4%/2)3 n 



1 + 



1 + 



■N 



ji < i2 < ja < ji 

where J4 equals the maximal spin jmax- Note the ranges of the intermediate recoupling spin 02 = a2 
to ()7.5p . Now without loss of generality assume 



(7.35) 



(7.36) 



{rnin) 



k according 



by 



{min) 



f3 Then the dimension of the matrix is given 

(7.37) 



airnQD = ^2 — 0-2 + 1 = 2ji + 1 
which imphcs M < N < ji + ^. So finally by using the definitions for x^"^'^"* above with 



(mm) 



J2 - Jl 



a:(i'2) ^ 2ji + 1 a;(3^4) ^ ^-^ ^ ^-^ _ +J1 + I 



2j2 



4'''^ = 2(J2-Ji) 



''1 
''2 

J3,4) 



i3 + 34 + 32 - 31 + 1 
- j3 + k + j2 - jl 
35 - j4 + j2 - 3l 



(7.38) 



31 Note that we might equivalently choose Oj"'*"'' = 34, — js here, as 117.2711 is symmetric with respect to interchanging (j'l, ^2) (jstji)- 
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we can estimate using (|7.34p 



N-l „(1,2) p(3,4) 

p(l,2) p(3,4) ^ 
l=M ^21-1 ^2l~l 



< 32 



moreover with ((7351) : 



(2 + - 27V) (2 + xt^^^) 27V) 



(2 + a:(i^2) _ 2M) (2 + a;(3>4) - 2M) 
(3 + 2ji - 27V) (3 + j3 + J4 - J2 + ]i - 2N) 



(3 + 2ji - 2M) (3 + j3 + u - h + ii - 2M) 



(7.39) 



N-l (1,2) ..(3,4) 3 



< (32)= 



(2 + 4^^^^ + 2M) (2 + xf^^ + 2M) 



(2 + xl,'''^ + 2iV) (2 + ^■''"^ + 2iV) _ 



(3,4) 



(3 + 2j2 + 2Af ) 
(3 + 2j2 + 2iV) 

(l+i2-Ji + M) 



(l+j2-Ji+iV) 



(3 + j3+i4+j2- ji + 2M) 

(3 + j3+j4+j2-jl+27V) 

(2-j3+j4+j2-Ji + 2Af) 



(2-.?3+j4+j2-.n+27V) 



(1 + M) 



(1 + 7V) 



(2 + j3-j4+j2-Ji + 2M) 
(2 + j3-i4+j2-ji+2iV) 



(7.40) 



Combining (|7.39p and (|7.40p we finally get an upper bound for the last product on the right hand side of (|7.27p by 
taking into account that < Af < iV: 



JV-l N-l (1,2) p(3,4) ..(1,2) ..(3,4) 

n-rOn - TT -"2; ^21 ^2i ^2i 
11 n(l,2) p(3,4) ..(1,2) ..(3,4) 
;=M l=M ^21-1 ^21-1 ^21+1 ^21 + 1 



< (32)^ 



(3 + 2ji ~ 2jV) 
(3 + 2ji - 2A'/) 

(3 + 2j2 + 2M) 



(3 + 2j2 + 27V) 



(3 + i3+j4-j2+ji-2A^) 

(3 + j3+j4-j2+jl -2Af) 

(3+j3+j4+j2-Jl+2Af) 
(3 + j3+i4+j2- jl +27V) 



(1+J2- ji+M) 

(l+j2- jl +iV) 



(2 - js + 34 + 32 - 31 H 


-2Af)' 


1 
2 


[(n 


- A/)' 


1 
2 


(2 + j3 - 34 + 32 - ji H 


h2Af) 


(2 - J3 + 34 + 32 - ji - 


^2N)_ 




[(1- 






. (2 + js - j4 + j2 - jl - 


f 27V) 



(7.41) 



By construction we have already seen that 



N-l 



n ^(20 < 1 



l=M 



(7.42) 



however the reason for the explicit calculation of (|7.4ip will become clear in a moment. Using (|7.42p we give a first 
upper bound on \L^^-^)\^ in ^^lYh 



j^(M,N) 



2 (24™")+4Af-3)(24™") 



<- 



4iV- 1 



..(1,2) ..(3,4) p(l,2) p(3,4) 
'^2Af-l"^2M-l-^2Ar-l-^2Ar-l 



(2a^"") + 4A^ - 3) (2a^™"^ + 4iV - l) 



{31 + 32 + a^""^ + 2A^)(-ji + J2 + af"^ + 2Af - l)(ji - + 4""'^ + 2Af - 1) 

V ' 

Ti 

1 

(J3 + J4 + 4"""^ + 2Af )(-j3 + 34 + a^"^") + 2A^f - l)(j3 - 34 + 4"""^ + 2M - 1) 

" V ' 

To 



(31 + 32 a^"^ - 2iV + 2)(j3 + 34 a^""") - 27V + 2) 



(7.43) 



T4 
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In order to get an estimate for the smallest non-zero eigenvalue we will now perform a power counting in = jmax 
for all the different factors Ti . . .T4 in (|7.43p . that is while respecting (|7.36p we have to minimize the obvious inverse 
power of j'4 ~ jmax in (|7.43p . To do so we make the following ansatz: 



(7.44) 



Jl 


— Jmax 




M 


— CtM jmax " 






= Ck2 jmax 




N 


= CtN jmax H 




is 


— Jmax 


+ /33 








J4 


— Jmax 











which is justified if in the factor terms Ti,T2^Ti contained in (|7.43p the power of jmax should be lowered (raised 
in T3) while simultaneously respecting (|7.36p as jmax ^ 00 [3- Moreover by inspection we find that in (|7.4ip only 
M, N matter (and M < N) , since all other factors in the numerator and denominator are identical and it is therefore 
sufficient to look at (|7.4ip alone. 

Upon inserting (|7.44p into (|7.43p and again assuming aj™*"^ — j2 — ji 



one obtains 



j^(M,N) 



<- 



(2j2 - 2ji + 4M - 3) {2j2 - 2ji + AN - l) 



(ji + J2 + .h - 31 + 2M)(-ji + j2 + k - Ji + 2M - l)(ji - ]2 + j2 - 31 + 2Af - 1) 

1 

(ja + J4 + j2 - ]1 + 2M){-J3 + J4 + j2 - Jl +2M- l)(j3 - J4 + j2 - Jl + 2Af - 1) 

1 



(jl + .72 - j2 + Jl -2N + 2) (.73 + j4 - .72 + Jl - 2iV + 2) 



< 



2 [{a2 -ai + 2aM) ■ ji + /32 - /3i + 2/3m - |] 2 [(02 - ai + 2aAr) • j4 + /32 - /3i + 2/3, 



5Af 



2 [ci • j4 + /32 + /3Af ] 2 [C2 • j4 + /32 - /3l + PM 

1 



i]2[c3 •j4 + /3Af - ^] 



[c4-j4 + /33 + /?2-/3l + 2/3M- 

1 



[C5 • J4 - /33 + /32 - /3l + 2/3Af - 1 ^ ] [C6 • j4 + /33 + /32 - /3i + 2/3m - 1 ^ ] 



2 [C7 • j4 - /?1 - + 1] [C8 • j4 + /93 - /32 + /3i - 2/3iV + 2] 

where we have introduced the shorthands 



Ci = 0^2 + CiM C4 = 1 + a3 + a2 — ai - 

C2 = a2 — ai + ttM C5 = 1 - Q!3 + a2 — ai - 

C3 = CiM C6 = -1 + a3 + Q!2 — ai 



2aM 
2aM 
' 2aM 



C7 = ai ~ aN 
cg = 1 + as — a2 



ai — 2a 



(7.45) 



(7.46) 



where Cfe > for all fc = 1 . . .8 due to positivity of the according factors in (|7.45p . The set of shorthands (|7.46p can 
be seen as a system of 8 equalities for the 5 variables ai, a2, as, am, a^- Thus we can eliminate these 5 variables by 
using 5 equations which results in three remaining consistency equations among the constants ci, . . . , cg. 



2ci = 2 + Cg + 2C7 — Cs 
2C2 = -2C3 + C5 + C6 

C4 — 2 + cq 



and Cfc > V fc = 1 , 



(7.47) 



The question is now how one can choose the maximal number of c's equal to consistently. By inspection of (|7.47p 
we can simultaneously set at most 6 of the c's to zero if: 



Case 1) ci — 1 
and 



C4 



C2 = C3 = C5 = Cg = C7 = Cg = 



Then the set (|7.46p obeys a solution: 

ai = 1 

a2 = 1 

as = 1 

aM = 

aAT = 1 



Case 2) cg = 2 C4 
and 



Ci = C2 = C3 = C5 



C6 



C7 = 



Then the set (|7.46p obeys a solution: 

ai = 

a2 = 
as = 1 
aM = 
aN = 



^^Certainly the ansatz H7.44|l qualifies only to estimate the largest inverse power of the maximal spin in the limit jmax — > oo- We could 
in principle assume the spins to be arbitrary polynomials of jmax- However we want to look for the smallest inverse power of jmax, that is 



we want the spins to be configured such that any higher inverse power of jn 
only be achieved if the spins scale with the same power of j'max ■ 



is compensated in Il7.43| l. In the limit jn 



00 this can 
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Case 1 
Here (|7.45l) reads as 



N) 



<- 



2(/32 -(3i + 2(3m - 1) 2(2j4 + /32 - A + 2(3n - ^) 



^ 2 (j4 + /52 + /?m) 2 {(32 Pm - 5) 2 - i)(2j4 + /^s + /32 - /3i + 2/3m)(-/33 + /32 - /3i + 2/3a/ - 1) 

1 



ip3 + P2-P1 + 2pM - 1) 2 (-/3i - /3jv + 2)(/33 -P2+(3i- + 2) 



1 



(7.48) 



However (|7.4ip gives in this case 
Y[ x{2l) < (32)4 



l=M 



(3 + 2/3i - 2/3jv) 



(3 + 2.74 + 2/3i-2/3m). 
(3 + 2j4 + 2/32 + 2Pm) 



(3 + /33 - /?2 + /3i - 2/3^) 



(3 + 2j4 + /33-/32+/3i-2/3M)_ 

(3 + j4 + /33 + J4 + /92 - /3i + 2/3m) 



(3 + 2j4 + 2/32 + 2j4 + 2/3jv) 
(1 + /32-/3i+/3m) 



(l + /?2-/3i+j4 + /?Jv) 



(3 + j4 + 133+ ji + P2-P1+ 2j4 + 2Pn) 
(2-/33 + /?2-/3i + 2/3m) 



(2 - /33 + /32 - /3i + 2^4 + 2I3n) 



{1 + Pm) 



{l+3i+PN) 



(2 + /33 - /34 + /?2 - /3i + 2/3Af ) 
(2 + /33 - /34 + /32 - /3i + 2j4 + 2/3w) 



(J4)^ 



We therefore find that in Case 1 Il^^^'^^I is bounded from above by a quantity of the order 



(7.49) 
Moreover 



according to ()7.37p the dimension of the matrix Q is given by dim Q = 2ji + 1 = 2j4 + 2/3i + 1 which is proportional 
to jmax- Therefore upon summing up aU terms in the row/column sums in order to apply theorem 17.2) we could 
get an additional factor of jmax- However, taking this sum implies increasing M and taking M — Pm ^ J4. This 
contradicts our minimization of the inverse power of jmax above immediately. Despite this we still may assume 
Ji ^ j2 < ji < and if we want the dimension D of Qjj to scale with the maximal spin then (|7.5p forces us to take 
ji ^ 32 ~ J3 ^ ji = Jmax- But then in (|7.46p all constants but C7, cg are non zero because for all a's we have < a < 1. 
Inserting this into (|7.48p and (|7.49p gives 



j^(M,N) 



1 



N-l 



and TT x{2l) ^ — unless /3 
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M = J4 



and therefore taking the row sum will not change the inverse power of jmax, because lim |L(*^'^)| 
thus the maximal row sum will be of an order (jmax)"'^ where r > 1. 



(7.50) 

(imax)^^ and 



Case 2 




Using 


7.44|) 



Ji = /3i 

j2 = P2 

ja = jmax + (33 

J4 — Jmax 



M = (3m 
N = (3n 



(7.51) 



Here (|7.45p reads as 
2 



- (M,N) 



< 



2(/32 - /3i + 2(3m - 1) 2(/32 - /3i + 2(3n - ^) 



2 {(32 + Pm) 2 {(32 - /3i + Pm - I) 2 (/3m - ^)(2j4 + P3 + P2-P1+ 2Pm){-P3 +P2-P1 + 2/3m - 1) 

1 

(/33 +P2-(3i+ 2/3m - 1) 2 (-/3i - /3jv + 2)(2j4 + ^3-^2+ Pi- 2/3^ + 2) 
1 



(i4)= 



(7.52) 



By observation, since M, N do not scale with the largest spin, (j7.4ip is only sensitive to M, N and thus becomes a 
constant in the limit j4 — + 00. Moreover x{2l) = 1 for M = N. Therefore in this case we can obtain a row sum of an 
order (jmax)^^- 
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Since ji < j2 we must have f3i < /32- Moreover ja < and thus /Sa = 0. We will also benefit from the fact that 
1 < A/ < iV < ji + i which implies, together with ()7.5ip . that 1 < (Jm < f^N < + ^- Additionally we observe that 
the row/column sum contains at most N = terms. Then we can continue: 



L 



(M,M) 



< 



< 



2(/32 - /3i + 2/3i + 1 - I) 2(/32 - /3i + 2/3i + 1 - i) 



2 (2/3m) 2 {Pm - 2 (/3m - ^)(2j4 + 2/3m)(2/3m - 1) 2 (2 - 2/3i - i)(2j4 - /Sa + A - 2/3i - 1 + 2) 



2(/32+/3i-^) 2(/32 + 2/3i + i) 



4 (2^4 + 1) 2 (2 - 2^1 - i) (2j4 - /32 + /3i - 1 + 2) 



< 



4/32 (/32 



■4(2j4 + l)2(i) (2j4-/32 + i) 



where in the second line we have set /3 



M 



1, /?, 



AT 



/3i + ^ in the denominator and /3 



M 



/3. 



N 



numerator. Moreover we have to approximate (3i — ^ in the third line due to the positivity of 



L{m,n) 



Pi 

2 



(7.53) 



i in the 
on the left 



hand side. Moreover since 



> by definition for all j4 



1 1 3 

2 ' ' 2 ' ■ 



we must have < /?2 < 2 in order to give 



a well defined positive expression. Setting /32 = 1 finally gives 



<- 



10 



10 



■4(2j4 + l) (2j4 + i) 



16(^4)2 + 12j4 + 2 



< 



< 



ij4 



(7.54) 



Having set /3i = ^ above implies that, according to (|7.37p . the dimension of the matrix is 2, and thus the row sum 
contains only a single term jL^^'^^j. We conclude that an upper bound for the eigenvalues of Q~jj^ is given by (|7.54p . 
However the finite upper bound (j7.54p for the spectrum of the inverse matrix provides us with a lower bound for 
the smallest non-zero eigenvalue oi Qd in leading order of the maximal spin jmax- 
This lower bound is then given for the spins ji ^ j2 ^ 1, ja = J4 = jmax as 



(7.55) 



This may be compared to the numerical result A^-I"*"-* = 2 • •\/jmax(jmax + 1) obtained in 24J contributed by the spin 

configuration ji = j2 — \ and ja = ji = jmax- Remarkably our estimate approximately reproduces this combination 

of spins. Recall that a|.7'"' = > M^. Therefore we conclude that: 
V V V 

Theorem 7.3 The smallest non-zero eigenvalue A^™-* ^ of the volume operator V acting on a gauge invariant 
four valent vertex v with spin configuration ji < j2 < Ja < ji — jmax is bounded from below by 

^!/™-vcrtcx(j--) ^ 4Vl^|-k(123)|-j™a. (7.56) 

where Z is the regularization constant in i2.25\} . (t(123) — 0, ±2, ±4, and £p is the Planck length. 



8 Summary 

In this paper and its companion [25] we have presented a comprehensive analysis of the spectral properties of the 
volume operator V in Loop Quantum Gravity (LQG), which may serve as a starting point for a computer based 
analysis of the action of constraint operators in the full theory. We have shown how the volume operator V can be 
handled analytically in case of a 4- valent vertex, and numerically in the general case. In this context we have analyzed 
how the spatial diffeomorphism invariant properties of the graphs underlying the states of the kinematical Hilbert 
space, as encoded in the sign factors e{IJK), become relevant in a practical computation. This aspect of the theory, 
although known, has not been analyzed in detail before the present work. 

A further summary and outlook is given in section 6 of the companion paper [25j . 
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Appendices 



This paper and its companion [25j make heavy use of special properties of the spin network functions, in particular 
the fact that the action of the flux operators on spin network basis states can be mapped to the problem of evaluating 
angular momentum operators on angular momentum eigenstates, the latter of which is familiar from ordinary quantum 
mechanics. This not only results in a great simplification, but also provides a convenient way to work in the SU (2)- 
gauge-invariant regime by using powerful techniques from the recoupling theory of angular momenta. 

We have added these appendices in order to provide the reader with details on the conventions we use, for example 
in the construction of recoupling schemes. Moreover the reader not familiar with these techniques is given a brief 
summary, in order to be able to redo the computations presented here. 

The appendices are organized as follows: In appendix [SI we give basic properties of matrix representations of 
SU{2), whose matrix elements are used for the definition of spin network functions. Appendix [Bl then summarizes the 
theory of angular momentum from quantum mechanics, and discusses an extended notion of recoupling of an arbitrary 
number of angular momenta in terms of recoupling schemes, as we use them. Appendix [C] completes this presentation 
with the definition of so called 6j-symbols, along with an elaboration of their basic properties. This provides an 
explicit notion of recoupling schemes in terms of polynomials of quantum numbers. Finally lAIBICI are connected in 
appendix [DJ in which the correspondence between spin network states and angular momentum states is reviewed in 
detail. 

A Representations of SU{2) 

Irreducible matrix representations of SU{2) can be constructed in (2j + l)-dimensional linear vector spaces, where j > 
is a half integer number; j = denotes the trivial representation. Note that we write for the defining representation 
of a group element h S SU{2) [tt^^ i (ft,)] ^^^^^ — Hab and for general representation matrices [7rj(/i)]^^^; the rows m are 
labelled from top to bottom hy m = j . . . — j and the columns n from left to right by n — j . . . — j, according to (jB.ip . 
In the following an overline denotes complex conjugation, while --^ indicates the transpose. 



A.l General Conventions — Defining Representation 
A. 1.1 Generators 

As generators of SU(2) we use the r-matrices given by :— — icfc, with ak being the Pauli-matrices: 
CTi = (^J (72 = (^J '^3= (^l \ ^ for which [CTi,(Tj] = 2ie*J'=CTfc 



(A.l) 



Additionally we use 



e = — T2 = [ ^2 ) with the obvious properties e ^ = e-^ = — e = f ^ ^ ] (A.2) 



A. 1.2 SU{2) Representations 

In the defining representation of SU(2), for a group element ft e G, we have 

h=(^_^^^ for which detft = lap + = 1 (A.3) 

where a,6 G (D. Moreover 

^ b a 

and we have the additional properties that 

ft= [ft-^]^= [eh^ e^f = ehe^ = ehe-^ (A.5) 

As mentioned above, we use the following convention for the matrix elements of ["'^=1 (^)] with A being the 
row and B the column index 

ftii=ftii=a hi 2 — hi 1 — b 

2 2 _ '-^ 2 2 _ (A 6) 

ft? 1 = h_i 1 = —b fto? = h_i _i = a \ ■ ) 

2 2 2 2 

For the r^'s we additionally have 

-T^^=Tfc = -(T,)-l (A.7) 
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A. 2 General Conventions for (2j + l)-dimensional S'[/(2)-Representation Matrices 
Here we follow [31 [H US]. 

A. 2.1 General Formula for SU{2) Matrix Element 

The {2j + l)-dimensional representatfon matrix [7rj(/i)] of /i e SU{2), given in terms of the parameters of the defining 
representation (|A.3p . can be written as 



^^^1 ' 1 ^^ + ~ + " ^77V-™-« ,m-rn'+l fT^l 



(A.8) 



where £ takes all integer values such that none of the factorials in the denominator gets a negative argument. By 
construction every representation of SU{2) consists of special unitary matrices [7rj(g)] (such that det [7rj(f;)] = 1 and 

[7rj((7)] [7rj(5)]^ = [7rj(g)]^[7rj(5)] = [7rj(l)] (with [7rj((7)]^ = ["'5(5)] ))• Using (|A.8|) we obtain a generahzation of 
(TOI) : 



A. 2. 2 Generators and e-Metric 

Upon applying the representation matrix element formula (|A.8|1 and the ansatz, 



(A.9) 



(A.IO) 



where in the defining two dimensional representation the exponential can be explicitly evaluated, one obtains|fl using 
A{j, m) = + 1) — m(m — 1)] ^ with the obvious property that A{i^ — m) = A{j^ m+1), 

['^j('^i)]„i„ ^ ~i A{j,m) 5ran+i ~ i A{j,m+ 1) (5„„i_i 
['^j('^2)]„j„ = - A{i,m) 5m n+i + A(j,m + 1) (5™„_i 

K(t3)]„„ = - 2 i m (5,„„ (A. 11) 

Moreover one finds by plugging (|A.2|) into (IA.8P that 

bMmn - (-l)'""" and h(6-^)]_ = (-1)^+" (A.12) 



B Angular Momentum Theory 
B.l Basic Definitions 

In this section we will summarize the results given in |29| . We have an angular momentum orthonormal basis 
u{j,m;n) — | j to ; n ) of a general {2j + 1) dimensional representation of SU{2) where a general state can 
be represented as a column vector with components 



component in 
component in 

component in 
component in 



j m—j ;n ^ direction 
j TO=j — 1 ;n ) direction 

j m = —j + l ;n ^ direction 
i m — ~j ; n ) direction 



(B.l) 



Note that we will adapt here the general notation where the components are labelled from top to bottom hy j, j — 1, 
• • • , —j + 1, —j. The index n stands for additional quantum numbers, not affected by the action of the angular 
momentum operators fulfilling the commutation relations 

[J\J^] = i e'^^ J'' (B.2) 

We can formulate ladder operators given as J+ = + iJ^ and J~ = — iJ^ where, using the shorthand 
A[j,m) = \J jij + 1) ~ to(to — 1) with the obvious property that A(j,—ra) — A{j,m+1), we have 

J~ I J to ;n ^ = j4(J;TO)| j to— 1 ;n ^ 

J^l j m — A{j,m + 1)\ j m+l ;n ) (B.3) 



33 See [22l[23] for details. 
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The I j m ; n ) simultaneously diagonalize the two operators: the squared total angular momentuno {JY and the 
magnetic quantum number |29| : 

{Jf\ J m ;n ) = {{J^fHJ'fHJ'f)\ ]rn,n)^ (j- /+ + J+ J") + ( j3)2) | ^ ) 

- + \jm-n) (B.4) 
J'^ I J TO ; n ) = TO I J TO ; n ) 

That is, I j TO ;n ) is a maximal set of simultaneous eigenvectors of (J)^ and . 
One then finds the following commutation relations 

[j3,(J)2] = [J+,(J)2] = [J-,(J)2] = o [J^J+] = ,/+ [j\j-] = -J- [J+,J-]=2J3 (B.5) 

such that we obtam for the (2j + l)-dimcnsional matrix representation with arbitrary weight j 

{J )m'm — ^ ^m' m 

{J^)m',n = -4(j, TO + 1) Sm' m+1 {J^)m'm = + J")m'm = ^^(j, TO + 1) (5„' ,„+i + ^A{j, to) 5„/ „i_i (B.6) 

{J'')m'm = ^(j, "l) (J^)m'm = ^ -'"")m'm = ^(j, W + 1) (5m' + f to) 5m' m_l 

In the defining representation {j = ^) we have 
to' 

i TO — > 

ji=i(j++j-)=(''; t") j2=-i(j+-j-) = ('o 5 



\ 2 0/ vi 

where we see the general construction principle of J'^ = [f^ — 1), the Uk being the Pauli matrices as indicated in 
(|A.ip . The indices (m', m) labelling the row and column respectively of the matrix elements run as indicated, from 
to' = TO = j = ^ on the upper left corner until m' = m = —j = — ^ in the lower right corner. 

B.2 Fundamental Recoupling 

Equipped with a small part of representation theory we can easily understand what happens if we couple several 
angular momenta. For that we first repeat the well known theorem of Clebsch & Gordan on tensorized representations 
of SU{2): 

Theorem B.l Clebsch & Gordan 

Having two irreducible representations TTj-^, ttj^ of SU{2) with weights ji and j2, their tensor product space splits into 
a direct sum of irreducible representations tt^j^ with \ji — j2\ < ji2 < ji + j2 such that 



Equivalently we can write for the resulting representation space = 'H(Di) ® "^(Da) (where Di ~ 2ji + 1, £'2 = 

2j2 + I, D = Di ■ D2 denote the dimensions of the Hilbert spaces): 

ji +32 

= ^(Di) ^(D2) = ^ ^(2ji2 + l) (B-8) 

ii2=|ii-j2| 

In other words, if we couple two angular momenta ji,j2, we can get resulting angular momenta ji2 varying in the 
range \ ji — ^2! < ji2 < ji + 12- The tensor product space of two representations of SU{2) decomposes into a direct 
sum of representation spaces, with one space for every possible value of recoupling ji2 with the according dimension 
2jl2 + 1. 

^-^Note that we denote single components of the angular momentum by J', i = 1,2,3, whereas the total angular momentum is denoted 

by (J). 
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B.3 Recoupling of Two Angular Momenta 

According to (jB.Sp we can expand, for each value of ji2, the elements | ji j2; ii2(ji,j2) M ) € "^(2^12+1) (called 
'coupled states') into the tensor basis | ji wi ) (g) | j2 m2 ) of H^'^^: 

I ii 32] ji2{ji,j2),M ) = ^ ( ji mi; j2 1712 I ji 32; jl2{jl,j2),M ) I TOi; j2 ) (B.9) 

mi+m2=Af " ' 

Here the expansion coefficients Cmim^ S IR- are the so called Clebsch-Gordan coefficients. On the right hand side we 
have introduced | ji mi; j2 m2 ) :— | ji mi ) (g) | j2 "I2 ) 

If we change the order of coupling then Cmim2 changes its sign: 

{ ji mi] j2 1712 I ii2(ii,j2) M = mi+m2 ) = (-l)Ji2-Ji-J2 / I j^^lj^^ji) M^mi + m2 ) 

(B.IO) 

= j2 TO2; ji mi I ji2ij2,ji) M^mi + m2 ) 

Since =Fji2 i (ji + J2) is an integer number we are allowed to flip the overall sign in the exponent. The coupled states 
again form an orthonormal basis: 

( jl j2; jl2{juj2),M I ji j2; jl2{jlj2),M ) = 5~ S^j^ (B.ll) 

In (jB.lip comes from the orthogonality of the W(2ji2+i) i'^i (jB.Sp . Sj^j^^ is caused by the orthogonality of the 

single I jl mi ), | j2 1712 ), since always AI = mi +m2. Normalization of the recoupled states (|B.11|1 implies, according 
to the orthogonality of the angular momentum states: 

|( -^'l -^^i 312{J1.J2),M I Jl J2; 312[3UJ2).M )|' = \Crn,ra,? = 1 (B.12) 

mi+m2— A/ mi+^^2 — A/ 

Furthermore the Clebsch-Gordan coefficients are all real, which is not obvious, but a result of two conventions one 
usually requires p5] : 

1) I jl 32; ji2{ji,j2) = ji+j2 M = ji+j2 ) = I Jl mi; j2 m2 ) 

2) All matrix elements of J^jj-^-^ which are non-diagonal in | ji j2] ji2{ji,j2) M ) are real and nonnegative. 

The maximal set of simultaneously diagonalizable (that is commuting) 2x2 operators (IB.4|) of the single Hilbert 
spaces 7i(£ij), 7^(02) is then in 'H(d) replaced by four operators: total angular momentum (J(_d))^, total projection 
quantum number Jf^j^y and single total angular momenta {J(Di))'^^ {J(d-,))'^'- E3 

{J(D)f ^ {J(Di) + J{D2)f J(D) ^ J(Di) + J(D2) iJ{Di)f [J (D^)?' (B.13) 

which are simultaneously diagonal in the new basis: 

{J(D)f I jl j2] jl2{jl,j2),M ) = jl2(jl2 + 1) I jl j2] jl2{jlj2),M ) 

Jfo) I 31 32; jl2{jl,j2),M ) = M I jl j2; jl2{jl-j2),M ) 

{J(Di)f I jl j2; jl2{jl,j2),M ) = ji(ji + 1) I jl j2; jl2{jl,j2),M ) 

{J{D2)f I 31 32; il2(jl, j2),M ) = J2(j2 + 1) I 31 32; jl2(jl,j2),M ) (B.14) 



^^Here we slightly abuse the notation. More correctly we should write: 
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B.4 Recoupling of Three Angular Momenta — 6j-Symbols 

In this way we can expand the recouphng of three angular momenta in terms of the Clebsch-Gordan coefficients. 
I ji2(ji,i2), i(ji2,j3) ) = I ji2(ji,j2) ja; jUuJs) M = mi+ni2 + ni3 ) 

= X! ( -^12 mi2; js 1713 I ji2 js; j mi2 + ) ■ 

I jl2 TO12 ) (Xl I J3 'TIS ) 

= ^ ( ji2 "I12; J3 TO3 I ji2 js; 3 mi2 + 7713 ) • 

( jl "T-i; 32 m2 I jl j2; jl2 "712 = 7"! + 7712 ) • 

mi 7712 

I jl 7771 ) (g) I j2 "72 ) (Xl I i3 "73 ) 
= X! ( -^12 "71 + "72; js M -1711-1712 \ jl2 js', j M ) ■ 

mi 7712 

( jl 7"i; J2 "72 I jl j2\ jl2 "71 + "72 ) • 

I jl 7771 ) I J2 "72 ) «) I js M -mi -1712 ) (B.15) 

As we can see, as we couple angular momenta successively, the order of coupling plays an important role. Different 
orders of coupling will lead to different phases of the wavefunctions (see (jB.lOl) ). Because of this it is useful to have 
a transformation which connects different orders of recoupling. The transformation between the two different ways of 
coupling 3 angular momenta ji, j2, js to a resulting j are the 6j-symbols; see appendix [Cl 

B.5 Recoupling of n Angular Momenta — 3nj-Symbols 

As mentioned earlier the successive coupling of three angular momenta to a resulting j can be generalized. For this 
purpose let us first comment on the generalization principle before we go into detailed definitions. 

Theorem IB. II can be applied to an n-fold tensor product of representations ttj^ (g) iTj^ (g) . . . (g) ttj^ by reducing out 
step by step every pair of representations. This procedure has to be carried out until all tensor products are reduced 
out. One then ends up with a direct sum of representations, each of which has a weight corresponding to an allowed 
value of the total angular momentum to which the n single angular momenta ji, j2, ■ • ■ , jn can couple. However, there 
is an arbitrariness in how one couples the 77 angular momenta together, that is the order in which tTj^ (g) TTj^ (g) . . . (g) tt^^ 
is reduced out (by applying IB. ip matters. 

Consider a system of n angular momenta. First we fix a labelling of these momenta, such that we have ji, j2, ■ ■ ■ , jn- 
Again the first choice would be a tensor basis | j rn ) of all single angular momentum states | jk TUk ), k — 1 . . .71 
defined by: 

n 

I J "7 ) = I (jl, j2, • ■ ■ , jn) {mi, 7172, ■ ■ • , Wn) ) | jk 7"fc ) (B.16) 

fc=l 

with the maximal set of 2r7 commuting operators {JiY , Jf, (/=!,..., n). 

Now we proceed as in section rB.3l finding the commuting operators according to (|B.14|) . that is a basis in which 
the total angular momentum {Jtot)^ = (J)"^ = (Ji + J2 + • • • + Jn)'^ is diagonal (quantum number j) together with the 
total magnetic quantum number J^^^ — J'^ ~ + J2 + ■ ■ ■ + (quantum number M). As (J)^ and are only two 
operators, we need 2{n — 1) more mutually commuting operators to have again a maximal set. We choose therefore 
the 71 operators {JiY, I = 1, ... ,71 of total single angular momentum (quantum numbers (ji, . . . , j„) := j). Then we 
are left with the task of finding an additional n — 2 commuting operators. For this purpose we define: 

Definition B.l Recoupling Scheme 

A recoupling scheme | g{IJ) j j tti ^ is an orthonormal basis which diagonalizes, besides (J)^, , and {JiY (I = 
1, . . . ,n) the squares of the additional n — 2 operators G2, G3, . . . , G„_i defined as: 

Gi := Ji, G2 ■= Gi + Jj, G3 := G2 + Ji, G4 := G3 + J2, ■ ■ ■ , 

Gi G/_i + J/-2, G/+1 := G/ + J/-i, G/+2 '■— G/+i + J/+i, G/+3 := G/+2 + Ji+2-, ■ ■ ■ , 
Gj G,7_i + J,/-i, Gj+i := G,7 + Jj+i, Gj+2 ■— Gj+i + Jj+2, ■ ■ ■ , 

Gn-l ■— Gn-2 + Jn-1 
36Note that formally G„ := Gn-i + Jn = Jtotal- 
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The vector 

diU) ■= { 92iji,jj), 53(52, ji), 9i+i{9iji-i), 9i+2igi+i,ji+i), gj{9j-i,jj-i), 

5n-l(5n-2, jn-l)) 

carries as quantum numbers the n — 2 eigenvalues of the operators (02)^, • . • , (Gn-i)'^ ■ 

So we recouplc first the angular momenta labelled by /, J where I < J and secondly all the other angular momenta 
successively (all labels are with respect to the fixed label set), by taking into account the allowed values for each 
recoupling according to theorem IB. II 

Let us define furthermore the so called standard recoupling scheme or standard basis: 

Definition B.2 Standard Basis 

A recoupling scheme based on the pair (/, J) = (1,2) with 



K 



L=l 

is called the standard basis. 



Using definition IB. II with the commutation relations (jB.Sp and the fact that single angular momentum operators 
acting on different single angular momentum Hilbert spaces commute!^, one can easily check that for every recoupling 
scheme 



(i) the G/'s fulfill the angular momentum algebra TB. 51 

(ii) ( J)^, (JiY: {GiY, J'^ V/ = 1 . . . n commute with each other 

Note that it is sufficient to prove these two points in the standard basis 5(12), because every other basis g{IJ) is 
related to it by simply relabelling the n angular momenta. 

We have thus succeeded in giving an alternative description of a system of n angular momenta by all possible 
occurring intermediate recoupling stages G/, instead of using the individual magnetic quantum numbers. Every 
recoupling scheme | g{IJ) j j m labeled by the index pair (/J), gives a distinct orthonormal basis. We are then in 
need of a transformation connecting the different bases, i.e. we wish to express one basis, belonging to the pair (/J), 
in terms of another basis belonging to the pair (KL). This leads to the following definition. 

Definition B.3 3nj-Symbol 

The generalized expansion coefficients of a recoupling scheme in terms of the standard recoupling scheme are called 
3nj -symbols. ■ 

I g{IJ) i j m) ^ ^ {g' (12) j 3 m \ g{IJ) j j m) \ (f (12) j j m ) 

all 3' (12)"" " ' 

3nj — symbol 

The summation has to be extended over all possible values of the intermediate recouplings 

g'(12) = (52(^17^2), 53(52;i3)i • ■ • I 5n-i(5n-2i Jn-i); ^^^^ valucs of cach Component g'f. allowed by theorem \B.l\ 

In calculations we will suppress the quantum numbers j,j,m, since they are identical in every expression, and 
write the 3nj-symbols as ( g{IJ) \ 5'(12) ). Note additionally the following properties of the 3nj-symbols: 

(i) They are unitary and real, due to the the fact that they can be expressed in terms of Clebsch-Gordan- 
coefficientfl^l: 

( g{IJ) I 5'(12) ) - ( .g'(12) I g{IJ) ) 

(ii) They arc rotationally invariant, i.e. independent of the magnetic quantum numbers nik occurring in ()B.16|) . 

B.6 Properties of Recoupling Schemes 

In this section we will briefly review the properties of recoupling schemes as defined in section IB.5I In what follows 
we will frequently use 



TTlfc := mi + TO2 + • ■ ■ + ITlk 



S^That is [Jj, J'j] = whenever I ^ J. 
^^Which are unitary and real. 
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B.6.1 A General (Standard-)Recoupling Scheme 

A general standard recoupling scheme is defined as follows: Fix a labelling ji, . . . , of the N spins to recouple. Then 
one constructs 

I a(12) J M ; n) = 

= I a,2{j\ 32) 03(0-2 33) ■ ■ ■ aK-i{aK-2 Jk-i) a,K{aK-i Jk) a,K+iia.K Jk+i) ■ ■ ■ a]V-i(ajV-2 Jn-i) Jio-N-i Jn) M ; ni ...tim) 
= ^ ( ii m-i ; 32 m2 I a2(ii J2) m2 ) 

mi + ...+mjv=M { a2 lfl2 ; js "13 | 13(12^3) "13 ) 

( aK-2 fhK-2 ; m/f-l I a,K-i(aK-2 3K-l) fhK-i ) 



a.K-1 rriK-i ; 3K rriK 
a-K rriK ; 3K+1 rriK+i 



aN-2 'm.N~2 ; iiV-i niJV-i | ajv-i (a]v-2 jjv-i) miV-i ) 
o-iV-i m-iV-i ; iiV miV | JiaN-i3N) M ) 



N-1 factors 



I ji mi ; 711 ) ® I j2 "12 ; "2 ) 



jjv_i m]v-i ; Jijv-i ) 



Jif TTiK ; ) ® I iif+i m^+i ; n/c+i ) ' 



(B.17) 



B.6.2 Orthogonality Relations Between Recoupling Schemes 

For the scalar product of two recoupling schemes we have: 

(a2«3 ••■ a'N-l J' \'^2 ■■■ a^-l J M) ShlM'Sa^a'^Sa^a'.^ ■ ■ ■ 5af,_^a'^_^5jJ> 5^^^ . . . Sy^j^S^'^m ■ ■ ■ Sn'j^UN 

This result can be easily understood by recalling the definition of a recoupling scheme 

I ci2(ji j^) 03(02 jz) ■ ■ ■ aN-i{o.N-2 Jn-i) J{o-n-i Jn) M ) as the simultaneous eigenstate for the operators (62)^ = 

(Jl + J2)2, (G3)2 = (G2 + Ja)^, . . . , [Gn-I? = (^^-2 + Jn-1?, = {Gn-I + Jiv)' = (Ji + . . . + J^)2 with 

eigenvalues 02(02 + 1), 03(03 + 1), . . . , gN-i{gN-i + 1), ^(^ + !)• 

B.6.3 Partial Orthogonality Relations Between Recoupling Schemes 

The same argument can also be applied to cases where we have to calculate the scalar product of two recoupling 
schemes of different recoupling order. For illustration let us consider two recoupling schemes 

\d J M) = \a2{jij2) a.3{a2 js) . . . aj(-i{aK-2 3K-1 aK{aK-i 3 k) ■ ■ ■ o-LiaL-i Jl) Ql+i ("l il^i • ■ ■ CiV-i (a]V-2 iiV-i) ^(ajv-ijjv) M) 

I g J Af ) = I 32 (ii j2) 93(a2 is) • ■ ■ gK-i{9K-2 3k-i giiigii-i 3p) ■ ■ ■ 5l(9z^i 3r) giA-iidL jL+i) ■ ■ • 9]V-i(3jv-2 3n-i) J{gN-i 3n) M) 

Here from 2 ... if — 1 the spins ji, j2 • • ■ jK-i are coupled in a and g in the same order. Then Jk ■ ■ ■ jh are coupled in 
the standard way for a but in a different order for g. After that jl+i ■ ■ -Jn are successively coupled to each scheme 
again. Now it is clear that o and g simultaneously diagonalize not only {G2f' — {Ji + 1 (G'a)^ = (G2 + Js)^, . . . , 
{Gk-iT ~ {Gk-2 + Jk-i)^ but also [GlY . . .{Gn-iY , — ( Ji + . . . + JnY- Therefore we can write down 
immediately 

{a J M \g J M) = (axiaK-ijK) ■ • ■ ol_i(ol_2 jz^i) aL{aL-ijL) \ gxiaK-ijp) ■ ■ ■ gi^iigL-2 jq) cLLigi^ijn) ) 

'xSa2g2 ■ ■ - SaK-igK-l ^ '^ai.ffi • ■ • (^ajv-igjv-l X ^JJ' X ^MM' (B.18) 

For a more detailed derivation of (jB.lSp see lemmas 5.1 and 5.2 of jT5] . 

C Properties of the 6j-Symbols 

In this appendix we will give an overview of the 6j-symbols, because they are the basic structure we will use in our 
recoupling calculations, as every coupling of n angular momenta can be expressed in terms of them. For further details 
we refer the reader to [5U].[55]. 
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C.l Definition 

The 6j-symbol is defined in p 92, as: 

:= [(2ji2 + l)(2j23 + 

x( jl2(jl, j2), j(jl2,j3) I j23{j2,.h),3{jl,323) ) 
= [(2ji2 + 1)(2J23 + 

X ^ { jl mi; j2 7712 I jl J2 ji2 mi + m2 ) 

mi m2 

X ( ji2 mi + m2; ja m - mi - m2 \ ji2 33 3 m ) 
X ( j2 m2\ 33 m-mi-m2 \ 32 33 J23 m - mi ) 

x( jl mi] j23 m - mi I ji ^23 3 ^ ) (C.l) 
The factors in the summation are Clebsch-Gordon coefficients. 

C.2 Explicit Evaluation of the 6j-Symbols 

A general formula for the numerical value of the 6j-symbols has been derived by Racah [29], p. 99: 

jl j2 il2 1 ^ , . . . , ^ , . . . . ^ , . . . ^ ^ , . . . ^ I jl i2 jl2 

= ^{31, 32, 312)^(31, 3, 323)^(33, 32, 323)^[33, J, , Jl2)W < 
33 3 323 I [ 33 3 323 

A(a,5,c) : 



'(a + 6 - c)!(a - 5 + c)!(-a + 6 + c)! 



X 



(a + fc + c+ 1)! 

jl j2 jl2 I .r-^ 

u;< V 2^(-l)"(n + 1)! X 

j3 j j23 J « 

X [{n - jl - 32 - jl2)!(»^ - jl - j - j23)!("- - j3 - j2 - j23)!('^ " j3 ^ j - jl2)!]"^ 
X [(jl + j2 + j3 +j~n)\ (j2 + jl2 + j + j23 - n) ! (ji2 + jl + j23 + j3 - !] (C.2) 

The sum is extended over all positive integer values of n such that no factorial in the denominator has a negative 
argument. That is: 

max[ji + j2 + jl2, jl + j + j23, j3 + j2 + j23, j3 + j + jl2] < < min[ji + j2 + j3 + j, j2 + jl2 + j + j23, jl2 + jl + j23 + j3] 



Remark From (|C.2p we are provided with some additional requirements the arguments of the 6j-symbols must fulfill: 
Certain sums or differences of them must be integers to be proper (=integer) arguments for the factorials: 

From A(a, h, c) one gets: 

• a, b, c must fulfill the triangle inequalities: (a + & — c) > 0, (a — + c) > 0, (—a + b + c) > 0, 

• (±a ±b ± c) must be an integer 

From the w-coefficient one gets: 

• jl + j2 + j3 + j, j2 + ji2 + j + j23, ji2 + jl + j23 + j3 are integers. 

The following (trivial but important) relations are frequently used in calculations involving 6j-symbols: 

(-1)" = (-1)-" VzeZ 
{-if'- = 1 Vz e Z 

(-1)^'= = (-l)-*^ ^k^^ with z el (C.3) 
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C.3 Symmetry Properties 

The 6j-symbols are invariant under 
• any permutation of the cohimns: 




• simultaneous interchange of the upper and lower arguments of two columns, e.g. 




C.4 Orthogonality and Sum Rules 
Orthogonality Relations 



323 [ J3 J J23 

Composition Relation 




in 


h 


jl2 1 


II 




33 


323 1 




f 33 31 


331 \ 




3 


J23 J 






3 


331 J 




[h 3 


312 J 



Sum Rule of Elliot and Biedenharn 



31 32 312 I I 323 31 Jl23 
33 jl23 ^23 J [ 34 3 314 

(^_^yi+j2 +j3 +ji +il2 + j23 +il4+ jl23 +j 



33 32 323 [ j 32 31 312 [ j 33 312 3123 
314 j 3124 I I 34 Jl24 314 I I 34 j jl24 



(C.4) 



(C.5) 



^(_l).23+,3.+,.2(2j,3 + 1) .... ..3 . J3 ^^^^^ 



xE(-i)-- i2n,, + i) { - - - - - \ { - ) (C.8) 



D Spin Networks and Representation Theory 

Concrete calculations performed in Loop Quantum Gravity heavily rest on the fact that the representation theory 
of SU (2) is closely related to the theory of angular momenta familiar from quantum mechanics. Employing this 
correspondence we are able to use the powerful techniques provided by recoupling theory of angular momenta in Loop 
Quantum Gravity. However, as we will describe in this section, there are certain subtleties which have to be worked 
out carefully in order to fully establish the indicated correspondence. 

D.l State Vectors and their Transformation Properties 

Let be a linear vector space equipped with a basis e^, where a uniquely labels the basis statej^. With respect to 
the given inner product (•, •) on 7i the basis vectors e^ should be normalized, orthogonal 

{&a'^^a) = ^a' a (D.l) 

and complete 

lH = E(-,eo )(e„ •) (D.2) 



^®The index a can well be a multi-label, as in the case of a quantum system described by more than one quantum number. 
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A general vector 4' e 7i can then be expanded into components with respect to this basis as 

M/^J^cWe^ (D.3) 

a 

with components Cq*'' = ( Gq,^). Of course the choice of basis vectors is not unique: 5* might equivalently be 

expressed with respect to a different set of basis vectors e^. Let O be a hnear invertible map acting on the components 

(*) J. 
Ca as on a row vector: 



Then we have for fixed ^': 



a a' [3 



a a'f3 

=--T.^Pe, (D.5) 

13 

Identity (jP.Sp nicely illustrates the transformation property of the basis states ep = ((C'""'^)"^)^^, with the 

inverse transpose of O. This is referred to as a contragredient transformation behavior [32 with respect to O. Moreover 
if O is unitary then (|D.ip , (|D.2p remain valid for as well. 



Active and Passive Transformations 



Let the state vector vjf itself now undergo a transformatiorri: Using again convention (jD.4p a transformation A acts 
on 5* as: 

ivf^vl/^^^cW i„,„e„, (D.6) 



Now there are two interpretations of the transformation (|D.6p : 

Active Interpretation The basis e^/ is not changed, but the components Cq*-* of 5* change according to 

c(*) := ^ ci*) and $ = ^ c^V e^, (D.7) 

a a' 

Passive Interpretation The basis ea> is changed, the components of A are the expansion coefBcients of the new 
basis vectors e^ in the old basis vectors ea'. The components cL*^ of ^' stay unchanged as 



D.2 Spin Network Functions 

D.2.1 Defining a Basis of Angular Momentum Eigenstates 

Due to the Peter & Weyl theorem we have [T] : 



/ 

JSU{2) 



1 



d^lH{h) [^j'ih)]^^,^, [7rj(/i)]^„ = -;^—^5fj5„,,,n5n'n (D.9) 



which states that the representation matrix element functions [7rj(/i)] are orthogonal with respect to the Haar 
measure dfinih) on SU{2). Theorem (|D.9p enables us to introduce an orthonormal basij^ on the total Hilbert space 

7i = 0W(^^ (D.IO) 
j 

by 

{h\jmn) -.^ ^/2jTl[7rj{h)]^^ and {jmn\h) :- + 1 h W]„„ - ^2^' + 1^ (D.ll) 



'^''e.g. a rotation 

^'^We will use Dirac's bracket notation here. 



where by construction 



{jmn\fm'n') / dfiH{h){jmn\h){h\fm'n') = ^(2/ + l)(2j + 1) / dtiH{h)[n,{h)]^Jny{h)]^,^, 

JSU{2) JSU{2) 

= 5f.j5rn'm5nn (D.12) 

Therefore the meaning of the representation matrices is twofold: On the one hand the rescaled representation matrix 
element functions \/2j + 1 [7rj(-)] =: | j to n ) provide an orthonormal basis, whereas the representation matrix 
[7rj(/i)]^^^^ of a specific group element h G SU{2) should act as a unitary transformation matrix on the basis states 
according to (ID.8|) as 

m' 

D.3 Right Invariant Vector Fields 
D.3.1 Action on Group Valued Functions 

The action of right invariant vector fields (^Xkf){g) ~ lt-o group valued functions f{g), g £ SU{2) is 

evaluated as: 

(X,/)(5) = |/(e*^'=5) 



t=o 



t=o 



ABC ^ H)^o 

where we have used the defining [j — i) representation of SU (2) in order to define the derivation of the group element. 
In the last line we then express the sum over the defining representation indices as a trace. Note, however, that using 
the defining representation here is a matter of choice, we are free to choose any other representation and accordingly 
the trace in (jD.14|) is then defined with respect to the chosen representation. From (jD.14|) the commutation relations 
of the Xk are found to be 

([X,,X,]f){g) = -2^^\Xkf){g) (D.I5) 
The action of the right invariant vector fields on representation matrix element functions can be evaluated as: 



t=0 



dt 

r 

= 1] hK-)]™,h(5)]™ 

r 
r 

where we use the representation propertjQ a-nd the definition of the Tk 's as a basis for the tangent space of the group 
at the identity element lsu(2)- Note that we have changed the notation in the last line to indicate that the resulting 

(2j + 1) dimensional matrix [Tfc]["'L is an element of the Lie algebra su(2) rather than a group element. Moreover as 
mentioned above we choose here the representation of weight j in order to evaluate the trace. 

However (ID. 161 ) yields a possible obstacle: while we would like to write the action of any transformation on the 
basis states (|D.11[) to be according to (|D.13p . the action of the right invariant vector fields (|D.16|) is defined as a 
matrix multiplication, that is it transforms the basis states (jP.lip as vectors rather than as basis states. Therefore 
we have to take a closer look at the group multiplication properties. 



[^3(9192)] = [i-jCsi)] [-^3(92)] 
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D.3.2 Group Multiplication — Introducing a New Basis 

For the representation of a general group multiplication we find 

= H [^j(^'^)]m"r ^ ^ ) 

r s m" 

r s m" 

By observation, (|D.f 7|) suggests a solution to the problem above: If we slightly modify our preliminary basis (|D.11|) 
by absorbing the according e components, we can introduce a new basis | j to ; n ) as 



J TO ; ri 



= (-i)^-™v/2jTTK-(.)]_„ 

( i TO ;n I := ^ J> " I = E J ^ ^ I = (-ir™( J 



= {-iy-^^2jTl [^.(•)]_„ (D.18) 



which is orthonormal as well by the unitarity of [vTj (e)] Pn . Note that there is an ambiguity in the definition (|D.18|) 
coming from the fact that due to (jA.9|) [7r(e^^ g^^ e)] = [7r(e g^^ c^^)] > ^"^^ therefore we can exchange ^ e 
consistently in (jP.lSp . 

Now continuing from (|D.17p . and using (jP.lSp . group multiplication can be rewritten as 

= /2747t E K-(^'')]m. h(.9"')]rs (^1 i ^ > 

We have now formally achieved a form analogous to the transformation property (ID.13|) on the redefined basis (ID.18|) . 
D.3.3 Rewriting the Action of the Right Invariant Vector Fields 

Let us evaluate the action of the modified right invariant vector fields Xk on the representation matrix clement functions 
[7rj(-)]^^. By using (|A.7p as — t|" — ^ "'''k^ can evaluate the action of the Xk on the states | j m ; n ) as 
follows 

Xk {h\ J TO ;n ) = V27TT(-l)-'-™(x4^,(.)]^ J W 
= V27TT(-ir™|([.,(cB*-/^)]_ J 

^E(-i)'^^'"™^l(h(-*"i)];t 



{h \ j r ■,n) 



d 



{h \ j r - n) 



= E(-i)''''"™^rf^([-^('^"*"i)]. 

r 
r 

where the action of the right invariant vector fields may now be interpreted as a basis transformation. The prefactor 
(_l)2(i-m) = I because J — TO is an integer number. Here we again have changed the notation in the last line to 
indicate that the resulting (2j + l)-dimensional matrix [t^J^^'L is an element of the Lie algebra su{2) rather than a 
group element. However, we will abuse the notation a tad by writing [T/c]r''m as [7rj(Tfc)] 



Also this can easily be seen from llm2l l and the fact that (-1)2^-™) = 1. 
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D.4 Correspondence to Angular Momentum Theory 

Now wc introduce the modified vector fieldfQ Y'^ := —^Xk, and their commutation relations as evaluated on group 
valued functions, which easily follows from (jD.15|l : 

[Y\Y^ {h\jm-n)= ie^^^Y^ {h\jm;n) (D.21) 

These are the commutation relations of angular momentum operators. Moreover using (jD.20p the Yk act on the basis 
states (ID.18|) as 



(D.22) 



Thus we can identify the action of Y^ on the new basis states | j m ; n ) with the action of usual the angular momentum 
operators Y^ = on basis states of an abstract spin system |j to)^^^^F^: 



Y^ \j m-n) ^ j'' \j m) 



(") 



(D.23) 



and obtain the usual angular momentum and ladder operator algebra and the algebra of ladder operators as can be 
found in [5^. This of course implies 



Y^ I j to' ; n ) 



J 



Y I J m' ; n ) 



= m' + l) |j to' + 1) 

t 

= [^J-(^2)]^„,J| jm;n)=A(j,TO' + l)| jTO' + l;n) 



A{j, m') I j to' - 1) 



(«) 



t 

1 J2 [h(^i)]mm' - i [^J" (^2 I J TO ;n ) = A{j,m')\ j m'-l ;n ) 

m=-j 



J3 



j in'). 



(n) 

j m! ;n 



t 



i' \j to') 



(n) 



r3| j m' - n) I J2 [^j(^3)]„„,| j m;n) 

m=-j 

j{j + 1) \j to') 



TO 



'I j to' ;n ) 



j to') 



(») 



'(") 



t 

j m' ;n) = j{j + 1)\ j m' ;n) 

where (y)^ ^ (Y^)^ + {Y^ + {Y^f = ^Y-Y+ + ^Y+Y- + (Y'^f. Upon defining the action of the F's in this 
way there is a correspondence Y^ <-> and therefore one realizes the following commutation relations (which are 
compared to the commutation relations of the usual angular momentum operators acting on spin states | j to )). 



Y\Y^] ^ie'-i^Y^ 



(J)2,J'=]=0 



Y^,Y- 



Y- 



Y+,Y-] = 2Y^ 



j+ 

-j- 

-"1 - 2 J3 



(D.24) 



D.4.1 Rewriting Right Invariant Vector Field Expressions 



We can now give the general prescription for translating the action of the modified right invariant vector fields Y^ on 
holonomies into the action of usual angular momentum operators on abstract spin systems. 



Y'^ V'2jTl[^,(-)]_ = Y'' I -TO ;n ) 



4^ 



(D.25) 



**Note that we can raise and lower sii(2)-indices with the Cartan metric on sn(2). Since this is simply &ij = 5''^ , there is no difference 
between upper and lower indices. We thus place sii{2)-indices as is convenient throughout the paper. 
"^^Here n denotes an additional quantum number which is irrelevant for the action of j''. 
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D.4.2 Recoupling Schemes 

Due to (|D.25[) we can therefore apply recoupling theory to the states | j m ; n ) as follows: 



a J M ; n j ) := ^ { jimi ; j2 m2 | a2 m2 ) x . 

m 1 + "Ti 2 + ■ ■ ■ + "Ti iv — 

X I ji mi ; ni ) 



D.5 Historical Remark 



X ( aN-irriN-i ; jNrriN 



JM ) 



JN ruN ;nN ) 

(D.26) 



For the sake of completeness we would like to add a remark on the development of correspondence (jD.25[) here. The 
necessity of (|D.13I) leading to 



Y\h\jm;n)^\Y.[^,{r,)l,^^{h 



] r \n 



) 



(D.27) 



was realized for the first time in [551 US] j where the authors evaluated the detailed matrix elements for the Yk ■ In 
calculations in the spin network formalism before [551 US] one always used the preliminary basis (jP.lip and directly 
evaluated the action of the Yk resulting from the group multiplication (jD.16p as 



This was justified because the commutation relations of the Y'^ resulting from (|D.28p 

[Y\Y^] {h \ j mn) = i^'^Y^ {h\jmn) 
are identical to those of (jD.2ip . As a result one directly identified 

Y^ \j mn) ^ \j m )^^^ 

However, comparing the operators obtained from (|D.28p on a representation of weight j we find: 

i A{j, n) 5„^ n-i 



(D.28) 

(D.29) 
(D.30) 



^m] = 

U J -1 mn 




1 

^2 


A(j,n + 1) (5„ „+i 


y2 ] _ 
(j)J mn 




i 
2 


A(j, n + l) 5yn „+i 


y3 ] _ 
{j)\mn 


2 L J ^ -^^J mn 







- A{j,n) n-i 

(D.31) 

One can already see an overall minus sign compared to (jB.6p . From this one can formulate ladder operators 



Y, 



^A{j,m+ 1) Srnn-l 
^A{j,m) Smn+1 



(D.32) 



but notices the difference in their definition as compared to (|B.3p . due to the assumed matrix multiplication (instead 
of the contragredient action (|D.6[) ) that 



^[Y(j)]mr 

r 


j 


r n ) 


— ~A{j, TO + 1) 1 J TO+1 n ) 




^[Y(^j)]mr 

r 




r n ) 


= -A(j,77l) j TO-1 n ) 




^[Y(])]mr 

r 


j 


r n ) 


~ —TO 1 J TO 71 ^ 








r n ) 


= .j(J + 1) j TO 71 ) 


(D.33) 



where (Y)^ = (Y^)^ + (Y^)^ + (Y^)^ = 1^^^+ + ^Y+Y- + (Y^)^. Upon defining the action of the F's in this way 
there is a correspondence Y^ ^ and therefore one realizes the following commutation relations (as compared to 
the commutation relations of the usual angular momentum operators acting on spin states | j m )). 

= -Y+ 



Y\Y^] =ie'^''Y'' 



(J)2, J'^l ==0 



Y^,Y+ 
Y^,Y- 
Y+,Y-'] 



= Y- 



-2Y^ 



= j+ 

= -j- 

= 2J3 



(D.34) 



Although the Y^s obey the usual angular momentum commutation relations on the left hand side of (|D.34|) . there is a 
difference in the signs of the ladder operator commutator algebra, providing possible phase convention obstacles when 
applying techniques of recoupling theory of angular momentum. 
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D.5.1 Recoupling Schemes 

In (jD.25[) we have found a matching between the basis states | j m ; n ) as defined in (jD.18[) to the states \ j rn ^ 
of an ordinary angular momentum system. Using (jD.18|l . where 



(") 



I j m ;n ) = (—1)"' ™\ j ~m n ) and ( j m ;n | = (—1)"' j —m n | (D.35) 
we can relate the states | j —m n ) to the abstract spin states | j m ) ^^^^ and can relate the recoupling scheme ()D.26|) 

I a J M ; n j ) := ^ ( ji mi ; ^2 ^^2 | 02 ) x . . . x ( apf-i niN^i ; Jn rriN \ J M ) 

X I ji mi ; ni ) (g) I j2 TO2 ; "2 ) ® • ■ • I JJV mjv ; riN ) 
■= X! ( ii TOi ; J2 TO2 I 02 TO2 ) X . . . X ( apf-i ; jN \ J M ) 

X (_l)Ji + ---+J"(_;L)-mi-...-m„ ^ | m ) (g) | ^2 - ^2 ^2 ) (g . . . | jjV -mjv riJV ) 



Introduce /ifc = —mki j^k ~ —"rhk, M = ^ 



= X! ( -^1 ' -M2 I 02 -Ai2 ) X ... X ( ttN-l -fj-N-l ] jN - ^J'N | J ) 

X + X I ji ^1 m ) (g I J2 M2 "2 ) «) • • • ® I JAT MiV J^TV ) 

(D.36) 

Using a symmetry property (See e.g. [^S], P-42) of the Clebsch-Gordan coefficients with respect to permutation of 
their argumentsHj 

( ii -mi ; 32 -1712 I 02 -r7i2 ) = (^-iyn-j2+a2 ^ . | 02 m2 ) (D.37) 

we can continue 

\d J M ; n j ) ^ ^ ( Ji Ml ; ^2 M2 I 02 /i2 ) X . . . X ( a^v-i Mtv-i ; jN tJ-N \ J fJ- ) 



Ml+M2 + ---+A'lv=A' 



X X I ji ^1 ni ) (g) I j2 M2 «2 ) «) • ■ • «) I JAT UN riN ) 

in complete analogy to (|D.35p . A similar calculation gives for the bra state 

{a J M ] n j \ = ^ ( ji Ml ; i2 M2 I 02 M2 ) X . . . X ( ajv-i Mjv-i ; Mw | M ) 

Ml+M2 + ---+A'lv=A' 

X X ( ji /ii m I g) ( J2 M2 «2 I g) . . . g) ( jiV Mw nN I 



(D.38) 



(D.39) 



We can thus conclude that for gauge invariant states J = Q, M = n — the recoupling schemes built from the 
preliminary basis \ j mn^ and the final basis | j m ; rt ) are identical and thus as long as one works with gauge invariant 
states the "historical" approach coincides with the construction implemented in the previous section. Moreover, since 
J + fj, = J — Misan integer number, the expectation values of operators diagonaO with respect to J, M will also 
be identical using the conventions (|D.11|1 . (jD.28|) with respect to our conventions ()D.18p . ()D.22|) . This is in particular 
the case for the volume operator. 



^The sum j\ + j2 — 02 is an integer number. 

'^By diagonal we mean that these operators require identical J, M upon evaluation in a recoupling scheme basis. 
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